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1 INTRODUCTION 


1 . 1 SYSTEM  BACKGROUND 

The  following  report  describes  a manpower  planning  system  that 
was  developed  to  assist  the  Army  in  improving  its  recruiting  programs. 

This  system  performs  three  basic  functions: 

1.  The  system  identifies  factors  that  significantly  influence 
the  number  of  supply  limited  accessions  that  the  Army  will  recruit, 

and  identifies  the  estimated  relationship  that  exists  between  these  factors 
and  the  estimated  accessions. 

2.  Given  restrictions  on  its  accession  resources,  the  system 
indicates  the  amount  of  resources  that  should  be  devoted  to  various 
budget  activities  (e.g.,  it  determines  the  amounts  of  money  that  should 
be  allocated  to  advertising,  recruiters  and  recruiter  assistants) . 

3.  Once  the  number  of  recruiters  and  canvassers  that  are  to  be 
made  available  during  the  year  are  estimated,  the  system  proceeds  to 
optimally  allocate  them,  by  DRC,  throughout  the  country. 

The  system  described  in  this  report  is  composed  of  three  types 
of  models.  Detailed  summaries  of  these  models  are  provided  in  the  next 
section  of  the  current  chapter,  and  detailed  descriptions  are  provided 
in  Chapters  2,  3 and  4.  The  current  section  is  devoted  to  a brief  over- 
view of  this  material. 

The  first  class  of  models  consist  of  five  nonlinear  optimal 
budget  allocation  models.  Choice  of  a model  will  depend  on  the  analyst's 
view  of  which  model's  assumptions  most  closely  describe  the  situation 
at  hand.  Subject  to  budgetary  constraints,  which  are  exogenously 
determined,  all  models  will  indicate  how  available  resources  should  be 
allocated  among  competing  programs,  to  yield  maximum  returns  on  investment. 
In  general,  the  budget  allocations  obtained  will  differ  with  the  model 
chosen. 

The  second  type  of  model  is  concerned  with  optimally  distributing 
recruiters  and  canvassers  throughout  the  United  States,  by  DRC.  This 
model  can  be  used  to  study  how  accessions  of  supply  limited  cohorts 
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(i.e.,  mental  category  I-IIIA,  high  school  graduate  males)  vary  with 
changing  economic  and  social  conditions,  and  with  changes  in  the  service's 
recruiting  activities.  The  number  of  recruiters  and  canvassers  avail- 
able for  allocation  will  be  determined  by  the  optimal  budget  allocation 
model  discussed  above. 

The  third  type  of  model  is  a stepwise  time  series  regression 
model  that  seeks  to  estimate  the  number  of  supply  limited  contracts 
that  the  service  can  expect,  given  past  and/or  present  values  for  such 
variables  as  the  number  of  recruiters  on  station,  the  national  unemploy- 
ment rate  (deseasonalized) , the  existence  of  enlistment  incentive 
programs,  and  the  existence  of  relevant  policy  constraints  on  service 
accessions.  This  model  is  useful  in  identifying  factors  that  significantly 
influence  yearly  accessions,  and  can  be  used  in  formulating  the  objective 
functions  df  the  optimal  recruiter  and  budget  allocation  models. 

1.2  SYSTEM  COMPONENTS 

Optimal  Budget  Allocation  Model  (OBAM) 

The  optimal  budget  allocation  model,  discussed  in  detail  in 
Chapter  2,  is  a nonlinear  optimization  model.  The  objective  of  the  model 
is  to  maximize  return  from  money  invested  subject  to  specified  budget 
constraints.  The  model  projects  accessions  on  a national  level  and  is 
constrained  by  limits  on  activities  such  as  recruiting  and  advertising. 

The  model  produces  reports  that  illustrate  the  expected  accession 
mix  of  high  school  and  non-high  school  graduates  and  mental  group  lVs. 

High  school  graduates  in  mental  groups  I-III  are  considered  explicitly 
by  the  model.  Non-high  school  graduate  and  mental  group  IV  accessions 
are  computed  on  the  assumption  that  they  are,  as  a matter  of  policy, 
demand  limited. 

This  OBAM  was  developed  in  the  course  of  other  studies  on 
recruiting.  It  was  determined  that  this  was  a prime  candidate  for 
inclusion  in  the  Market  Analysis  Subsystem.  Since  this  was  not  part  of 
the  study  entitled  Analysis  for  Management  of  Recruiting  Resources  and 
and  Operations,  it  is  possible  that  the  model  may  not  be  as  refined  as 
required  for  certain  recruiting  analyses. 
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Optimal  Recruiter  Allocation  Model  (ORAM) 

The  optimal  recruiter  allocation  model,  discussed  in  detail  in 
Chapter  3,  is  of  the  nonlinear  optimization  type.  Given  national,  regional, 
combat  arms  and  term-of-service  goals  for  the  accession  groups,  this  model 
searches  for  minimum  allocation  of  recruiters  and  canvassers  throughout 
the  country.  A full  range  of  optimal  solutions  are  produced  by  the  model, 
and  these  enable  it's  users  to  investigate  the  effects  that  factors,  such 
as  future  unemployment  rates,  have  on  their  planning. 

Output  produced  by  ORAM  consists  of  report  evaluations  which 
are  based  on  the  optimal  recruiter  allocations  that  it  generates.  These 
evaluations  are  conducted  on  the  DRC  level,  and  provide  the  following 
information: 

• DRC  identification; 

• DRC  number; 

• Number  of  recruiters  on  station; 

• "amber  of  canvassers  (where  applicable); 

• Number  of  accessions  (total); 

• Number  of  accessions  by  cohort  groups 

In  addition,  subtotals  by  region,  national  totals,  regional  percentages  and 
national  summaries  are  also  provided. 

MERGE/BMD  Time  Series  Region  Model 

The  MERGE  component  of  the  MERGE/BMD  system  has  the  task  of  handl- 
ing data;  the  BMD  portion  of  this  system  consists  of  a stepwise  regression 
model  for  processing  MERGE's  output.  By  combining  MERGE  with  BMD,  large 
scale  time  series  regression  problems  can  be  formulated  with  a minimum 
of  program  preparation.  Objectives  for  the  time  series  regression  model 
are  twofold: 

• To  examine  recruiting  factors  which  influence  service 
enlistments,  based  on  a month-by-month  analysis  of  enlistments. 

0 To  forecast  enlistments  on  a monthly  basis  into  the 

future. 

Model  output  is  on  a regional,  rather  than  national  level  and  can  involve 
either  numbers  of  accessions  or  enlistment  contracts. 


( 
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2 BUDGET  ALLOCATION  MODEL 


2. 1 INTRODUCTION 

A major  concern  in  any  market  is  the  investment  required  to 
stimulate  the  "purchase"  of  a product.  The  U.S.  Army  is  concerned  with 
the  potential  demand  for  Army  careers  that  may  be  expected  from  its 
investment  in  recruiting  and  advertising  programs.  The  question  that 
one  must  be  able  to  answer  is  the  allocation  of  budget  dollars  in 
various  programs  that  maximize  the  return  for  each  dollar  invested. 

This  chapter  considers  five  possible  models,  and  the  methodology  for 
each  form.  One  model  (multiplicative  exponential)  has  been  included 
in  the  market  information  subsystem. 

2.2  METHODOLOGY  FOR  OPTIMAL  BUDGET  ALLOCATION  MODEL  (OBAM) 

Several  models,  all  of  which  have  decreasing  rates  of  return 

for  large  budgets,  were  considered  as  potential  candidates  for  part  of 
an  optimization  study.  The  general  character  of  the  optimal  returns  per 
total  budget  is  represented  by  Fig.  2.1. 


Optimal 

return 


Total  budget 


Fig.  2,1 — Variations  in  Optimal  Return  with  Changes  in  the 
Bud get 


All  curves  appearing  in  Figs.  2.1  through  2.4  are  optimal  in  the 
sense  that  the  budget  constrained  problem  (P)  is  solved.  The  statement 
of  this  problem  is  as  follows: 


Problem  P: 


maximize  F(X. ,...,X  ) 

1 P 

Xi Xp 

Subject  to  constraints 


(return  function) 


P 

Z X.  - 3 (total  budget  constraint) 
i-1  i 

X > 0 i » l,...,p  (no  negative  budgets) 

where 

X^  is  the  budget  allocation  to  program  i and  p is  the  total 
number  of  programs. 

The  Lagrangian  conditions  for  a constrained  maximum  of  ? are: 

3L  ? t 

r - = 0 (i™l,...,p)  where  L = F+  X(B-  IX. ) 

1 1=1 

provided  the  total  budget  constraint  is  binding  on  the  optimal  solution. 


Five  different  models  for  F have  been  considered  and  the  associated 
optimal  solutions  have  been  written  down.  For  reference,  F is  indexed 
by  ^ according  to  the  following  convention: 
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Models  F 

n 

Model 

Abbreviation 

Name  of  Model 

Mathematical  Form 

Fi 

CD+ 

Cobb-Douglas 

Additive 

5 3i 

F = Z A.  X. 

1 i-1  1 L 

F2 

CDX 

Cobb-Douglas 

Multiplicative 

P 3 

r,  - a i x. 

2 i-i  1 

V 

‘3 

EXP+ 

Additive 

exponential 

F « I Ai(l-e_3iXi) 

i*l  1 

‘4 

EXPX 

Multiplicative 

exponential 

? _ -3  y 

? = A TT  (1-e  i i) 

4 i=l 

F5 

EXPXS 

Square  root 

multiplicative 

exponential 

P -B. 

F = A n 1 - (1+B  1 

3 i=l  11 

Remark : The  expression  in  the  parentheses  was  used  as  part  of  an 

objective  function  in  reference  (1)  and  is  inserted  here. 

Two  of  the  models  CD+  and  EXP+  are  additive,  i.e.,  returns  from 
each  program  i are  considered  separately  and  then  together  to  form  the 
return.  The  remaining  three  models  are  multiplicative  and  represent  the 
type  of  return  which  depends  on  a contribution  from  every  supporting 
program. 


Parameters 

The  parameters  A^,...,A^;  are  estimated  in  this  paper 

foi  purposes  of  illustration  using  point  values  for  returns,  budgets, 
and  marginal  productivities.  The  minimum  amount  of  information  needed  to 
"fit"  the  models  (i.e.,  determine  the  parameter  values)  is  3p  for  the 
additive  models  and  2p+l  for  the  multiplicative  models. 


Determining  the  Model  Parameter  Values  (using  minimum  information) 
Additive  models 


F = E A.f,  where  f.  depends  only  on  X.  and  the  value 
j=l  j J J J 

of  the  parameter  B. 


Input  information 


base  point  budget 


X. 

1 P 


Base  point  marginal  productivities  v,,...,v  ( 3 values 

1 P ' p 


base  point  return  by  program 


al»* •' ,ap 


Solve  the  following  2^  equations  for  A^  B^  (the  equations  are  in 
in  general  transcendental). 


a_F  _ A.  a£i 
ax.  " 1 ax.  vi 

i i 


A.fj  = a.  provided  X.  = X. 

ill  li 


combining  the  above  equations  yields  an  equation  to  be  solved  for  B. 


a f i f i / vi  \ 

- 1 ] = ° solve  for  B.  (i=l,...,I 


followed  by  the  equation  in  A. 


A.  = i where  f.  is  evaluated  at  X.,  B. 
iy  i i’  i 

i 


nple  - Additive  exponential  model 
Let  F = Aj  (l  - e + A2  (l  - 


-B2X2 


Solve  for  B^,  A^  the  equations 
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-B  X.  -b,x.\  /„  \ 

B e 1 - ( 1 * e ) / _1  j=  0 

\ } v ai  / 


i = 1,2 


-bitm 


Mul t iplioative  Models 


F = A t r f.  where  f.  depends  only  on  X and  the  value 

-is.  i J j ^ 

J of  the  parameter  B. 


Input  Information 


base  point  budget 


X,,...,X 
1 P 


base  point  marginal  productivities  v^,...,v 


base  point  return 


Solve  the  following  p+1  equations  for  A,  B^...,B 


dF  F _ where  X.  = X 

ax.  ~ F?  ax,  “ Vi 

l i i 


equivalently  since  a = A u f.  (X.) 

j = l J J 


a£i  Vi 


A = a 


i = 1 > • • • > P 


where  f.  is  evaluated  at  X B 
l i i 


Example  - Cobb-Douglas  Multiplicative  Model 


B,  B 


B1^B2 


Let  F 


A X1  1 X2  2 i.e.  = xj  ^ 311(1  3 = A Xi  X; 
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Solve  for  B.,  A the  equations 


since  ^ 0 

Bi  = Vi*i  i = 1,2 
a 


and  thus 
A = a 


Note  in  this  example  the  equations  in  are  in  closed 
solution  can  be  found. 


Solving  for  the  Lagrangian  Stationary  Conditions 

Solve  qJC  = 0 where  = F + X MB  - L 

' i=l 

Additive  Models 


dc 


Example 
F ■ 


Additive  Exponential  Model 


form  and  a direct 


(Lagrangian  stationary 
conditions) 
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A.B.e 
1 i 


-B.X. 
1 i 


X = _±_  2n 

1 B.  A.B. 

i li 


Multiplicative  Models 


ax.  f.  X. 

l l i 


or  equivalently 


Provided  A 


=£  0 i.e. , X.  s 0 


3f. 

-T7T-  - k f = 0 

o X^  l 


Example  - Square  Root  Multiplicative  Exponential  Model 


F = A n 

J- 


2 f -B  v/xA 

T ^1  - (1+B.  v/Xj)  e 3 3 J 


7 -kfi-  i *Biv/Fi  - + 


Note  that  in  this  model  the  solution  X^  for  specified  k is  not  explicit. 
Recursive  methods  are  necessary  to  solve  the  equations. 
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Sketch  of  Optimal  Solutions  to  the  Five  Models 

A brief  sketch  of  the  optimal  solutions  to  problem  P for  each  model 
follows: 

Optimal  Solutions  to  Models 

Optimal  values  of  the  models  are  denoted  by  F 

Optimal  solutions  of  the  models  are  denoted  by  X*  i-l,...,p 

Optimal  productivities  of  the  models  are  denoted  by  X* 

The  change  in  optimal  program  solutions  with  total  budget  is  denoted 
i<. 

by  dX. 

i 

dB  * 

Cobb  Douglas  Additive  Model  F 


X \3.-l 

Vi  1 


i . e . , X = A . 3 . X . 

ill 


i*l» • • • ,p 


P * B 
Z A.  X. 
i-1  1 1 


P * 
Z X. 
i-1  1 


* 

dX  — - 

i _ dX_ 

dB  “ dB_ 
dX 


Cobb-Douglas  Multiplicative  Model  F, 


* Bi 
X.  - — 
i k 


where 


* * 
X = k F_ 


A it  B i 
* i-1  1 


it  k i 
i-1 


(1*1 p) 


Z B. 
i-1  1 


dX . B . 

i i 


Additive  Exponential  Model  F. 


*11  1 
X.  - - — In  A.B.  - — — — In  X 

1 i Bi  1 1 ei 


* r r -I 

Ft  - E A.  - X E 
3 1 i-i  3i 


if  X < A.  3. 

i i 

otherwise 


X = A. B . e 

l l 


•k 

-B.X. 

i l 


• 1 v \ 
3 = E — In  A 3.  - In  X [ E ~ 

i-i  3i  1 1 In  3i 


i=l i • • • > p 

* 1 

dX.  3 

i _ i 

dB  “ E 1_ 

iel  3. 


Multiplicative  Exponential  Model  F, 

Zl 


= bt  ln  (Lir)  whare 


A » k F. 

4 


* 4 

F = 

AK 

1=1  \ i 


i=*i 


ki+B . 
i 


i-1  • i 


{i|x*  >0} 
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Operational  Modes 
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Square  Root  Multiplicative  Exponential  Model  F 


Solve  for  X. 

1 

/ \ —B . X. 

k ( 1+Bi  \ J + e 1 1 - k = 0 


★ * 
where  \ = k F<. 


. * _ A /■  p D\  2 i=l 

mnZiV  e 


2 B.  vx. 

■ 1 1 


dX  is  not  available  in  explicit  form. 
dB 

APPLICATION 

An  exercise  which  compares  the  four  models  using  a base  point 
described  by  reference  (2),  pages  54  and  55,  is  presented  here.  T 
input  information  is  as  follows:  (Tables  2.1  and  2.2). 


Table  2.1 
Input  Values 


BASE  POINT 


Program 


Recruiter  Assistants 
and  Canvassers 

Print  Media 

General  Recruiting 


Accessions 

Productivity 

Budget 

5C000* 

700 

20.4 

50000* 

710 

13.0 

50000* 

210 

90 

150,000  was  used  for  the  multiplicative  models 


V'.  ,'Kf 


The  parameter  values  defined  for  the  four  models  are  as  follows: 

Table  2.2 
Parameter  Values 


\ Model 

Paraae  cers  NT ^ 

Values  | \ 

Y \ 

CD+ 

COX 

EXPX  | 

EXPXS 

Ai 

21131.95264 

— 

56676.906 

— 

— 

31141.25146 

— 

53327.027 

— 

— 

9125.703117 

— 

61057.220 

— 

-- 

A 

— 

54529.97061 

— 

162061.6 

172456.7643 

1 

31 

.2356 

.0952 

.10433833 

.130407655 

1.078554806 

32 

.1346 

.06153333 

.21341253 

.326369351 

1.533629790 

33 

.373 

.125 

.013985365 

' .0367115143 

.465569754 

Table  2.3 

Model  Solutions  at  3udgec  Level  123.1  — $M 


Accessions  ? 


3ase  Point 
Model  I 

T 

20.4 

13.0 

90.0 

123.4 

150,000 

CC+ 

53.39 

30.73 

39.22 

123.36 

160,952 

CDX 

41.57 

26.37 

55.02 

123.46 

157,760 

EXP+ 

29.43 

17.50 

76,46 

123.40 

152,901 

EXPX 

25.32 

15. 73 

32.17 

123.23 

151,263 

EXPXS 

28.16 

12.44 

77.78 

123.38 

152,621 

Table  2.4 


MODEL  DIFFERENCES  FROM  BASE  POINT  ALLOCATIONS 


^2 

AX3 

^Accessions 

Model 

CD+ 

32.99 

17.73 

-50.78 

10,952 

CDX 

21.17 

13.87 

-34.98 

7,760 

EXP+ 

9.03 

4.50 

-13.54 

2,901 

EXPX 

4.92 

2.78 

- 7.83 

1,268 

EXPXS 

7.76 

4.44 

-12.22 

2,621 

RESULTS 

The 

OF  MODEL  COMPARISON 
optimal  solutions  of 

the  base  point 

appear  in 

Table  2.3.  Table  2.4 

gives  the  optimal  readjustment  to  program  budget  allocations  at  the  base 
point  budget  of  123.4  million  dollars. 

CONCLUSIONS 

Considerably  different  program  budget  allocations  result  from  using 
different  models.  Except  for  the  Cobb-Douglas  additive  model,  all  models 
allocated  the  most  funds  to  general  recruiting  (program  3)  followed  by 
recruiter  assistants  and  canvassers  (program  1)  and  print  media  (program 
2).  The  Cobb-Douglas  additive  model,  however,  allocated  more  funds  to 
program  1 than  program  3. 

The  types  of  optimal  solutions  that  each  of  the  models  produces  should, 
perhaps,  be  considered  as  a factor  in  the  selection  of  the  particular  model 
to  be  used  in  a given  situation.  For  example,  the  Cobb-Douglas  multipli- 
cative model  yields  program  solutions  that  are  always  a fixed  percent  of 
the  total  budget  allocated,  no  matter  what  the  total  budget  is.  The  Cobb- 
Douglas  solutions  were  more  clustered  than  the  exponential  solutions  in  the 
exercise  here.  Table  2.5  gives  the  limiting  percent  of  program  budgets  as 
the  total  budget  increases.  The  multiplicative  exponential  had  the  greatest 
divergence  in  program  percent  allocations. 
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Ic  should  be  noted  chac  the  base  point  allocation,  although  not 
optimal,  was  least  displaced  by  the  exponential  type  models  (Table  2.4). 


LIMITING  VALUES  FOR  OPTIMAL  PROGRAM  ALLOCATION 


As  the  tocal  budget  B is  made  large,  the  optimal  budget  allocations 
* 

to  the  programs  approach  limiting  fixed  percentages  of  the  tocal 
budget.  The  percentages  differ  for  each  model  and  are  functions  of  Che 
parameter  values.  The  limiting  percentages  are  as  follows: 


B. 


CDX 


100 


P 

Z 3. 
i-1  1 


1-3. 


CD+  X. 

1 


100 


(BiV 


E (BiAi}  nr 

i=i  1 1 1 3i 


"CP+  and  EM?” 


B, 


100 


\ r 


-2 


EXPXS  X. 

l 


100 


R 

i 

i=l 


-2 


For  the  values  of  the  parameters  of  the  five  models  the  limiting 
percent  allocation  of  the  tocal  budget  is  shown  in  Table  2.5. 


Table  2.5 

Percent  of  Total  Budget  Allocation 


a 

* 

Model 

X2 

X3 

CD+ 

27.06 

5.63 

67.29 

CDX 

33.67 

21.76 

44.56 

EXPH- 

15.46 

8.53 

75.99 

EX?X 

9 . 22_ 

3.32 

87.45 

EXPXS 

14.58 

7.16 

78.25 

2.3  DYNAMIC  PROGRAMMING  APPROACH  TO  OPTIMAL  BUDGET  ALLOCATION 

Throughout  the  discussion  in  this  paper  focus  has  been  directed 
to  optimal  allocations  of  budgets  to  programs.  If  an  optimal  allocation 
is  greatly  displaced  from  a base  point  allocation  then  a reallocation 
of  program  budgets  may  "jar"  the  system.  In  such  cases  a gradual  read- 
justment of  total  budget  may  be  made  such  that  restoration  to  optimal 
behavior  proceeds  along  reasonably  smooth  lines.  Reference  (3)  addresses 
this  question  for  the  case  of  the  multiplicative  exponential  model. 


2.4  INPUT  REQUIREMENTS 

The  input  requirement  for  the  budget  vary  according  to  the  number 
of  programs  to  be  considered  in  any  run  of  the  model.  The  information  is 
required  to  define  basepoint  budget  considerations  and  is  discussed  below. 

Input  Data  Description 

The  first  two  input  cards  are  a descriptive  title  with  a maximum 
length  of  130  characters,  80  on  card  1 and  50  on  card  2.  Both  cards  are 
required,  even  when  the  title  is  less  than  80  characters  in  length.  The 
remaining  cards  depend  upon  the  type  and  number  of  budget  programs  to  be 
considered  with  a minimum  of  one  and  a maximum  of  five. 

The  first  six  columns  of  each  succeeding  card  is  a control  parameter 
that  may  have  any  of  the  following  values: 

BUDGET  - indicates  baseline  budget 
ACCESS  - indicates  baseline  accession 
BUDLOW  - minimum  budget  to  be  considered 
BUDUP  - maximum  budget  to  be  considered 
BUDINC  - budget  increments 

blank  of  program  name  - these  cards  are  required 

for  entering  productivity  values  for  each  program 
under  consideration 

END  - end  of  budget  definition 

In  the  case  of  program  productivity,  two  cards  are  required, 
identified  by  the  following  control  words  in  columns  11  uhrough  16. 
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PRRATE  - productivity  rate  for  budget  programs 
BUD-$M  - basepoint  budget  for  this  program  is 
millions  of  dollars 

A maximum  of  five  program  budgets  may  be  entered.  The  remaining  entry  on 
each  card  is  the  numerical  value  corresponding  to  the  identifying  control 
words.  These  values  are  in  columns  21-30  and  two  considerations  must  be 
considered.  First,  the  decimal  point  must  be  in  the  field  and,  second, 
all  budgets  are  expressed  in  terms  of  millions  of  dollars. 

Two  cards  follow  the  'END'  card  of  the  budget  definition  series. 
These  cards  provide  data  concerning  demand  limited  supply  groups  under 
the  following  formats. 


Card 

Columns 

Format 

Definition 

1 

1-10 

A10 

Service  name  (ARMY,  NAVY,  USAF,  USMC) 

11-20 

F10. 1 

Maximum  number  of  female  enlistments 

21-30 

F10.1 

Maximum  number  prior  service  enlistments 

31-35 

F10.5 

Expected  proportion  GED  enlistments 

36-40 

F10.5 

Expected  proportion  HSG  enlistments 

41-50 

F10.5 

Expected  proportion  MC  IV  enlistments 

2 

1-10 

A10 

Service  name 

11-20 

F10.1 

Maximum  number  of  MC  IV  HSG  enlistments 

21-30 

F10.1 

Maximum  number  of  NHSG  enlistments 

31-40 

F10.1 

Total  accessions  objective 

If  the  demand- limited 

reports  are 

not  desired  the  service  name  on  each 

must 

be  replaced  with 

the  keyword 

NO  in  columns  1-2. 

3 RECRUITER  ALLOCATION  MODEL 


3.1  INTRODUCTION 

This  chapter  describes  the  methodology  developed  for  determining  the 
optimal  allocation  of  recruiters  and/or  canvassers  throughout  the  United 
States,  and  the  resultant  allocations  by  DRC.  In  addition,  a complete 
description  of  input  data  requirements  is  provided  along  with  operational 
considerations. 

The  model  optimally  allocates  recruiters  and  recruiter  assistants  to 
stations  distributed  across  the  United  States.  Specifically,  several  classes 
of  supply-limited  groups  of  first-term  enlistees  are  defined,  e.g.,  high 
school  graduates.  Mental  Category  I-IIIA.  These  supply-limited  classes  then 
become  the  focus  for  the  minimal  cost  allocation  of  recruiters  and  recruiter 
assistants  under  a broad  range  of  objectives  which  include  national  goals, 
regional  goals,  combat  arms  goals,  and  term-of-service  goals. 

Model  input  parameters  which  are  evaluated  using  cross-sectional 
analyses  include  unemployment  rates,  regional  attitudes  toward  the  military, 
historical  propensities  to  enlist,  advertising,  and  other  influencing  factors. 
The  parameter  elasticities  of  the  recruiters  and  recruiter  assistants  for 
specific  types  of  supply-limited  groups  are  also  evaluated. 

A full  range  of  optimal  solutions  is  produced  by  the  model,  which 
enables  the  user  to  interact  and  ask  questions  about  the  effects  of  factors 
such  as  future  unemployment  rates  on  planning.  Sensitivity  analysis  is  also 
available,  again  pointing  out  the  effects  of  special  requirements  on  total 
minimal  cost. 

A system  overview  of  the  general  characteristics  of  the  inputs  and 
outputs  to  the  model  is  illustrated  in  Figure  3.1.  A more  detailed  discus- 
sion of  these  elements  will  be  presented  later  in  this  chapter. 
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3.2  THE  OPTIMAL  RECRUITER  ALLOCATION  MODEL  (ORAM)  METHODOLOGY 

This  section  contains  a precise  formulation  of  the  mathematical  opti- 
mization problems  which  are  solved  by  the  computer  algorithm.  A statement 
of  the  problem  is  first  presented, followed  by  the  Lagrangian  formulation 
of  the  equivalent  problems.  The  Lagrangian  stationary  conditions  are  then 
derived. 

Five  optimization  problems  (denoted  by  the  index  p;  p = 1,...,5)  will 
now  be  formulated.  The  notational  conventions  used  in  this  discussion  are 
as  follows: 


Notation 

The  following  notation  will  be  used  for  the  optimization  problems. 
Problem  size 

• Number  of  DRCs  N N s 64 

• Number  o-f  regions  M M = 5 (currently) 


Number  of  different 
supply  limited 
cohort  types 

J 

J <:  8 

Number  of  choices 
of  term  of 
service 

K 

K = 2 

Number  of  choices 
of  occupational 
specialty 

L 

Ls  2 

Number  of  independent 
variables 

N 

If  recruiters 
variables 

are 

only  independent 

2N 

If  recruiters  and  canvassers  are 
independent  variables 

Number  of  dependent 
variables 

NJ 

NJ<  512 

Independent  variables 

X, 


‘i 

: the 

i 

: the 

* 

: the 

* 

ip 

: The 

district  i 
district  i 

problem  p 
problem  p 


(i=  1 N) 

(i=l , . . . ,N) 
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Historical  percentages  or  fractions 


Cijk2 


the  fraction  of  total  enlistees  within  district  i of  co- 
hort type  j with  term  of  service  k and  occupational 
specialty  2 

(1-1.....N),  (j=l,...,8),  (k=l  ,2)  , (2  = 1,2) 
where 


H 

n 

* 

: 3 year  term  of 

service 

2 

: 4+  year  term  of 

service 

2=1 

: combat  arms 

2 

: non-combat  arms 

t 


ljk2 


ij  ki 


~ijk2 


J 

E 

j=l 


'ijk2 


Note : 


tijk2  = 1 
jk2  J 


Costs 

y : annual  cost  of  placing  an  additional  recruiter  in  a district 

6 : annual  cost  of  placing  an  additional  canvasser  in  a district 


Upper  limits 

Y : upper  limit  of  total  number  of  canvassers 

Primary  accession  requirements 

Tg  : total  accessions  across  supply  limited  cohort  types  j, 

(j=l J) 

Tg  is  accession  requirement  for  optimization  problem  1 
(p=l) 

: term  of  service  requirements  across  supply  limited 

cohort  types  j,  (j=l,...,J) 

k=l,2  corresponds  to  optimization  problems  2 and  3 

(P«2,3) 
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<p  T 

s 


occupational  specialty  requirement  across  supply  limited 
cohort  types  j, 

i=l  corresponds  to  optimization  problem  4 (p=4) 

total  accessions  across'  supply  limited  cohort  types  j, 

in  region  indexed  by  m,  (m=l,...,M)  where 

< o is  the  fraction  of  T which  is  required  to  be  in 
^m  s 

region  m. 

Note:  generally  £ <pm  = 1 

m 

The  set  of  regional  requirements  corresponds  to  optimization 
problem  5 (p=5) 


Dependent  variables 


number  of  accessions  of  cohort  type  j enlisting  in  district  i 


Qalj  wtt2j  “3J  Ga4j  a5j 
Qij  ij  ij  Gij  ij 


a 


(x, 


6j 


a 

+ Y. 


7j 


Statement  of  the  five  optimization  problems 


Problem  1 


min 


Y 


i 


N 


£ (yXi  + 6Y  ) 
i=l  1 


(cost) 


k 

The  terms  , A , are  defined  in  Chapter  5,  Analysis 

for  Management  of  Recruiting  Resources  and  Operations,  Paul  S.  Souder,  Jr., 
et  al,  OAD-CR-166. 
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subject  to 


N»J  a,.  a~.  a.,  a.,  a,..  /cl,  a,,\ 

£ Q w 2j  u.?J  g.4j  a 5j  (x.6:i  + Y.7j  )v  = 

i=l  ij  ij  ij  ij  ij  v 1 1 / j 


(total  accession  requirement) 


S Y £ Y (upper  limit  on  canvassers) 

i-1  1 


Problems  2,  3 


min  2J  <yx  +6Y.) 

X1,Yi  i=l  1 


(cost) 


subject  to 


NSJ  ( s t'  Vu  «y  y «4j  «5j  ^ a6j  + «73\ 
i-1  \i=l  iJk7  ij  ij  ij  ij  \Xi  + Yi  > j 


£ Y £ Y 
i-1  1 


(term  of  service  requirements,  k = 1,2) 


(upper  limit  on  canvassers) 


where  k = 1 corresponds  to  three-year  term  of  service 
requirement  (problem  2)  and  k = 2 corresponds  to  four- 
or-more-year  term  of  service  requirement  (problem  3) 


Problem  4 


min  £ O^X  + 6Y  ) 
Xi,Yi  i-1  1 


(cost) 


subject  to 


k 

£ t 


a,,  a. 


a, , at 


g v-  iiuJ " ;ij  v ^ 


(occupational  specialty  requirement) 
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I 


L Y < Y 
i=  1 1 


(upper  limit  on  canvassers) 


where  l = 1 corresponds  to  combat  arms 


Problem  5 


min  L OX  + 6Y  J 
Xi,Yi  i-1 


(cost) 


subject  to 

N,J  a,.  a~.  a,,  a,,  a, , / a,  a7,\ 

ij  ij  ij  ij  ij  V i i / J m s 

1 (m=l,...,M) 

(regional  accession  requirements) 

S Y.soY  (m=l,...,M) 

J6I  1 m 

in 

(regional  upper  limits  on  canvassers) 
where  I = { i | district  i is  in  region  m ] 


I 


The  Lagrangian  Formulation  of  the  Optimization  Problems 

Several  notational  simplifications  will  be  made  before  proceeding 
to  the  Lagrangian  formulations  of  the  five  optimization  problems.  The 
simplifications  are  as  follows: 


a.,*  a'>A  «•»-!  a,  . a,, 
c = o ^ w ^3  v ^3  r ^ a 
cu  qij  "u  vij  Gij  Au  vj 

bJ  ' “6J 
d3  ma7J 


'll 


(A  ‘ «“) 


iJP 


( & C'i]2t)  °i 


( 


t-1  tlJkl) 


'ij 


if  p = 1 
if  p = 2 

if  p = 3 

if  p = 4 
if  p = 5 


The  five  optimization  problems  now  are  rewritten  using  the  concise 
notation. 

Problem  p,  (p=l,...,4) 

N 

min  £ 0%  + 6Y.)  (cost) 


j 


(primary  requirement  constraint) 
N 

E Y.  s Y 
i-1  x 


N J 

subject  to  L E h 
i-1  j=l 


ijP 


Cxi 


if  p = 1 
if  p = 2 
if  p = 3 
if  p » 4 


(canvasser  constraint) 
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Problem  5 (p=5) 


min  £ (yY  + 6Y  ) 
Xi,Yi  i-1 


subject  to: 


J b d 

s s h X^+Y^  =cpmTs  (m=l M) 


i£i  j=i 

in 


(primary  requirement  constraints) 


s Yi  * <PnY 

i£I 

m 

(canvasser  constraints) 

Lagrange  Multipliers 

Lagrange  multipliers  corresponding  to  the  constraints  of  problem 
p (p=l,...,5)  are  denoted  as  follows: 

For  p=l, ... ,4 

X : the  Lagrange  multiplier  corresponding  to  the 

P 

primary  requirements  constraint  of  problem  p, 

(p=l, . • • ,4)  Xp  k 0 

u : the  Lagrange  multiplier  corresponding  to  the 

P 


For  p=5 


canvasser  constraint  of  problem  p,  (p=l,...,4) 


u 2 0 
pP 


the  Lagrange  multiplier  corresponding  to  the 

five  regional  accession  requirements  of  problem 

5.  X S 0 (m=l .... ,M) 
m 

the  Lagrange  multiplier  corresponding  to  the 

five  regional  canvasser  requirements  of  problem 

5.  a 0 (m=l,...,M). 

m 


Lagrangian  Functions 

* 

Five  Lagrangian  functions  are  now  formulated. 

For  problems  (p=l,...,4)  the  Lagrangian  function  is  defined  by  the 
equation: 


For  problem  5 (p=5)  the  Lagrangian  function  is  defined  by  the 
equation 


» « _ / J ^ bi  dl\ 

L (yx  +6Y.)  - L A f L h ' X J + Y J ) - 

i i m=l  m ( i€  I ijpV  1 1 J 

\j=lm 


pT 

s 

m 


The  Kuhn-Tucker  Conditions  for  the  ORAM 

The  five  Lagrangian  functions,^ , (p=l,...,5)  are  convex  in  X^,  Y^^ 

provided  0 < b^  £ 1,  0 < d s 1,  (j=l,...,J).  Thus  the  Kuhn-Tucker 

necessary  conditions  are  also  sufficient  for  an  optimal  solution. 

The  Kuhn-Tucker  conditions  that  are  satisfied  by  optimal  solutions 
* * 

Xip,  Y , solved  by  the  ORAM  are  for  p=l,...,4  as  follows: 


The  following  references  are  fundamental  to  the  study  presented  here: 

• Kuhn,  H.W.,  and  Tucker,  A.W.  , "Nonlinear  Programming,"  Proceed- 
ings of  the  Second  Berkeley  Symposium  on  Mathematical  Statistics  and 
Probability. " University  of  California  Press,  Berkeley,  1951. 

• McCormick,  G.P.,  "Second  Order  Conditions  for  Constrained  Minima," 
SIAM  Journal  Applied  Math.  15  (3),  pp.  641-652,  1967. 

• Fiacco,  A.V. , and  McCormick,  G.P.,  Nonlinear  Programming,  Sequential 
Unconstrained  Minimization  Techniques,  New  York,  John  Wiley  and  Sons,  Inc., 
1968. 
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I 


a=f. 


£_ l 


X* 

ip 


o,  a-^ 


Y* 

ip 


0 (i=  1 , . . . ,N) 


d*T 


x* 

p 


°’ 

p 


<;  0 


(p=l,. . . ,4) 


u*  (y  - Z Y*  ) = 0 
P \ i=l  iP / 


X*p  a 0,  Y*p  2 0 (i-l,...,N) 


X*  i 0 , n*  s 0 

where  X*.  u*  are  the  values  of  the  Lagrange  multipliers  at  the  optimizing 
P P 

point  X*  , Y*  . 
ip  IP 

For  problem  5 (p=5)  the  Kuhn-Tucker  conditions  are 

3X! p 


IX* 


m 


- 0 

*p  . 

o 

(i=lj • 

• • jN) 

- u, 

- 0 

o 

(m=l, . 

• • >M) 

- u J 

V 

K (v  ■ 

z yO= 

161  P 
m 

0 

(m=l, . 

• • »M) 

xip  2 °* 

y*pso 

bV 

IV 

o 

(m=l, . 

. . jM) 

where  X*,  u*  (m=l,...,M)  are  the  values  of  the  Lagrange  multipliers 
m m 

at  the  optimizing  point  X*  , Y*  . Points  satisfying  the  Kuhn-Tucker 

ip  ip 

conditions  are  called  Kuhn-Tucker  points. 
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A Procedure  for  Solving  for  Kuhn-Tucker  Points 

This  section  is  very  detailed  and  describes  the  algorithm  which  is 
used  to  solve  the  optimization  problems  by  means  of  a process  which  con- 
verges to  the  Kuhn-Tucker  points  just  described.  The  algorithm  seeks 
solutions  to  the  following  set  of  equations  and  inequalities: 

For  problems  p,  (p=l , . . . ,4) 


J h 1 V J V1 

E b h X°j_1  , E d.h,  . Y . 6 + M 

j = l 3 1 = Xp  j=l  J 


(i=l, . . . ,N) 


N J A d \ 

E E h.  X.J  + Y.J  ) = R . Y - E Y.  Jfc  0 
i=l  j = i iJP  V 1 i/  P i-J  i 


/y  - E Y V 0 

\ i=  1 V 


X.  a 0,  Y.  s 0,  X s 0,  n s 0 (i=l N) . 

i l p p 


For  problem  5,  (p=5)  the  following  set  of  equations  and  inequalities  is 
solved. 


J b -1 

ViiPXi 


J d -1 

y s d.h..  y,j 

j 


6 


m 


E E h X j = - E Y4  * 0 

i^I  j-1  X3P  1 


m 


m S’  vm*  " ‘i 
m 


^m  W ' S Yi 

v i£I 

m 


X > 0,  Y > 0,  X ^ 0,  M a 0 
i l m m 


where  i£I  , (m=l,...,M) 
m 
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- A three  part  decision  rule  is  used  to  solve  for  Kuhn- Tucker  points. 
Each  part  is  called  a branch  and  uses  a Fibonacci  search  technique  to 
determine  the  solutions  to  its  particular  set  of  equations.  A Newton 
method  is  used  to  solve  equations  (la),  (2a)  (3a)  (following)  within 
branches  1,  2,  3 as  functions  of  the  right  side  values.  The  Fibonacci 
search  is  made  to  satisfy  equations  (lb),  (2b),  and  (3b). 

The  decision  branch  rules  for  problem  p,  (p=l,...,4)  will  now  be 
discussed. 

Decision  Branch  Rules  (problem  p,  p=l,...,4) 


« Branch  1.  Solve  for  Y (i-1, . . . ,N)  and  v the  (N+l)  equations 
J d.-l 

(la)  2 d.h  Y 3 = v (i=l,...,N) 

J j-jp  1 


where  v is  a candidate  value  for 


such  that 

N 

(lb)  2 Y.  = Y 
i=l  1 


If 


N J 
2 2 
i-1  J=1 


‘ijp‘i 


J „ 


N J < 

S E h.,  Y. 
i-1  j-1  ljp  1 


= Z <;  R 


j . 


= Z > R 


then  retain  the  value  v and 
go  to  branch  2 


go  to  branch  3 


• Branch  2.  Solve  for  X^(i=l , . . . ,N)  and  w the  (N+l)  equations 


(2a) 


J 

E 

J-l 


b.h,  . 

J iJP 


j 

i 


such  that 


(2b) 


N 

E 

i-1 


w 


where  w is  a candidate  value  for 


R - Z 
P 
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From  the  relations  v = 


, w = 


evaluate  /i 


If  (ip  2 0 then  X^ , Y^i^l N)  and  satisfy  the  Kuhn-Tucker 

* * * * 

conditions  and  thus  X.  = X.,Y.  = Y,  , A = A , u = u (i=l,...,N)  and 

ip  i’  ip  i’  p P *P 

problem  p is  solved.  If  < 0 go  to  branch  3. 

» Branch  3.  Solve  for  X^ , Y^  (i“l, . . . ,N) , A , the  (2N+1)  equations 


J b -1  J d -1 

(3a)  L b.h.  . X.J  =^52  d.h.  . Y.J  = r^- 
j-1  J 1 Xp  j-l  J ijP  1 Xp 


(i=l N) 


(3b)  such  that 


N J / b d 

L L h,.  (X.J+Y.J 

IIP  \ i i 


i=l  j=l 


Brief  Discussion  of  Branches 

Branch  1 considers  the  contribution  of  the  total  number  of  canvassers 
Y to  the  primary  accession  requirement  R . If  the  number  of  accessions  Z 
attained  by  Y canvassers  exceeds  the  primary  requirement  R^  then  clearly  the 
total  number  of  canvassers  is  too  large  and  a smaller  number  will  be  calculated 
in  branch  3.  If,  however,  the  accessions  attained  by  canvassers  is  less  than 
the  total  requirement  then  the  associated  values  for  recruiters  will  be  calcu- 
lated in  branch  2 with  the  accession  requirement  now  being  the  value^R^  - Z ^ 
Branch  2 completes  the  candidate  solutions  for  recruiters  X^  and  canvas- 
sers Y^  and  then  tests  the  value  of  the  Lagrange  multiplier  ^ for  non- 
negativity. If  is  negative  the  cost  of  the  Y canvassers  is  too  large 
for  the  accessions  attained  and  branch  3 then  calculates  uhe  appropriate 
new  value  of  total  canvassers.  Notice  that  the  Kuhn-Tucker  conditions 


are  tested  by  branch  2. 
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Branch  3 is  used  if  there  are  no  canvassers  in  the  model  or  if  there 
are  canvassers  and  they  either  cost  too  much  or  are  not  all  required 
because  they  are  able  to  exceed  the  accession  requirements. 


Decision  Branch  Rules  Problem  5 

For  problem  5 the  branch  rules  are  applied  to  each  region  m (m=l,...,5) 
separately.  The  rules  are  the  same.  The  following  modifications  in  notation 
should  be  made.  Primed  relations  denote  changes  to  corresponding  unprimed 
relations  of  the  preceding  section. 

For  (m=l,...,M),  p = 5 

• Branch  1'.  Solve  for  Y.,  i6  I and  v the  equations 
* i m m ^ 


J d -1 

(u> ' & VupV  -v 


m 


such  that 


(lb/  S Y = ip  Y 

i£I  1 m 
m 


If 


16 1 j = l ijP  1 

m 


j 


S £ h, . Y, 

161  j-1  ijP  1 

m 


j 


Z <;  o T_  then  retain  the  value 
m ^m  S _ 

v and  go  to  branch  2 
m 


Z > T.  go  to  branch  3 ' 
iq  mb 


Branch  2'.  Solve  for  X..  i6I  and  W the  equations 
m m ^ 


T b -1 

(2a)  S ^h.  . X. J = W where  W is  a candidate  value  for 

j-l  J iJP  1 " X 


such  that 


(2b) ' £ h X j = jj  T.  - Z 
ijP  i mS  m 

ra 
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r 


From  the  relations  v = ^"^m  , W = ^ 

m m 


evaluate  u 

' fl 


If  ii  s 0 then  X.  , Y.  i£I  and  X , u satisfy  the  Kuhn-Tucker 
'm  1 i m m 


conditions  and  thus  X. 


— * — * — 


X Y.  = Y. , X - X , /x  =M  161 
i*  ip  i ra  m 'm  'm  m 


If  u < 0 go  to  branch  3 ' 
m 


• Branch  3'.  Solve  for  X.,  Y.  i€I  , X , the  equations 
i l m m 


J b .-1 

(3a)'  £ b.h.  . X.3 

j=l  3 *JP  1 


J d.-l 

; S d.h.  . Y.J 

X J-l  J 1JP  1 
P 


such  that 


(3b)  £ 


ijP  \i 


'bj  di 

X.J  + Y3 
. i 1 . 


= <P„Tc 


The  Fibonacci  Search  and  Newton  Iteration 

All  "b"  branch  conditions;  namely  (lb),  2b),  (3b),  (lb)  (2b)  (3b) 

of  the  decision  branch  rules  for  problem  p,  (p=  1 , . . . ,5)  are  solved  indirectly 
by  means  of  finding  the  proper  value  of  the  right  side  values  of  the  "a" 
branch  conditions;  namely  (la),  (2a),  (3a),  (la)',  (2a)  (3a)  The  right 

side  of  the  appropriate  "a"  branch  condition  is  a variable  for  the  Fibonacci 
Search.  Within  any  given  search  the  associated  solutions  in  terms  of  the 
independent  variables  as  functions  of  the  right  side  are  solved  by  a Newton 
method. 
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A detailed  examination  of  the  Fibonacci/Newton  iteration  algorithm 
as  applied  to  problem  p,  (p=l,...,4)  in  branch  -1  follows.  The  iteration 
is  similar  for  all  problems  and  branches  hence  no  further  examination 
of  additional  problems  and  branches  is  needed  to  understand  the  general 
process . 

The  problem  solved  by  Fibonacci  Search  is 


( N 

min  I E Y . - 
l < v < uy  i=i 

subject  to  the  N equations 


f(v) 


J 

£ 

j=l 


d.h. . Y. J 
J IIP  i 


-1 


v (i=l 


. ,N) 


where  L and  U are  the  lower  and  upper  values  of  the  interval 
over  which  the  search  is  made. 


Next,  denote  the  optimizing  value  of  v by  v , the  Fibonacci  Search  iteration 

index  by  k and  the  Newton  iteration  index  by  JL  . The  Fibonacci  search  at 

k k k * k 

iteration  k produces  an  interval  (U  ,L  ) such  that  L < v < U where 
ljk  - L,k  w .618033983^'  (U-L) . This  means,  for  example,  that  if  the 
initial  interval  were  100  units  in  length,  that  after  24  iterations  the 

•k 

length  of  the  interval  containing  v would  be  less  than  .001  units  in 

* 

length.  This  means  v is  known  within  an  accuracy  of  .001  units  after 
24  iterations.  Corresponding  to  each  Fibonacci  iteration  k is  a set  of  N 
complete  Newton  iterations.  The  Newton  iterations  will  be  discussed  after 
the  Fibonacci  process, although  they  actually  occur  within  the  Fibonacci 
iterations . 
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The  Fibonacci  Search  Technique 


L°  - L,  U°  = U,  C2  = ( v/5-1) /2  Cl  = 1 - 


1.  v 


*‘-0,  (</■-!*). 


k+1 

V2  ' C2 


(f  - Lk; 


2.  Compare  f ^v^^jand  using  the  Newton  method  to  evaluate 

the  necessary  function  values  (usually  only  one  of  the  pair) 


3.  If  f 


(vik>  £C 


k+1  _ k k+1  _ k k+1  _ ( v+\  k+l\  k+1  k 

then  L -L,U  “v2,vl  Cl\  ~ U j ' v2  =V1 

Go  to  2. 


4.  If  f v 


'!>  £0lk) 


then  Lk+1  = vxk,  Uk+1  = Uk,  v^1  = v2k,  v2k+1  = C2  (uk+1-  Lk+l 


Go  to  2 . 


I£  4^)-  f(v2k) 

then  if| v^k  - v2k|  is  sufficiently  small 
* / k k \ * 

approximate  v by  (v^  + v2  ) /2  and  f(v  ) a 0 and  the  process  is 

1 k k 1 

finished.  If  | v^  “ v2  I no£-  sufficiently  small  then 

tk  k k k 
L = vx  , U = v2 


Go  to  1, 
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The  Newton  Iterations  V.'ithin  the  Fibonacci  Iteration  k 


Corresponding  to  each  value  vk  (dropping  the  subscript  on  vk)  of 
the  Fibonacci  process  is  a set  of  nonlinear  equations  which  is  solved 
using  a Newton  method.  The  Newton  method  determines  the  value  of  f(v  ) 
in  conjunction  with  the  solution  to  N equations.  The  Newton  method  is 
used  in  step  2 of  the  Fibonacci  Search  to  evaluate  either  f(v  ) or 
f(v£  ).  The  Newton  method  solves  the  following  system  of  N nonlinear 
equations : 


J 

E 

j=l 


d -1 

d h . Y.  2 
J IIP  i 


k 


v 


(i=  1 , . . . ,N) . 


For  each  equation  i the  following  iteration  on  l , the  Newton  iteration 
index,  is  performed. 


a+D  _ Y U) 

i i 


J 


E d.h 
1-1  J 


ijP 


d.-l 

1 


J 

E d (d.-l) 
j = l 2 2 


k 


v 


d.-2 

J 


(i-l,...,N) 


The  Newton  iteration  process  stops  when  I Y.  - Y, 

k 1 i 

sufficiently  small.  The  value  of  f(v  ) at  this  point  is 


a) 


is 


f(vk) 


E Y.  (^  + 1)  - Y | . 
i=l  1 


Control  is  returned  to  step  1 of  the  Fibonacci  Search. 


Figures  3.2,  3.3,  and  3.4  represent  graphically  the  interplay  between 
the  Fibonacci  Search  and  the  Newton  iterations. 
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f axis 


N U+l)_ 
| S Y - Y 

j i 


k+1  k+1 

v,  v. 


iteration  k 
iteration  k + ] 


Fig.  3.2 — Fibonacci  Search 


Figure  3.2  depicts  the  decreasing  interval  length  when  passing  from 

iteration  k to  k+1.  Figure  3.3  corresponds  to  step  3 of  the  Fibonacci 

Search  technique.  It  should  be  noted  that  f (v2k+^=  ^ so  'that 

only  one  value,  namely  need  be  calculated  for  the  iteration  k+1. 

/k+1  k+l^\ 

The  new  interval  (L  , U y is  approximately  .618  times  the  length  of 


the  interval  i L 


iteration. 


The  value  f 


( -r) 


is  obtained  by  using  a Newton 
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3.3  ORAM  OPERATIONAL  DOCUMENTATION 


3.3. 

1 General 

Input  Data 

As  the 

input  routine  for 

the  ORAM  model  is  currently  specified. 

the 

first  11  data  cards  contain  general  information  regarding  processing 

options  desired. 

data  adjustment, 

and  boundary  conditions.  The  informa- 

tion 

required  on 

these  cards  is  detailed  below.  The  user  may  find  it 

helpful  to  refer 

to  the  annotated 

listings  for  the  subroutine  READIN  (App  D). 

Card 

Columns 

Entity 

Remarks 

1-2 

1-130 

ITITLE 

This  title  appears  on  all  reports. 

3 

1-5 

IRUN 

If  IRUN  = 1,  the  program  will  undertake 
a full  optimization.  Otherwise,  only 
an  evaluation  of  the  initial  recruiter/ 
canvasser  placement  will  be  performed. 

4 

1-5 

MAXSTA 

Number  of  DRC's  to  be  considered  (<65). 

6-10 

JTYPES 

Number  of  supply  limited  types  to  be 
considered. 

1 

11-15 

KTERMS 

Number  of  different  terms  of  service 
to  be  considered. 

! 5 

1 

1-5 

ISWTCH 

If  ISWTCH  = 1,  only  recruiters  are  con- 
sidered. If  ISWTCH  = 2,  both  recruiters 
and  canvassers  are  considered. 

| 

6-10 

ISWPCT 

If  ISWPCT  = 1,  a recruiter’s  marginal  pro- 
ductivity with  respect  to  the  supply- 
limited  categories  is  determined  by 
historical  accession  profiles.  Other- 
wise, the  model  will  accept  estimates 
of  these  productivities  by  type.  Data 
indicated  by  ISWPCT  is  input  as  the 
variable  PCTVL. 

6 

1-10 

FUNEMP 

FUNEMP  is  a factor  (near  1,  e.g.,  1.5) 
which  modifies  the  basic  unemployment  rates 
at  which  the  cross-sectional  analysis  was 
made. 
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F.’  CT0R 

FACTOR  adjusts  the  data  for  incompleteness: 

e.g.,  if  data  on  hand  represents  8 months 
of  an  annual  total,  setting  FACTOR  equal 
to  1.5  would  adjust  the  incomplete  data 
to  an  annual  value. 


Card 

Columns 

Entity 

Remarks 

7 

1-10 

SUPLY 

SUPLY  is  the  desired  level  of  supply- 
limited  cohorts. 

11-20 

CANTOT 

CANTOT  is  the  total  canvasser  limit. 

8 

1-10 

11-20 

RECCST 

CANCST 

RECCST,  CANCST  are  the  relative  recruiter 
and  canvasser  unit  costs.  For  best  con- 
vergence properties,  set  to  (e.g.)  1,1 
where  1 may  actually  represent  1 unit 
of  $25,000. 

9 

1-10 

11-20 

FSERV(l) 
FSERV (2) 

FSERV(K)  is  the  fraction  of  total  supply 
with  term  of  service  type  K. 

etc. 

etc. 

10 

1-10 

FCOMBT 

FCOMBT  is  the  fraction  of  total  supply  in 
combat  arms. 

11 

1-10 

11-20 

etc. 

FREGN(l) 

FREGN ( 2 ) 
etc. 

FREGN(K)  is  the  fraction  of  total  supply 
in  region  K. 

3.3.2 

DRC  Specific  Input 

Data 

As  the 

program  is 

currently  structured,  DRC  specific  data  is  input 

beginning  with 

card  number 

12.  The  types  of  data  currently  accepted  by  the 

program  are  listed  in  Table  3.1. 

Table  3.1 


VARIABLES  ACCEPTED  FOR  ANALYSIS  BY  ORAM 


Identifier 

QME 

WAGE 

UNEMP 

ATT 

ADV 

RECR 

CANV 

PETVL 


Definitions 

Qualified  military  populations  by  cohort  type 

Ratio  of  military  to  civilian  range  by  cohort  type 

Unemployment  by  cohort  type 

Attitude  measure  by  cohort  type 

Advertising  measure  (e.g.,  expenditures)  by 
cohort  type 

Recruiters  on  station  by  DRC 
Canvassers  on  station  by  DRC 

If  ISWPCT  = 1,  PCTVL  is  the  historical  accession 
file  described.  If  marginal  productivities  for 
recruiters  by  cohort  type  are  to  be  estimated  in 
some  other  way,  the  estimates  are  input  as  PCTVL 
with  ISWPCT  f 1. 


3-24 

LJ  J 
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Currently  ORAM  assigns  the  same  value  of  QME,  WAGE,  UNEMP,  ATT, 
ADV,  RECR,  and  CANV  to  each  cohort  type.  This  value  is  the  first  and 
only  value  input  for  each  DRC.  The  program  is  easily  modified  to  allow 
for  input  differentiated  by  cohort  type  (see  annotated  program  listing 
for  READIN  lines  100-175).  The  general  input  structure  for  the  variables 
QME,  WAGE,  UNEMP,  ATT,  ADV,  RECR  and  CANV  is  as  follows,  where  as  above 


MAXSTA  is 

the  number 

of  DRC's  to 

be  considered. 

Card 

Columns 

Entity 

Remarks 

1 

1-5 

IDENTIFIER 

Left-justified  selection  from  the 
following:  QME,  WAGE,  UNEMP,  ATT, 
ADV,  RECR,  CANV 

2 thru 

1-5 

DRC  identi- 
fier 

MAXSTA+1 

6-7 

DRC  index 

DRCs  are  assigned  consecutive  indices 
beginning  with  1. 

8-9 

Filler 

10 

Region  num- 
ber 

Currently  either  1,3, 4, 5 or  6 

11-18 

VAL(l) 

Value  of  the  variable  for  cohort  type  1 

19-26 

VAL (2) 

Value  of  the  variable  for  cohort  type  2 

etc . 

etc. 

etc . 

MAXSTA+2 

1-3 

’END’ 

Signifies  end  of  data  for  the  variable 

under  consideration 

As  noted  above,  at  present  VAL(l)  is  assigned  to  all  cohort  types 
and  is  accordingly  the  only  value  input.  A data  card  must  appear  for  each 
of  the  MAXSTA  DRCs  under  consideration.  The  above  format  is  repeated  for 
each  variable  to  be  input. 

The  marginal  productivity  for  a recruiter  by  cohort  type  is  cur- 
rently estimated  from  the  historical  accession  file  described  in  App  C. 
These  data  are  input  under  PCTVL  as  follows. 
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Card 


Columns 


Values 


First  1-5  1 PCTVL ' 

Following  Historical  accession  data 

file 

Last  1-3  'END' 

(In  the  event  that  a direct  input  of  marginal  recruiter  productivities 
by  cohort  type  is  required,  some  reprogramming  of  READIN  may  be  desirable.) 

3.3.3  Inputting  Model  Parameters 

For  convenience,  the  following  definitions  will  apply  in  this 

section: 

X^:  The  K independent  variables  will  be  denoted  X^  for 

1 = 1,2,...,K.  E.g.,  in  the  sample  problem  X^  corresponds  to 

QMA,  X^  to  wages,  etc. 

E(I,J):  The  elasticity  of  variable  X^  for  type  J is  denoted 

E(I,J).  In  the  sample  problem  E(1,J)  = E(2,J)  = 0 for  all  J. 

CON(J):  The  value  of  the  scaling  constant  for  the  model 
representing  cohort  type  J is  denoted  CON(J). 

Y(J):  The  number  of  accessions  arising  from  cohort  type  J 

is  denoted  Y(J). 

In  the  sample  problem  and  for  cohort  type  J,  ORAM  considers  the 
following  model  (described  in  standard  FORTRAN  notation): 

Y(J)  = EXP**CON(J) *X(1) **E (1, J) *X(2) **E (2 , 1) * 
X(3)**EU,J)*X(4)**E(4,J)*X(5)**E(5,J) 

*X(6)**E(6, J) 

where 


X(1) 

is 

QMA, 

X(2) 

is 

the  ratio  of  military  to  civilian  wages, 

X(3) 

is 

the  general  unemployment  rate. 

X(4) 

is 

an  advertising  measure. 

X(5) 

is 

an  attitude  measure 
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X(6)  is  recruiters  on  station, 

X(7)  is  canvassers  on  station, 

CON(J)  and  E(I,J)  are  as  defined  above. 

Values  for  E(I,J)  are  input  as  follows: 


Card 

Columns 

Entity 

Value 

Remarks 

1 

1-5 

Header 

'ELAST' 

2 

1-5 

Identifier 

Optional 

6-7 

Variable 

I 

number 

8-11 

Filler 

Spaces 

12-19 

VAL(l) 

E(I,1) 

20-27 

VAL (2) 

E(I, 2) 

etc. 

etc. 

etc. 

3 thru 
8 

Elasticities 

for  the  remaining  variables  as 

per  format 

9 

1-3 

Trailer 

’END' 

In  its  present  form  READIN 

expects  to  be  provided  with 

sets  of 

elasticities 

(corresponding  to  X ( 1 ) - X(7)). 

Values  for  CON(J)  are  input 

as  follows: 

Card 

Columns 

Entity 

Value 

Remarks 

1 

1-5 

Header 

'CONST' 

2 

1-5 

Identifier 

Optional 

6-1 1 

Filler 

Space 

12-19 

VAL(l) 

C0N(1) 

20-27 

VAL(2) 

CON (2) 

etc. 

etc. 

etc. 

3 

1-3 

Trailer 

'END' 

In  the  sample  problem,  constants  are  provided  for  each  of  the 
eight  cohort  types  considered.  For  every  type,  the  constant  is  = .601068; 
but  as  noted  above,  the  constants  need  not  be  identical. 
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1.4  The  Nonlinear  Programming  Model — Optimal  Recruiter  Allocation 
Model  (ORAM) 

Given  all  the  established  input  conditions,  the  ORAM  proceeds  to 
solve  for  the  optimal  allocation  of  recruiter/canvassers  and  objective 
allocations  to  each  DRC  for  a set  of  primary  accession  requirements.  The 
ORAM,  in  effect,  solves  five  complete  optimization  problems  per  computer 
run.  Specifically,  given  the  input  conditions  of  the  model,  what  is  the 
most  cost-effective  allocation  of  recruiters/canvassers  if: 

1.  A total  national  goal  of  X number  of  accessions  is  required? 

2.  A total  national  goal  of  X number  of  accessions  with  3-year 
terms  of  service  is  required? 

3.  A total  national  goal  of  X number  of  accession^  with  4-year 
and  greater  term  of  service  is  required? 

4.  A total  national  goal  of  X number  of  accessions  with  combat 
arms  specialty  is  required? 

5.  A regional  goal  of  0 percait  of  the  X number  of  total  accessions 

m — 

is  required  in  region  m(m=l , . . . ,m) ? 

3.3.5  ORAM  Model  Outputs 

The  optimal  allocation  solutions  to  the  ORAM  are  used  as  a basis 
for  a series  of  output  report  evaluations.  Each  evaluation  is  fairly 
detailed  and  provides  the  following  information  by  DRC: 

• DRC  Identification 

• DRC  Number 

• Number  of  Recruiters  on  Station 

i«  Number  of  Canvassers  (where  applicable) 

• Number  of  Accessions  (total) 

■ • Number  of  Accessions  by  Cohort  Group 

In  addition,  subtotals  by  region,  totals  nationally,  percent  regional  and 
national  summary  tables  are  printed. 
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All  of  the  above  information  is  evaluated  for  each  of  five  optimal 

solutions  and  for  each  of  five  accession  content  situations  within  each 

optimal  solution.  Figure  3-5 , Output  Report  Guide,  depicts  the  range  of 

evaluations  which  are  made.  Each  cell  represents  a particular  set  of 

conditions  for  which  a complete  evaluation  is  made.  Figure  3-6  presents 

one  of  the  25  output  reports  printed  after  the  five  ORAM  optimizations 

are  completed.  The  sample  report  depicts  the  combat  arms  content  for 

CAT123AHSG  black  and  white  when  the  overall  strategy  is  to  meet  a total 

* 

national  requirement  (p-1)  of  54568  CAT123AHSG  accessions.  Using  the 
notation  of  Fig.  3-5,  the  sample  report  is  R^.  The  values  of  the  para- 
meters determined  from  the  cross  section  analysis  are  as  follows: 


recruiters 

°61  = 

“62  = 

.835798 

unemployment 

a31  = 

a32  = 

.41938 

atti tudes 

a41  " 

a42  = 

.724 

constant 

n 

fH 

> 

V2  " 

1.5  exp  (-.908979)  = .604403 

These  values  of  the  parameters  yield  the  following  accession  equation 
for  DRC  i: 

For  black  content: 

CAT123AHSGBi  = (%  black^**  (.604403)  (Unemployment^ ' 419138 
. (Attitudes^)  . (Recruiters^) 

For  white  content: 

CAT123AHSGW±  = (%  white^**  (.604403)  (Unemployment ) * 419138 

^ .. 72482  , . 835798 

. (Attitudes^  . (Recruiters^ 

For  any  particular  DRC  the  appropriate  values  of  the  environmental 
factors  determine  the  coefficients  of  the  recruiter  accession  equation. 

For  example,  the  enlistment  supply  equation  for  DRC  1A,  Albany,  in  region 
one  is  evaluated  as  follows  for  CAT123AHSG  black  and  white. 


* 

This  is  based  upon  CAT  I-II1A,  HSG  being  the  supply  limited  group. 
If  the  mix  of  the  supply-limited  group  changes,  ORAM’s  coverage  must 
also  be  changed  to  encompass  all  of  the  supply-limited  categories. 

Historical  percentage  breakout  by  DRC. 
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v Evaluation  of 
\ Accession 
\ Content (q) 

\ ^ q=  1 

q=2 

q=  3 

q=4 

q = 5 

Optimal  \ 
Solution  (p\ 

National 

3-Yr  Term 

4-Yr+  Term 

Combat  Arms 

Regional 

p=l 

National 

National 

Optimal 

1 

j 

p=2 

3-Yr  Term 

* 

R 

pq 

p=3 

4-Yr+  Term 

1 1 

p=4 

Combat  Arms 

Combat  Arms 
Optimal 

p=5 

Regional  i 

‘ 

Regional 

Optimal 

R is  the  evaluation  report  of  accession  content  q based  on  the 
optimal  allocation  of  recruiters/canvassers  to  primary  requirement  p. 

Example : is  the  evaluation  report  of  the  number  123AHSG  acces- 

sions enlisting  for  four  or  more  years  if  the  recruiters/canvassers  are 
allocated  to  achieve  a specified  three-year  term  of  service  objective  at 
minimum  cost.  (See  Fig.  3-6  for  an  example  of  output  report  R^.) 

Fig.  3-5 — Output  Report  Guide 
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STRAT'OY  - “IN  COST  SUOJEST  TO  NATIONAL  TOTAL  ACCESSION  RrOLIRFMENT  54  56  A. 

SUPPLY  LINITEO  TYPES  1,2-  CAT12TAMSG  16U0J  TOTAL  RFO.  WITH  3T.A  PERCENT  «UPPLY  RECR.  ELASTICITY  « .815  . UNENPL 


Fig.  3-6  — Sample  Output  Report 
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Fig.  3-6  (continued) 


DRC  1A,  Albany,  N.Y.,  Region  one 


Unemployment  rate  8.3% 

Attitude  factor  40.0% 

Portion  black  2% 

Portion  white  98.0% 

Then  for  Albany: 


CAT123AHSGB  = (.02)  (.604403)  (8. 3) *419138. (40) * 72482. (Recruiters  ) -835798 
= .42  (Recruiters)'835798 

CAT123AHSGW  = (. 98) (. 604403)  (8 . 3) ' 419138 . (40) ' 72482 . (Recruiters  )' 835798  i 

= 20.84  (Recruiters)'835798 

The  cross-section  analysis  average  unemployment  rate  was  7.1%  for  FY75. 

If  unemployment  were  increased  to  7.7%  or  1.0845  times  the  average,  then  the 

coefficients  are  multiplied  by  (1.085) '4*9*38^  = 1.0345.  Thus  the  equations 

used  for  DRC  1A  in  region  one  for  an  average  unemployment  of  7.7%  are: 

835798 

CAT123AHSGB1  = .44  (Recruiters^  ' 

835798 

CAT123AHSGW1  = 21.56  (Recruiters^ 

Given  46  recruiters  the  expected  CAT123AHSGB  content  at  DRC  1A  is  about  10 
and  for  the  equivalent  white  content  529. 

These  enlistment  supply  equations  for  all  other  DRCs  are  evaluated  in 
a similar  manner.  Specific  content  such  as  CAT123AHSG  combat  arms,  for 
example,  is  evaluated  using  historical  proportions  to  modify  the  basic 
enlistment  supply  equations. 

A brief  discussion  of  Fig.  3-6 — Sample  Output  Report  R follows: 

The  overall  optimization  strategy  (National)  is  indicated  in  the  head- 
ing. The  minimum  cost  in  this  case  is  the  number  of  recruiters  since  unit 
cost  is  directly  related  to  recruiters.  The  national  total  accession  require- 
ment for  CAT123AHSG  was  54,568  which  constituted  almost  34%  of  a total  require- 
ment of  161,000  accessions.  The  recruiter  elasticity  and  average  unemployment 
rate  are  announced  at  .835  and  7.7%,  respectively. 

A detailed  breakout  of  the  CAT123AHSG  combat  arms  accessions  is  given 
DRC  by  DRC  with  subtotals  by  region  and  grand  total  nationally.  Type  one 
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1 


is  CAT123AHSG  black  and  type  two  is  CAT123AHSG  white.  Notice  that  a total 
recruiter  force  of  A, 423  is  required  to  meet  the  national  goal  of  54568 
total  CAT123AHSG  accessions. 


r 


A percent  summary  table  completes  the  report.  This  section  of  the 
output  report  describes  the  complete  Regional  distribution  of  the  optimal 
national  recruiter  allocation  and  associated  breakout  by  type.  It  can 
be  noted  that  if  the  objective  was  to  find  the  optimal  recruiter  force 
allocation  to  meet  a specified  CAT123AHSG  combat  arms  requirement  (Report 
R^)  that  fewer  recruiters  would  be  required  to  attain  the  13132  combat 
arms  accessions.  In  fact,  a requirement  of  12400  CAT123AHSG  combat  arms 
accessions  can  be  met  with  a recruiter  force  of  3421.  However,  this 
recruiter  allocation  only  yields  a national  total  of  39393  total  CAT123AHSG 
accessions. 

There  are  several  different  uses  of  the  ORAM — some  of  which  are  men- 
tioned briefly  here.  Perhaps  foremost  on  the  list  of  uses  of  the  ORAM  is 
assessing  the  effects  that  a range  of  requirements  for  specified  accession 
groups  have  on  determining  the  recruiting  force  size  and  mix.  Another 
important  use  is  concerned  with  the  effects  that  changing  environmental 
factors  have  on  the  recruiting  force  and  the  accession  groups.  Changes 
in  unemployment,  for  example,  modify  the  productivity  of  the  recruiter.  The 
output  reports  discussed  above  assumed  an  average  unemployment  rate  of  7.7%. 
If  the  unemployment  rate  were  higher  a slightly  smaller  recruiter  force  would 
have  been  required.  Finally,  mention  of  the  potential  use  of  the  ORAM  for 
budgeting  uses  should  not  be  overlooked.  Departures  from  or  modifications 

budgets  for  the  recruiter /canvasser  programs  can  be  investigated 
ORAM. 


of  basepoint 
by  using  the 


I 

4 THE  MERGE/BMD  TIME  SERIES  REGRESSION  MODEL 

4 . 1 GENERAL 

The  MERGE/BMD  time  series  regression  model  has  evolved  to  its  pre- 
sent form  through  a series  of  studies  concerned  with  the  accession  programs 
of  the  military  services,.  It  is  structured  around  an  extensively  modified 
version  of  the  UCLA  Biomed  stepwise  regression  program.  These  modifications 
have  been  concerned  with  the  application  of  regression  analysis  to  time 
series  data  and,  in  particular,  to  military  accession  data.  The  model  dis- 
cussed in  this  chapter  is  presently  operational  on  the  USAREC  UNIVAC  1108 
computer  system.  The  operational  aspects  of  the  documentation  pertain  spec- 
ifically to  that  system. 

4.2  METHODOLOGICAL  CONSIDERATION  FOR  MODEL  APPLICATION  TO  U.S.  ARMY 

ACCESSION  ANALYSIS 

The  regression  model  that  has  been  developed  for  accession  analysis 
is  a nonlinear  regression  model.  The  analysis  assumes  the  dependent  variable 
time  series  is  initially  in  unadjusted  or  "seasonalized"  form  and  the  inde- 
pendent variables  do  not  reflect  the  seasonal  variation  apparent  in  the 
observed  values  of  the  dependent  series.  Yearly  data  which  are  subject  to 
periodic  fluctuations  are  of  this  type.  The  model  is  solved  yielding  the 
adjustment  factors  (seasonal  factors)  simultaneously  with  the  parameters  of 
the  ordinary  linear  model.  Two  models  result  from  the  above  analysis:  the 

first  looks  upon  the  solution  as  a deseasonalization  of  the  dependent  variable 
and  yields  a least  squares  fit  to  the  deseasonalized  variable  (this  is  the 
model  solved  by  the  program);  the  second  considers  the  adjustment  factors  as 
weights  for  the  linear  portion  of  the  model  and  is  used  for  prediction 
purposes. 

4.3  THE  MODEL 

(Dl)  A dependent  variable  Y (a  (n  x 1)  vector)  which  is  typically 
a time  series 

(D2)  An  independent  variable  data  set  X (a  (n  x p)  matrix) 

(Al)  The  number  of  observations  of  the  (seasonalized)  dependent 
variable  Y is  considerably  longer  than  one  cycle  of  the  seasonal  period 
(say  at  least  two  cycles) 

Dixon,  W.  J. , ed. , "BMD  Biomedical  Computer  Programs,"  Health  Sciences 
Computing  Facilities,  University  of  California,  Los  Angeles,  1964. 
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(A2)  The  data  set  X does  not  reflect  the  seasonal  variation  con- 
tained in  the  observed  values  of  Y 

Two  models,  Ml  and  M2,  with  separate  objective  functions  can  be 
developed  from  the  above  assumptions.  They  are: 


(Ml)  minimise  f (S'1,  p)  = jlS-1^  - X3II  2 

where  S “is  a (n  x n) diagonal  deseasonalization  matrix 
3 is  a (p  x 1)  parameter  vector 


(M2)  minimize  f (S,  =)  = j|Y  - SXp||  d 

S,  P 

Ordinary  linear  regression  analysis  assumes  S is  the  identity  matrix 


and  one  solves  for  the  parameter  vector  3.  This  is  not  the  case  here 
however.  Model  Mi  assumes  that  Y is  deseasonalized  by  3 “ and  yields  a 
least  squares  fit  of  Xg  to  the  deseasonalized  vector  S “Y.  Model  M2  con- 
siders S to  he  an  adjustment  matrix  which  in  turn  modifies  the  linear 
portion  X;  to  form  SX3  as  a least  squares  approximation  to  Y.  Both 
models  are  nonlinear  since  and  are  nonlinear  functions  over  the 
full  (S,g)  space. 


For  simplicit;.  , the  diagonal  elements  (seasonal  factors) 

S..  of  S reolace  an  earlier  notation  and  are  defined  by 
11 


sii"n  exp  {Vij} 

j-i 

where  5 . = ( 1 t , and 

lj  ( O irj 


is  the  seasonal  coefficient  for  month  i 


INPUTS 


Dependent  tine 
series  with 
seasonal  variation 


Independent  time 
(data)  series  - not 
necessarily  ccntainin 
seasonal  variations 


Stepwise  regression 
on  deseascnalized 
series  and  independent 
variables  - Model  Ml 


Further  deseascnali- 
zation  of  dependent 
time  series 


Convergence 

criterion 

satisfied 


OUTPUTS 


1.  Deseascnali cation  factors 

2.  Parameter  values  for  these 

variables  in  the  model 
3-  Residual  errors,  etc. 

4.  Predictions  using  Model  M2 


Fig.  4-1  — The  Deseasonalized  Step  Regression  Model 

Basic  Cycle 
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4.4  THE  SOLUTION  PROCEDURE 


A sequence  of  deseasor.alizafcion , and  stepwise  regressions  is  used 

to  solve  model  Ml  . The  solution  procedure  is  iterative.  Iteration  i 

, i.-l  i 

generates  the  deseasonalization  matrix  (S  ) and  parameter  vector  g 

1 . . 

which  alternately  satisfy  the  stationary  conditions  G“(f. ) = 0 and 

P 1- 

G~  i ( f . ) = 0 [G  (f)  denoting  the  gradient  vector  of  f in  the  w subsrace]  . 

S~J-  l . w 

The  full  space  joint  condition  GX_i  ( f . ) = 0 does  not  hold  during  the 

S31. 

course  of  the  iterations  since  in  general  if  G1_t( f , ) = 0 then  G1(f1)  * 0 

i * /i\-1  W ir“I 

and  conversely.  However  as  g approaches  p and  (S  ) approaches  5 , 

*•- 1 * 

G 1 (f. ) approaches  G , (f.,)  which  vanishes  at  (S  , g ) . A complete 

S"  g S*'  0 - 

stepwise  regression  analysis  is  performed  on  the  deseasonalized  variable 
/ U'1 

(S  ) Y during  each  iteration  i.  Independent  variables  which  are  ir.  the 
solution  at  iteration  i are  not  necessarily  in  the  solution  at  iteration 


i + 1. 


•The  iteration  eouations  for  model  Ml  are: 


(11)  p1  = (X'X)“  X'(Si_1)’_Y 

(12)  InQd  r Diag{[(5i)  Y - XpTjY  } 1 = 0 (matrix) 


where  S = I 


—f 

1 is  a column  vector  of  ones 


I , is  a n matrix  of  rank  r containing  [n/r]  full  id; 
mod.  r s'-/-! 

matrices  of  rank  r and  a residual  identity  submatrix  of 
n - [n/r]r 

( r is  the  period  of  cr.e  cycle. 

The  constraint  (C)  below  is  arbitrarily  imposed  on  the  deseas; 
elements  to  allow  uniformity  of  solutions  for  comparitive  pur: 


ilizaoion 


k+r 

(C)  Z S = r 
j=k 

Constraint  (C)  means  that  the  sum  of  the  diagonal  elements  of  S over  any 

cycle  is  normalized  so  that  the  parameter  vector  3 vill  absorb  the  major 

part  of  the  solution,  and  the  average  value  of  3 ever  one  cycle  is  1/r. 

J J 

The  iteration  process  is  terminated  when  condition  (T)  below  holds. 

(T)  | ^ - IQ|<_Diag  ^element  by  element 

where  £ is  prescribed  (currently  .001). 

Table  4.1  contains  the  gradient  vectors  and  Hessian  matrices  for 
and  fg  of  models  Ml  and  M2  respectively,  f has  a positive  semi- 
definite  Session  for  non-negative  data  and  thus  model  Ml  can  be  solved 
as  a convex  problem.  This  is  not  the  case  with  model  M2. 
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4.5  REGRESSION  PROBLEMS  WITH  FIXED  ELASTICITIES 

If  each  variable  appearing  in  the  regression  model  is  normalized  by 
dividing  each  variable  (including  the  dependent  variable)  by  its  mean 
value,  then  the  following  corresponding  coefficients  e^  are  called  the 
elasticities  for  the  kth  variable.  The  deseasonalized  regression  model 
becomes 

S-IY  - eQ  + ex  _X^  + . . . + el  XP 
S~^Y  X1  ~X? 


where  the  superscript  K on  X indexes  the  vector  X with  n components 
(observations)  and  S ^ is  the  nxn  diagonal  deseasonalizing  matrix. 

Then,  if  e 4 o 


where  E is  the  expected  value  operator,  thus 


The  relation  between  e^  and  the  ordinary  coefficients  *-s  as  follows: 

For  S-1Y  - b + £ & ^ 

° k-1  * 


if  3^o  then 
^ o 


(k  =*  o...  .,p) 


3 S-LY  e,  3.  ^ 

and  a k or  =*  : k 

x*  s-1y 

Frequently  it  is  desired  to  fix  one  or  more  of  the  elasticities  of  a 

set  of  variables  and  formulate  a regression  problem  on  the  remaining  set 

of  variables.  Such  problems  are  considered  in  the  following  manner. 

If  we  assume  q elasticities  e^  (k=l,...q),q  < p, are  prescribed,  then 


'h  - j3q  + (”ei  S-1Y  X1  + ... 


+ e S 

q 


i. xq) + ^q+ixq+1  +-  •• + 


is  the  deseasonalized  model  with  q fixed  elasticities  and  p-q  independent 
variables . 

The  coefficients  and  deseasonalized  matrix  S 1 are  found 

by  solving  the  equation 

min  F = L (Z.  - Y.)2  /3q+1>...  S_1 


where  Z.  = S . ^ Y.  -S  ^Y 
i ii 


*A 


X' 


e X? 
i o . q i 


* -<»„  +Vix,+1  + •••  *e/- 


4.6  MERGE /BMD  TIME  SERIES  REGRESSION  MODEL  OPERATIONAL  DOCUMENTATION 

General  system  flow  diagrams  for  PROGRAM  MERGE  (Figure  4„2)  and 
PROGRAM  BMD  (Figure  4.3)  are  presented  here  at  a level  sufficient  to  describe 
the  general  input-output  flow  of  information.  Tabular  descriptions  of  the 
I/O  units  used  by  PROGRAM  MERGE  (Table  4.2)  and  PROGRAM  BMD  (Table  4.3)  are 
also  included. 
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Time  series 
Daca  bank 


See  ,mit  descriptions  for 
details 


Fig.  4.2—  PROGRAM  MERGE  General  Flow  Diagram 
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4.7  BMD  CONTROL  LANGUAGE 

A detailed  description  of  the  BMD  CONTROL  LANGUAGE  and  its  use  is 
contained  in  this  section.  BMD  CONTROL  LANGUAGE  WORDS  are  listed  in 
Table  4.4. 


Table  4.4 

BMD  CONTROL  LANGUAGE  IN  ORDER  OF  OCCURRENCE 
IN  A PROBLEM  SPECIFICATION 


Card 

Type 

Remarks 

1 

PROBLM 

See  Table  4.5 

2 

LABELS 

3 

format 

(10X, 10F12.4/ (10X, 10F12.4) ) 

4 

SUBPRO 

See  Table  4.6 

5 

CONDEL 

See  Table  4. 7 

6 

blank  card 

Mandatory 

7 

IDXPLT 

Optional 

8 

FINISH 

Cards  4-7  are  repeated  for  each 

additional  subproblem. 

Specifying 

the  BMD  Control  Cards 

BMD  control  cards  will  be  entered  in  the  following  sequence,  either 
initially  or  following  the  control  card  RESTART  (key  word).  That  is,  these 
cards  are  read  once  initially  and  used  for  every  regression  problem  until 
a RESTART  card  is  read  (p.  4-31).  The  control  cards  are  unchanged  from 
their  BMD  specifications,  with  the  exception  of  column  70,  which  still 
applies;  however,  they  are  revised  automatically  for  every  run  by  MERGE 
to  incorporate  the  data  counts  and  label  counts.  The  actual  cards  generated 
are  listed  in  MERGE  output. 

PROBLM  Card 

This  should  be  completed  according  to  Table  4.5.  The  variable  count, 
observation  count,  orediction  count  and  labels  count  are  rewritten  to  re- 
flect the  contents  of  the  data  files  and  the  variables  used  in  this  run. 

The  prediction  count  reflects  the  remainder  of  the  54  points  in  each  time 
series,  which  occur  after  the  last  actual  observation.  All  data  is  stretched 
using  sample  and  hold  in  (p.  4-21),  and  residual  calculation  is  used  to 
predict  the  dependent  variable  for  the  values  of  the  independent  variable. 
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Table  4.5 


1 


SPECIFICATIONS  FOR  PROBLM  CARD 


Column 

Entity 

Remarks 

1-6 

'PROBLM' 

Mandatory 

7 

Report  Title 

Set  equal  to  1 

8 

Filler 

Blank 

9 

Print  control 

3 suppresses  intermediate  output; 
2 or  less  prints  details  of  each 
step  in  the  regression 

10-15 

Problem  name 

6 alphanumeric  characters 

17-20 

Sample  size 

Leave  blank.  Will  be  computed  by  MERGE 

24-25 

p,  number  of  original  vari- 
ables 

Leave  blank.  Will  be  computed  by  MERGE 

29-30 

m,  number  of  TRNGEN  cards 

0 < m < 99 

34-35 

q,  number  of  generated 
variables 

- 9 < q < 78 

39-40 

Unit  number  if  data  is 
on  an  input  medium  other 
than  cards 

blank  - if  card  input;  13  - otherwise 

44-45 

s,  number  of  subproblem 
cards 

1 < s < 99 

48-49 

Number  of  variables 

Leave  blank.  Will  be  computed  by 

labeled  on  LABELS  card 

MERGE 

51-53 

Output  switcn 

'YES' , if  means  and  standard  devia- 
tions are  to  be  printed;  blank,  otherwise. 

55-57 

Output  switch 

'YES' , if  covariance  matrix  is  to  be 
printed;  blank,  otherwise 

59-61 

Output  switch 

'YES' , if  correlation  matrix  is  to  be 
printed;  blank,  otherwise. 

63-65 

Output  switch 

'YES'  , if  zero  regression  intercept 
is  desired;  blank,  otherwise 

68-69 

Input  switch 

'NO',  if  tape  specified  in  columns  39-40 
is  not  to  be  rewound  before  this  problem; 
blank,  if  columns  39-40  are  blank  or  Y 
tape  rewind  is  desired. 

70 

Operational  mode 

See  Section  4.9 

■ 
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Table  4.6 


SPECIFICATIONS  FOR  SUBPRO  CARD 

Column  Entity  Remarks 

i-6  Identifier  'SUBPRO'  (mandatory) 

9-10  Dependent  variable  This  is  the  variable  number  assigned 

ID  in  MERGE  processing.  Ordinarily  the 

dependent  variable  should  be  the  first 
variable  in  the  variable  list  and  accord- 
ingly '1'  should  appear  in  column  10. 

13-15  Maximum  number  of  If  this  is  left  blank,  the  default  is  one 

steps  allowed  in  step  per  variable.  Both  original  vari- 

stepwise  regression  ables  and  variables  generated  by  lags 

or  'CREATE'  are  considered  in  computing 
the  default  value. 

2Q-25  F-level  for  inclusion  This  will  be  .01  if  left  blank. 

30-35  F-level  for  deletion  This  will  be  .005  if  left  blank. 

40-45  Tolerance  level  for  This  will  be  .001  if  left  blank, 

comparisons 

49-50  Number  of  variables  This  must  be  between  0 and  30. 

on  the  index-plot 
card 

53-55  CONDEL  switch  'YES'  if  CONDEL  cards  are  included 

58-60  Residuals  switch  'YES'  if  a list  of  residuals  is  to 

be  printed 

63-65  Summary  switch  'YES'  is  a summary  table  is  to  be 

printed 
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Table  4.7 


SPECIFICATIONS  FOR  CONDEL  CARD 


! 


Column 

Entity 

Remarks 

1-6 

Identifier 

'CONDEL'  (mandatory) 

7 

■k 

Control  value  for 

1st  variable 

8 

k 

Control  value  for 

• • • 

2nd  variable 

72 

* 

Control  value  for 

66th  variable 



. 

Control  values  are 

1 - Delete  variable  from  consideration 

2 - Free  variable 

3 - Low-level  forced  variable 

4 - Higher-level  forced  variable 


9 - Highest-level  forced  variable 

Values  3-9  allow  prioritizing  of  the  entry  of  variables  into 
the  regression 
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LABELS  Card 


A blank  LABELS  card  must  be  included  to  signal  where  the  labels  are 
required  in  the  deck.  The  appropriate  labels  are  generated  from  the  vari- 
able specified  by  using  the  first  6 characters  of  each  name. 

Format  Card 

The  format  card  is  fixed  by  MERGE  and  should  be: 

(10X,  10F12.4/(10X,  10F12.4) ) 
in  columns  1 - 27 . 

SUBPRO  Card 

The  SUBPRO  card  snould  be  completed  according  to  Table  4.6.  It  is 
convenient  to  make  the  dependent  variable  1,  but  this  is  not  a restriction. 

CONDEL  Card 

The  CONDEL  card  follows  Table  4.7  specifications.  The  card  also 
signals  the  creation  by  MERGE  of  2 additional  cards  defined  below,  both 
of  which  are  new. 

LIMITV  Card 

The  LIMITV  card  notes  the  variable  count  limits,  if  any. 

SIGNS  Card  (Coefficient  Sign  Restriction) 

The  SIGNS  card  records  the  sign  controls  for  the  variables,  if  any. 

TITLE  Card 

In  addition  there  is  now  a title  card,  which  may  be  blank,  and  permits 
step  print  suppression  and  provides  a run  title  for  all  runs.  If  the  card 
is  blank  the  independent  variable  title  is  used. 


Cols 

Entity 

Remarks 

1-2 

Print  control 

Specifies  step  to  store  printing 
(only  valid  if  print  not  suppressed 
on  problem  card,  Table  4.5) 

3-74 

Title 

Leave  blank  unless  desired  title 
different  than  dependent  variable 

Example:  (03  C123  - combat  arms  categories,  1,  2 and  3) 
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This  would  have  the  effect  of  starting  the  regression  printout  at  Step  3 
and  provides  a title  printed  bold-face  across  the  leading  page. 

IDXPLT  Card  (optional) 

♦ The  IDXPLT  card  follows  Table  4.8  specifications.  This  card  is 

required  if  columns  49-50  of  the  subproblem  card  are  greater  than  0. 

FINISH  Card 

The  file  is  terminated  by  a FINISH  card. 


Table  4.8 

SPECIFICATIONS  FOR  IDXPLT  CARD 


Column 

Entity 

Remarks 

1-6 

Identifier 

'IDXPLT'  (mandatory) 

7-8 

Number  of  1st  variable  to 
be  plotted 

9-10 

Number  of  2nd  variable  to 
be  plotted 

65-66 

• • • 

Number  of  30th  variable 
* 

to  be  plotted 

No  more  than  30  variables  may  be  plotted  per  subproblem. 
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4.8  MERGE  CONTROL  LANGUAGE 

The  MERGE  Control  Language  is  listed  in  Table  4.9.  A detailed  dis- 
cussion of  the  language  and  its  use  is  presented  in  this  section. 

Defining  a Variable  Time  Series  for  MERGE 

MERGE  recognized  any  appropriate  data  defined  by  the  rules  of  this 
section  as  a legitimate  time  series. 

Input  Data  for  the  Variable  Time  Series 

This  is  done  by  card  (or  card  image)  input  and  is  read  in  on  input 
Unit  5.  The  first  card  defines  the  first  data  occurring  in  any  time  series, 
and  a list/no-list  option. 

Time  Series  Initialization  Card 


Column 

Entity 

Remarks 

1-20 

Filler 

May  be  used  as  data  stream  identifier 
field 

21-25 

Year 

Starting  year  of  all  time  series 
(right  justified) 

26-30 

Month 

Starting  month  of  all  time  series 
(right  justified) 

31-34 

Filler 

36-40 

List  option 

Right  justified  values: 

LIST  - causes  listing  of  all 
time  series 

NOLIST  - suppress  listing 

41-80 

Filler 

f 

{ 


•- 
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Table  4.9 


MERGE  CONTROL  LANGUAGE 
(order  of  precedence) 


Time  Series  Initialization  Card 
Title  card  for  time  series  input 
Time  series  data  input 
END 

LAST  FILE 

(BMD  Control  Cards) 

BASE 

CREATE 

Variable  creation  cards 
TREND  (optional) 

END 

Dependent  variable  name  (Table  4.10)  1 
Independent  variables  (Table  4. 10)  ^ 
DUMMY  (Additive  optional  seasonals) 

SET 

Independent  variables  to  be  included 
in  subsequent  runs 

END 

NET;  DEPVAR  (optional) 

New  independent  variable  (Table  4. 10) 
END 

RESTART  (optional) 

New  variable  set  (dependent  and 

independent  variable)  (Table  4.10) 

END 

RESET  BASE  (optional) 

New  BASE  card 

END 

STOP 


Repeated  for  every  time 
series  (max.  99) 


Repeated  for  additional 
sets  of  variables 


4-20 


Each  time  series  consists  of  at  least  three  cards: 

(a)  The  first  card  is  a title  card  used  to  title  the  printout  if 
a listing  is  desired. 

(b)  The  second  and  subsequent  cards  each  define  the  data  point  from 
a specific  time  and  held  to  the  end  of  the  time  series  of  54  months.  This 
simplifies  data  definition. 

(c)  The  last  card  of  each  time  series  consists  of  the  key  word  END 
entered  in  columns  1-3.  An  example  of  a time  series  is  displayed  in  Fig.  4.4. 
The  name  on  the  last  data  card  entered  for  each  time  series,  names  the 
series,  a feature  which  makes  it  easy  to  change  a name  to  suppress  a file, 
for  example. 
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Fig.  4.4 — Xeroxed  Listing  of  Time  Series 
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Time  Series 

Definition  Cards 

(N  months  of 

data) 

Card 

Columns 

Entity 

Remarks 

First 

1-72 

Title 

Descriptive  title  of  time 

series 

2nd  and 

1-6 

Name 

Time  series  name 

following 

5-10 

Filler 

• 

11-15 

Month 

Month  of  data  - right  justified 

• 

16-20 

Year 

Year  of  data,  two  digits 
justified 

- right 

• 

21-30 

Data 

Numerical  value  of  data, 
point  must  be  present 

decimal 

Last 

1-3 

Name 

'END' 

Additive  seasonal  variables  are  stretched  periodically,  as  appropriate, 
if  the  first  four  characters  of  their  name  are  SEAS;  otherwise,  they  are 
stretched  using  the  above  sample  and  hold  technique. 

The  end  of  all  time  series  input  is  indicated  by  the  card  LAST  FILE, 
where  cols  1-4  - LAST  and  7-10  - FILE. 

Obtaining  a Data  Listing 

If  the  LIST  option  is  entered  on  the  first  data  card,  all  of  the  time 
series  are  tabulated  in  a form  suitable  for  Xerox. 

Selecting  Data  for  the  Regression 

This  is  done  by  specifying  the  names  of  the  time  series  desired,  in 
the  order  desired  with  additional  variable  control  information. 

BASE  Selection 

The  first  card  for  the  run  must  specify  the  first  and  last  data  points 
to  be  used  in  the  dependent  variable  time  series. 


Columns 

Entity 

Remarks 

1-' 

BASE 

may  be  blank 

5-10 

Filler 

11-15 

First  year 

Beginning  year  of  data  for 
analysis  (2  digits,  right 
justified) 

16-20 

First  month 

Beginning  month  of  data  for 
analysis  (right  justified) 

21-25 

Last  Year 

Last  year  of  data 

26-30 

Last  month 

Last  month  of  data 
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Variable  Selection 

This  is  performed  by  control  cards  specifying  the  file  names  in  the 

order  desired.  All  names  must  start  in  column  1 and  should  be  no  more 

than  6 characters  long. 

Example:  C123 

MILPAY 

CA-OPT 

BNS-CA 

This  will  set  up  a data  file  for  BMD  with  four  variables  as  specified, 
using  the  assumptions  that  the  first  variable  is  the  dependent  variable. 
Additional  runs  can  be  generated  by  further  "files"  of  control  cards, 
exactly  as  above,  in  which  each  variable  is  mentioned  explicitly.  In 
every  case,  the  data  used  will  extend  from  the  first  to  the  last  data 
point  specified  on  the  BASE  card. 


Repeated  Variable 

Repeated  variable  sequences  can  be  set  up  very  easily  for  use  in 
several  subsequent  runs  by  using  a SET  card  which  simply  preserves,  as 
base,  all  preceding  cards  in  that  file  of  control  cards. 

example:  C123 

MILPAY 

BNS-CA 

CA— OPT  run  1 

SET 

ADVT 

END 

PAID  TV 


RECR 

END 


This  sequence  will  set  up  three  five-variable  runs  in  which  the  first 
four  variables  are  always  the  same  but  in  which  the  last  variable  is 
changed  each  time.  The  point  of  this  is  to  reduce  the  number  of  cards 
required. 

The  SET  card  will  also  incorporate  additional  cards  in  a previous 
setting  which  allows  the  user  to  gradually  extend  his  base  case.  The 
setting  continues  until  either  a RESET  BASE  card  or  a RESTART  card  is 
encountered. 
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RESET  BASE  (Key  word) 

The  RESET  BASE  key  word  in  columns  1-10  clears  the  setting  of  SET 
and  returns  control  to  the  base  specification  point,  from  any  point  in 
the  control  card  sequence. 

Shifted  Time  Series  (Leads  and  Lags) 

MERGE  was  originally  constructed  to  permit  generation  of  multiple 
copies  of  time  series  with  various  shifts  so  that  regressions  with  leading 
and  lagged  data  could  be  performed  easily. 

If  the  shifts  are  listed  on  the  variable  cards,  they  are  generated 
automatically. 

Example:  MILPAY  -3  -2  -1 

will  include  four  variables,  all  copies  of  MILPAY  with  lags  of  -3,  -2, 

-1,  and  0 months  relative  to  the  dependent  variable.  0 is  always  included 
-nless  the  last  entry  is  99. 


example : 


MILPAY 


-2  99 


will  include  one  copy  of  MILPAY  with  a lag  of  2 months,  relative  to  the 
dependent  variable. 

Coefficient  Sign  Control 

Coefficient  sign  control  can  be  obtained  by  adding  a P,  N or  blank 
to  the  variable  card,  for  positive,  negative  or  free  sign  control 
respectively.  This  means  that  the  variable  and  any  of  its  copies 
generated  from  this  card  will  enter  the  regression  only  if  its  initial 
coefficient  sign  is  positive  or  negative  as  requested  on  the  card.  This 
is  done  by  checking  the  sign  of  the  partial  correlation  coefficient  of 
the  variolas  currently  out  of  the  regression.  This  overrides  CONDEL 
o«  :if icacion. 


- * "'jj'.u  jf  a variable  with  different  lags  in  the  regres- 

•- vers  lent  to  limit  the  number  which  can  actually 
" . t « joecifving  a limit  on  the  variable 

i ‘ . - • - - ' x.t  if  variables  -which  may 


I 


i 

i 

4 


Note  that  the  variables  are  entered  in  order  of  significance  at  that 
time,  so  that  3 of  10  does  not  give  necessarily  the  best  3 out  of  10  shifts, 
but  the  first  3,  Some  skill  is  still  necessary. 

This  restriction  overrides  the  CONDEL  card  specification. 

Table  4.10 

VARIABLE  CARD  FORMAT 


Columns 

Entity 

Remarks 

1-6 

Variable  name 

Left  justified 

7-10 

Filler 

11 

Initial  sign 

P = positive 

to  enter 

N = negative 
blank  = either 

12 

Variable  count 

Maximum  number  of  lagged 
variables  with  this  name 
allowed  to  enter  the 
equation 

14-15 

Shift  1 

Positive  of  negative 

16-20 

Shift  2 

Shifts  of  variable. 

21-25 

Shift  3 

99  restricts  unlagged 

• 

• 

• 

variable  from  equation 

61-65 

Shift  10 

Optional 

As  many  as  10  shifts  may  be  entered  on  one  card,  positive  or  negative, 
with  0 or  "no  lag"  added  on  the  end  unless  suppressed  with  a 99.  It  is 
the  user's  responsibility  to  ensure  that  shifted  data  is  still  within  the 
range  of  the  time  series  definition  by  adjusting  the  base  card  setting. 

STOP  Card  for  MERGE 

The  STOP  keyword  in  columns  1-4  ends  MERGE  so  that  BIOMED  can  be  run. 
At  this  state  BIOMED  receives  its  data  on  file  13,  its  control  cards  on 
file  20  and  certain  report  generator  data  on  file  11. 

DUMMY  (Additive  Seasonals) 

The  DUMMY  keyword  in  columns  1—5  serves  to  automatically  generate  12 
additive  seasonal  dummies  labelled  JAN  through  DEC.  MERGE  uses  a standard 
dummy  file  DMMY12  and  introduces  12  leads  of  0 through  11  to  generate  the 
months.  Thus  DUMMY  introduces  12  variables  and  the  C0NDEL  and  TRANGN  cards 
of  BMD02R  should  note  this  in  indexing  the  variables. 
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MERGE  Output 

MERGE  lists  the  cards  accepted  and,  if  the  LIST  option  is  taken,  lists 
the  data  retrieved  for  each  regression  run  setup. 

CREATE  - Creation  or  Modification  of  Variables 

The  key  words  CREATE  and  END  read  in  (A6)  format  signal  that  the  group 
of  cards  beginning  with  CREATE  and  ending  with  END  create  new  variables  or 
modify  an  old  variable  using  simple  scaling  of  variables.  This  is  equivalent 
to  transgeneration  cards  but  simpler  to  use.  The  formats  of  the  cards  are 
given  in  Table  4.11 

The  specific  options  offered  under  CREATE  are  detailed  below. 


Columns 

1-6 

7 

8 

9-15 

16-18 

19-24 

25-26 

27-33 

34-35 

36-41 

42-43 

44-50 

51-53 

54-59 

60-61 

62-68 

69-71 

72-77 

78-79 


Entity 

* 

New  name 


Remarks 

Name  of  variable  to  be  created 
Equality  sign 


Scale  factor 
* 

Variable  name 

Shift 

Scale  factor 
* 

Variable  name 

Shift 

Scale  factor 
* 

Variable  name 

Shift 

Scale  factor 
* 

Variable  name 

Shift 


Multiplicand  scaling  factor 

Indicates  multiplication 

Variable  scaled  by  preceding  multi- 
plicand 

Positive  or  negative  shift  of  value 
range  (-9,  +9) 


All  products  are  summed  to 
produce  the  new  name  variable 


If  new  name  has  not  occurred  before,  then  its  data  is  initialized  to 
zero,  a variable  of  name  "new  name"  is  created,  with  title  given  by  the 
creation  card  itself,  and  whose  data  content  is  defined  on  the  card. 

If  new  name  has  occurred  before,  then  the  modifications  are  applied 
to  that  time  series.  Thus  several  data  cards  can  be  used  to  define  a new 
time  series  or  modify  an  old  one. 
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Table  4.11 


FORMAT  FOR  CREATE  CARD 

Variable  scaling  and  linear  combinations  of  variables 


Columns 

Entity 

Remarks 

1-6 

Variable  name 

Variable  to  be  created  or  modified 

7 

Filler 

Blank 

8 

s 

Equals  sign  or  blank 

9-15 

Scale  1 

Scale  factor  for  first  variable. 
Decimal  point  and  sign  must  be 
included 

16 

* 

Multiplication  sign  or  blank 

17-22 

Variable  name 

Name  of  first  variable,  must  exist 
in  file  or  by  prior  CREATE  operation 

23-24 

Shift  1 

Shift  for  preceding  variable  range 
(-9,  +9) 

25-26 

Filler 

Blank  field 

27-33 

Scale  2 

As  per  scale  1 

34 

* 

35-40 

Variable  name  | 

41-42 

Shift  2 ' 

43-44 

Filler 

45-51 

Scale  3 

52 

* 

53-58 

Variable  name 

59-60 

Shift  3 * 

y As  above 

61-62 

Filler  1 

63-69 

Scale  4 ! 

70 

* 

71-76 

Variable  name 

77-78 

Shift  4 J 

79-80 

Filler  J 
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Table  4.11  (coat'd) 


Operations  on  CREATE  variable 


Columns 

Entity 

Remarks 

1-6 

Variable  name 

Names  of  variable,  predefined 

7-12 

Operation 

LIST,  list  variable 

CLEAR,  set  all  points  to  zero 

LOGN,  natural  log  of  variable 

Function 

in  CREATE 

Columns 

Entity 

Remarks 

1-6 

Variable  name 

Resultant  variable  name 

7 

Filler 

Blank 

8 

= 

Equals  sign  of  blank 

9-15 

Scale 

Scale  value  applied  to  function  result 

16 

* 

Multiplication  sign  or  blank 

17-24 

Function  Name 

Left  justified  - MAX,  MIN,  or  PRODUCT 

25-26 

( 

Left  parenthesis  or  blank 

27-33 

Scale  1 

Scale  value  applied  to  first  variable 

34 

Filler 

Blank 

35 

* 

Multiplication  sign  or  blank 

36 

Filler 

Blank 

37-42 

Variable  1 

First  variable  name 

43-44 

Shift  1 

Shift  of  variable  1 (-9,  +9) 

45-51 

Scale  2 

52 

Filler 

Blank 

53 

* 

Multiplication  sign  or  blank 

54 

Filler 

Blank 

55-60 

Variable  2 

61-62 

Shift  2 

63 

) 

Right  parenthesis 

64-80 

Filler 

Blank  field 
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Example:  (TEMP  = 1.74200  * MILPAY  - 1 + 1. 760Q  * CA-OPTS) 


which  creates  a new  variable  named  TEMP  from  MILPAY  (-  1)  and  CA-OPTS  (0) 
as  defined  above.  All  entries  are  right-justified. 

(a)  The  CREATE-END  packet  should  occur  immediately  after  a BASE 
card;  it  can  occur  anywhere,  but  must  occur  before  the  new  variable  or 
its  modification  is  referenced.  The  data  are  not  placed  in  the  BIOMED 
run  unless  explicitly  called. 

(b)  The  new  variables  have  all  the  properties  of,  and  can  be  called 
exactly  as,  permanent  variables.  It  is  anticipated  that  this  feature  will 
be  most  useful  for  generating  new  dependent  variables,  and  the  BIOMED  runs 
are  titled  with  the  mode  of  creation  of  the  variable  if  it  is  new  and  a 
blank  title  card  is  used. 

(c)  Constants.  If  the  variable  name  is  absent,  the  scale  value  is 
treated  as  a simple  constant. 

(d)  "TREND"  is  a reserved  variable  name  for  a ramp  function  taking 
values  1 to  54  in  months  1 to  54  in  CREATE  only. 

(e)  LIST.  The  reserved  name  LIST  in  cc  11  to  14  causes  a list  of 
the  variable  specified  using  the  report  format.  This  is  useful  for  spot 
checking  the  results  of  transformations  in  CREATE  only. 

Example:  C123  LIST 

Cf)  CLEAR.  The  reserved  name  CLEAR  in  cc  11  to  15  sets  the  specified 
old  named  variable  to  zero.  This  is  done  automatically  for  new  variables. 

Example:  TEMP  CLEAR  • 

(g)  Functions.  Simple  functions  are  permitted,  one  per  card  in  the 
following  format: 


Function  of  Variable 


Columns 

Entity 

Remarks 

1-6 

New  name 

Name  of  variable  to  be  created 

7-8 

= 

Equality  sign 

9-15 

Scale  factor 

Applied  to  result  of  function 

16-18 

* 

Multiplication  sign 

19-26 

Desired  function 

PRODUCT 

MAX 

MIN 

( 

Left  parenthesis 

27 

Scale  factor 
* 

Applied  to  following  variables 

Variable  name 

[ 

Shift 

> 

Comma  used  for  separation 

I 

Scale  factor 
& 

I 

Variable  name 

Shift 

) 

Right  parenthesis 

The  scale 

values  and  shifts 

can  be  omitted. 
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example: 

C123 

MILPAY 

BONUS 

END 

Run  1 

NEWDEP 

XNEW 

XNEW 

Run  2 is 

MILPAY 

END 

BONUS 

Caution: 

The  intended  primary  use  of  the  key  word  NEWDEPVAR  is  in  comparing  the 
effects  of  the  same  programs  and  policies  on  different  population  groups 
(dependent  variables).  When  using  the  special  multiplicative  mode  operation 
of  the  seasonal  model,  Section  4.9,  do  not  change  programs  and  policies  in 
conjunction  with  the  NEWDEPVAR  key  word. 

RESTART  (Key  word) 

The  key  word  RESTART  in  columns  1-8  returns  control  to  the 
BMD  CONTROL  LANGUAGE  at  read  in  point  PROBLEM  p.  4-13 . The  SET  setting 
(if  any)  is  also  cleared. 


4.9  BMD—  OPERATIONAL  MODES 

The  modified  regression  program  is  operated  in  four  different  modes, 
as  depicted  in  Table  4.12.  Note  that  operational  modes  3 and  4 require 
that  the  same  set  of  independent  variables  is  run  against  different  sets 
of  dependent  variables. 


Table  4.12 
OPERATIONAL  MODES 


Operational 

mode 

Special  controls 
required 

ddoRT  FM 

Special 

set-up 

instructions 

PROBLM  card 

CONDEL  card 

i ivuuLLil 

restrictions 

1 

Col  70  • blank 

none 

none 

none 

2 

Col  70  - M 

none 

none 

none 

3 

Col  70  - E 

Col  8 i 1 

Same  set  of  inde- 
pendent variables 

See  Note  1 
Section  4 . 5 

4 

Col  70  - E 

Col  8 - 1 

Same  set  of  inde- 
pendent variables 
with  the  second 
variable  having 
fixed  elasticity 

See  Note  2 
Section  4.5 

1 Additive  mode 

2 Multiplicative  seasonal  mode  with  no  restrictions  on 
problem  sets 

3 Special  multiplicative  seasonal  mode  with  no  fixed  elasticities 

4 Special  multiplicative  seasonal  mode  with  fixed  elasticities 
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Note  1:  The  MERGE  CONTROL  LANGUAGE  (key  word)  NEWDEPVAR  (see 

Section  III)  is  used  for  setting  up  mode  3 operation. 

Note  2:  The  special  instructions  required  to  set  up  a problem 

for  mode  4 operation  are  as  follows: 

Given  fixed  elasticities  e (k=l,...,q)  for  independent  variable  k 
k k 

with  mean  value  X , a new  variable  v^  is  defined  for  each  value  i of  the 
series  i using  the  CREATE  option  of  PROGRAM  MERGE  as: 

q e 

v “ 2 k i i = 1,...,N. 


v^  must  then  be  inserted  as  the  second  variable  of  the  regression  and  run 
together  with  the  remaining  p-q  variables. 


4.10  BMD  REPORTS 


Summary  Repore 

The  summary  report  summarizes  the  last  step  of  each  regression  and 
contains  the  following  information: 

a.  Dependent  variable 

b.  RSQ 

c.  Constant 

d.  List  of  all  independent  variables  used  in  a related  run  series  and 
their  time  shifts 

e.  The  coefficients  of  each  variable  in  the  regression  and  Fs  to 
remove  or  enter;  if  the  coefficient  is  non-zero,  the  variable  is 
in  the  regression  and  the  F-value  is  the  F to  remove;  if  the 
coefficient  is  zero  and  the  F-value  is  non-zero,  the  variable 

is  not  in  the  regression  and  the  F— value  is  the  F to  enter; 
if  the  coefficient  and  F-value  are  both  zero,  the  variable 
was  not  used  in  the  particular  run. 

It  should  be  noted  that  the  dependent  variable  is  assumed  to  be  the 
first  variable  and  unchanged  for  each  run  on  one  page.  Dependent  variable 
changes  are  precipitated  by  a RESTART  or  a RESET  Base  control  card. 

Residuals  Report 

A residual  is  computed  as 

A 

residual  »y-y  = y-  a - E /3.X 

i 1 1 

where 

y * dependent  variable 

y = a +L  /3.X.  is  the  prediction 
i 1 * 

a = constant 


1 


0^  = coefficient 
X * independent  variable 

i = number  of  independent  variables  in  the  regression 

b.  The  list  of  residuals  has  been  expanded  by  adding  future  data 
points  where  the  number  of  past  and  future  data  points  is  54.  The 
residual  computation  for  future  data  points  is  the  same  as  for  known  data 
points  except  for  the  setting  of  the  variable  to  zero. 

c.  Future  values  of  the  independent  variables  are  stored  by  automati- 
cally extending  the  last  value  of  each  time  series  to  the  full  54  points 
(Jan  70  to  Jun  74).  However: 

1.  If  future  policy  changes  are  known, these  can  be  added  to 
the  time  series  as  normal  data  points.  The  base  card  defines  the  cut-off 
between  past  and  future. 

2.  Seasonal  values  and  dummy  values  need  to  be  explicitly  extended. 

3.  The  PROBLM  card  is  automatically  revised  to  show  the  number 
of  future  points. 
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APPENDIX  A 


SETUP  INSTRUCTIONS  FOR  USAREC  UNIVAC  1108  COMPUTER  SYSTEM 

The  purpose  of  this  appendix  is  to  provide  instructions  on  utilization 
of  the  components  of  the  USAREC  Market  Information  Sybsystem. 

A.l  Optimal  Budget  Allocation  Model 

@RUN  RUNID,  ACCOUNT  NUMBER,  PROJECT,  TIME,  PAGES 
@TYPE  (this  card  is  required  by  systems  software) 

@PASSWD  Password  for  System 
@XQT  FILENAME.  OBAM 
Budget  Model  Input  Data 
@FIN 

A. 2 Optimal  Recruiter  Allocation  Model 

(§ RUN  RUNID,  ACCOUNT  NUMBER,  PROJECT,  TIME,  PAGES 
0TYPE 

0PASSWD  system  password 
@XQT  FILENAME.  ORAM 
Recruiter  Model  Input  Data 
@FIN 

A. 3 MERGE/BMD  Time  Series  Regression  Model 

@RUN  RUNID,  ACCOUNT  NUMBER,  PROJECT,  TIME,  PAGES 
@TYPE 

@PASSWD  system  password 
@XQT  FINENAME.  MERGE 
Time  Series  Data  (see  note) 

BMD  Control  Cards 
MERGE  Control  Cards 
@XQT  FILENAME.  BMD 
@FIN 

If  the  user  wishes  to  save  any  of  the  output  on  intermediate  results 
that  are  contained  in  the  work  files,  each  file  desired  must  be  explicitly 
created. 

MERGE  files  - the  following  cards  must  precede  the  @XQT  FILENAME.  MERGE 


card:  @USE  UNIT//,  @ASG,CP  NEWFILE,  F/l/TRK/64,  NEWFILE  where  NEWFILE 


1 


- 


I 


is  a user  file  name  in  correspondence  with  installation  standards,  and 
unit  # is  one  of  the  unit  numbers  presented  in  Table  4.2. 

BMD  files  - the  following  cards  must  precede  the  @XQT  FILENAME. 
BMD  card:  @USE  UNIT#,  @ASG,  CP  NEWFILE,  F/l/TRK/64,  NEWFILE  in  this 

case  unit  # corresponds  to  Table  4.3. 

Note:  The  time  series  data  could  be  contained  on  a card  image  mass 

storage  file  and  could  be  inserted  at  the  appropriate  location 
in  the  runstream  as  follows: 

@ADD  TSDATA. 

where  TSDATA  is  an  example  of  the  file  name. 

If  multiple  sets  of  time  series  exist,  they  could  be 
included  in  one  file  as  symbolic  elements.  Your  system 
software  should  be  able  to  clarify  procedures  for  creating 
this  type  of  file. 
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• MICROCOPY  RESOLUTION  TEST  CHART 
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SYSTEM  SPECIFICATIONS  FOR  USAREC  1108  COMPUTER  SYSTEM 

The  contents  of  this  appendix  present  the  system  resource  require- 
ments of  the  different  components  of  the  USAREC  Market  Analysis  Subsystem. 

In  addition,  the  subroutines  and  functions  used  in  each  component  are 
presented  with  a short  description. 

B.l  Optimal  Budget  Allocation  Model 
Resource  requirements 

Average  run  time  - minimum  8 CPU  seconds 
Core  required  - 12K  words 
Input  cards  - minimum  5 
Output  pages  - minimum  10 

Subroutines  (presently  compiled  as  one  program  with  internal  subroutines) 
MAIN  program  - the  main  program  reads  the  input  cards  and 
verifies  parameter  names. 

SEARCH  - performs  general  fibonacci  minimum  search. 

FVAL  - function  of  Newton  method  to  solve  for  budget. 

B.2  Optimal  Recruiter  Allocation  Model 
Resource  requirements 

Average  run  time  - approximately  50  seconds  CPU 
Core  required  35K  words 
Input  cards  - approximately  600 
Output  pages  - approximately  60 
Subprograms 

Main  program  - the  main  program  serves  as  a monitor  for  the 
execution  of  the  model  options  and  prints  solution. 

READ IN  - this  subroutine  reads  input  data,  preforms  preliminary 
calculations  and  performs  minor  error  checks. 

READOU  - this  subroutine  produces  the  evaluation  reports. 
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REQMTS  - calculates  adjusted  right  side  value  on  constraints 
to  account  for  recruiter  and/or  canvassers. 

START  - calculates  the  initial  value  of  Lagrange  Multiplier 
for  setting  Search  Interval  end  points. 

PARTSL  - calculates  canvasser  partial  solution. 

SEARCH  - performs  general  fibonacci  minimum  search. 

FVAL  - function  of  Newton  method  to  solve  iteratively  for 
recruiters  or  canvassers. 

SOLN  - subroutine  that  determines  final  solution  using  the 
Newton  method  as  in  the  function  FVAL. 

ENLST2  - subroutine  that  calculates  enlistments  by  cohort  group 

B.3  MERGE/BMD  Time  Series  Regression  Model 
Resource  requirements 

Average  run  time  - minimum  8 CPU  seconds 
Core  required 

MERGE  20K  words 
BMD  39K  words 

Input  cards  - including  time-series  data  usually  200-600 
Output  pages  - minimum  16  pages 
Subprograms 
MERGE 

Main  program  - reads  all  input  data,  BMD  controls  cards,  and 
MERGE  control  cards.  Performs  syntax  check  on  all  control 
cards. 

CREATE  - creates  all  new  user  defined  variables. 

MONTH  - sets  up  alph  name  abbreviation  of  months  for  use 
by  BMD  reports. 

NOTE  - write  output  message. 

FILES  - reads  time  series  data  check  into  core. 

IDENT  - verifies  that  all  variables  to  be  used  in  regression 
exist. 

PACK  - left  justifies  name  and  removes  embedded  blanks. 

i 
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FILE50  - creates  edited  control  card  file  for  input  to  BMD. 
LISTRP  - lists  time  series  data. 


Main  program  - the  main  program  reads  control  cards,  checks 
syntax,  and  performs  calculations  for  the  basic  statistics 
(mean,  standard  deviations,  correlations) 
function  AF  - orders  array  elements  by  index  pairs, 
subroutine  AOUT  - prints  correlation  or  variance  matrices, 
subroutine  CHECK  - checks  sign  of  coefficients, 
subroutine  RDLBL2  - reads  in  labels  cards, 
subroutine  TRANGN  - variable  transgeneration  program, 
function  CALUIN  - zero  divide, 
function  LSHIFT  - set  character  values. 

subroutine  RESIDS  - performs  residual  analysis  and  plotting, 
subroutine  STEP  - enter  variable  into  regression, 
subroutine  STEPRG  - performs  stepwise  regression, 
subroutine  REGRPT  - prints  stepwise  regression  report, 
subroutine  SENSE  - evaluates  sensitivity  of  coefficients  for 
ridge  regession  and  prints  report. 

subroutine  MULSEA  - computes  the  multiplicative  seasonals,  and 
prints  seasonal  report. 

subroutine  REPRT2  - report  produced  when  fixed  elasticities  are 
indicated . 

subroutine  REPORT  - summarizes  the  end  results  for  each 
regression  run. 

subroutine  MOVECH  - moves  characters  from  a character  string  to  a 
second  character  string, 
function  COMP  - word  comparison. 
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THE  HISTORICAL  ACCESSION  FILE 


The  Historical  Accession  File  was  developed  from  USAREC  accession 
data  for  FY74  and  FY75.  The  file  contains  the  fraction  of  total  acces- 
sions accounted  for  by  each  of  the  supply-limited  cohorts  under  con- 
sideration, as  depicted  in  Figure  C.l  on  the  following  page. 

Table  C.l  displays  the  current  definition  of  items  in  the  sort-field. 


Table  C.l 

Entity 

Columns 

Definition 

DRC 

1-2 

DRC  identifier 

DRC  IND 

6-7 

Sequential  DRC  index  (1-64) 

REGION 

10 

Sequential  region  index  (1-5) 

ENL  OPT 

11 

Enlistment  option  coded 

1 = combat  arms 

2 = non-combat  arms 

TERM  SER 

12 

Term  of  service  coded 

1 = 3-year  term 

2 = 4-year  term  or  more 

Cohort  Group  1 

13-17 

Proportion  of  enlistments  in 
first  cohort  group. 

2 

18-22 

3 

23-27 

4 

28-32 

5 

33-37 

6 

38-42 

7 

43-47 

8 

48-52 

Filler 

53-57 

TOTAL  ACCESSIONS 

58-67 

C-l 
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ANNOTATED  LISTING  OF  THE  OPTIMAL  RECRUITER 
ALLOCATION  MODEL 
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ANNOTATION 


SUBROUTINE 

STATEMENT 

NUMBER 


14 

15 

16 


20,21 

22 

24 

25,26,27 

28,29 

31,32 

33 

37 

38,39,40,4 

42 

44,45 

46 

47-51 

52 

53 

54,55 

56-59 


ROUTINE:  CONTRL 


This  is  the  central  control  unit  for  the  ORAM 

The  call  to  READIN  reads  all  the  input  data  and  switches  for 
the  problem. 

READOU(O)  gives  a full  evaluation  of  the  initial  recruiter/ 
canvasser  allocation. 

If  IPUN  J 1 then  only  an  evaluation  of  a particular  allo- 
cation is  made — otherwise  a full  set  of  optimizations  is 
made. 

ISWITCH  = 2 for  recruiter  and  canvasser  model  set  upper 
limit  for  canvassers. 

Set  constraint  limits  and  adjust  coefficients  according  to 
constraint  no.  M7. 

Calculate  guess  for  Lagrange  multiplier  w 
Branch  to  250  for  recruiter  only. 

Calculate  initial  guess  for  canvasser  multiplier;  set  end 
points  for  Fibonacci  Search. 

Search  and  find  candidate  value  of  multiplier  for  canvassers. 

Find  tentative  solutions  for  canvassers;  adjust  right  side 
constraint  value  for  recruiters. 

Calculate  candidate  multiplier  value  for  recruiters. 

If  canvassers  achieve  more  than  total  accession 
requirements  recalculate  recruiter  solution. 

L Calculate  candidate  value  for  multiplier  test. 

Calculate  canvasser  test  multiplier  for  final  decision. 

For  positive  canvasser  multiplier  set  associated  value  of 
recruiter  multiplier 

Canvasser  multiplier  set  to  zero. 

Set  final  recruiter  multiplier. 

Generate  solutions  for  recruiters  if  ISWITCH  = 1 or 
recruiters  and  canvassers  if  ISWITCH  = 2. 

Loop  end  for  all  constraints. 

Write  out  full  set  of  solutions. 

Write  out  matrix  for  recruiter  solutions  only;  not 
intended  to  be  a report  as  such. 
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ft'  5 E COi 


program  cuntrl 


7 3/73 


OPT  * 1 


FT  *4  4.  ?*h:; 


:j/i:/v  n.-.p.3’ 


PAGE 


3 3 


P^r:,r  A k CCNT?L  l INPUT  .OUTPUT  , PUNCH?  TAP*  CsINPUT  , TAP £6= CUTOUT , TAP £7 *P 
lU'.'CH) 

c . uf!r*'  /c  ize/:m<sta  , j ty  >ts , •;  r lrms.nor  egn 

wGM:';‘:/C«LCCr  / £ (/,  ')  , 0 (Ev  ,3)  ,R£C  (t5)  »CAN<65)  , cNL  C65?4) 

CC-MC*./ P-P.C'-r  / PC  T (65  « *< ) * P C T T r.  ( 1. 1>  * d » 3 1 ,p:'lTCA(  65,  3,21 
coim:  ./OkC  (.:•./  j wgrcio5> 

CO'*  -tCN/CONS  WI  : , 12, 13 , 14, 15,1 SWTCH, IR6GSW 
COlMCN/^TSICw/  R3IO£(50) 

c ./  *MCC  / I 'J  , I £ , I C CUN  T . Fl A G*  «( 5C  » , »S I9£R( 5 0 ) . RSI OEC 1 5 C ) 

CCMM  "N/ RESULT  / XR£C  (w>  , 15)  , X CAN  (5  5, 1 5)  , RECG  ( o5  ) , CA  N j ( 65  ) , RLAGRG  (50) 

1,  Cl  Af*  G(5C) 

co,,mcu/cos  t/f:  *cst  ,R£ccsr  ,c  angst 
ct»ms.u/qut2/:p.ln 

CALL  R £ A 0 1 N 
CUL  P £ A 0 0 U ( C ) 

:F(  IRu:  !.:»£•  i ) CALL  E * I T 
CO  1 G J 0 K T 5 1 , 1 5 

IF  ( ( 2SWTC)!.::i.i  ) .AND.  ( MT.E3.2  ) ) GO  TO  1C00 
Ir(  IOrtTCH,*.,  ! ) CALL  REQMTS(2,i) 

C AwL  : £ C t s ( ’ , 1 ) 

CAll  CCMTG  (NT  ,2) 

IF< ICw TCH, £0. 1)  call  STARTR(WW) 

» 3RANCM  IF  kcCP’J  I T£R  j Aftt  TH£  CNl  Y VARI AQLcS  •** 

IF  ( I i*»  T CH#  EG  • 1 ) 30  TO  25C 
CALL  STARTCIVVM 
TL  - / . / 1C, 

tu  * i : . • vv 

CALL  ''iArCH(^L,TJ,l.t-2,1.0,FFVAL»TT4,l,«T) 

V J - f T A 

■ CALCULATE  CwK’ENT  CANVASSER  PA  FT  IA  L SOLUTION 
CALl  *>•.*>  TSL  ( ■•'/) 

CAwL  f-£  TS  ("T  ,1) 

CAlL  3 T A r T R ( W W > 

► Cj-.REWT  CANV ACC£R  L.;cL  OiT&RMlNgS  RIGHT  SIQ£  LIMITS  ON  RECRUITERS  •* 


50 


IF  (RSI  JEM  f*T)  •L£,3«3)  GO  TO  250 
TU  = 10.  * WW 

TL  = WW  / 10. 

CALL  SEARCH  (TL  , TU»1.  £-2,1.0,  FF  VALUTA,  2,  MT) 

WW  = T T A 

CL  AG-.G  ( MT)  = VV  * RECCST  / WW  - CAUCST 
C •*  test  FOR  *:Cf:-%.£GATlVlTY  CF  CANVASSER  LAGRANGE  MULTIPLIER 
IF(CLACPGC'T)  .ST,  0.0)  RLAGRGCMT)  a RECCST  ✓ WW 
IF<C~AGRG<mT) ,GT. 0.0 ) GO  TO  500 
250  CLAGAC(MT)  = C. 

WW  = Pf CCS T / WW 
TL  = W W / 10, 

TU  * 10.  • WW 

CALL  S*AFCH(TL, TU, t.fe-2, 1.0,FFVAL.TTA,3,MT) 

FLAGRu(MT)  a T T A 
500  CALL  SCLN(MT) 

1C  3 0 CONTINUE 

:o  looc  hc«i?I5 

‘♦CCS  CALL  RSAJCUtKC) 

w-:T£;e , i: so  d,i=i,i5) 

103  3 FORMA  T < lrtl , 3’X  , • SCLUT  1 0 fS»  > 13 1 1 1/  I 


jc  z:;a  i=i,raxsra 

zz:o  h-iTKt,  i:  j-on^cdi , i,-.£c:  a i,  (xseci  i,mt»  xs> 

uii  FC’rar(i<,i»,i2,i3<it,Fio.i)) 
caLL  exit 
in  o 
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ANNOTATION 


SUBROUTINE 

STATEMENT 

NUMBER  SUBROUTINE:  ENLIST 


This  subroutine  calculates  the  enlistments  by  cohort 
type  j at  DRC  i for  a given  number  of  recruiters 
and  canvassers  on  station. 

8,9  Reset  recruiters  to  X(i)  value. 

10-12  Calculate  enlistments  due  to  recruiters. 

14-18  Calculate  additional  enlistments  due  to  canvassers 

if  ISWTCH=2. 
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L 


ocq  m : E COPY 

U J l « • • • 


w »•  ^ - i T . ‘-i  l 'L  CK 


/c/7?  C P T = £ 


F T N 4«y*R4*G 


C2/10///  ;9.45.J2 


PAGE 


13 


23 


33 


43 


45 


SlM°JtlT!f.:  “’LOTCD'C) 

COVM  n/P*  3ULT/  x C ( c 5 t 15)  . y c A *4  (65  • 15  I , c£.CJ  <o>  ) . GANG  (65)  • RLAG3G  <5C  > 


1 .31  t >PCTCM  65,9,2) 


L « C - A ij  ■' 

: ( c ) 

~ " *4  0 : 

/ODT/MA  fTl  ,J 

T»  «£5,  <Ti«r2;  ,N«E 

0 " *■*.  I*. 

/ : c If-'/ pct  (t>5 

, 3 ) , A 0 T T v ( O.'l.J) 

CCVVC  . 

/.'  eGi  V./  I--.*:  3.(5)  , 13  uGI M *3)  , UNO 

CC*< ~ V. 

/c  ' 

CO'-MC  . 

/ZilC"':/  : c / , 5) 

• C(o?i3)  f 4 i o 5 ) » i 

rcc  ‘ 

VO  - 14 

If  ( MC 

C. cQ.il  GO  TO 

103 

30  113 

1*1  . '-AX  OTA 

.<*  < : > 

= x~;c(i#mo 

:-i  re 

. 2J.:  ) -.jo 

= RECOU) 

00  no 

J=1 . JTYPES 

113 

E.-iL  ( : . 

j)  = ct r , j>  • 

RECCI) * • £ ( 6 > J) 

go  re 

115 

13  0 

ro  o o 

l»  1 •NURSGN 

LL  = L 

- 1 

1 J 1 = 

i’Eoritu 

1 “3  - 

:g\’3(l> 

CO  It 

Is ! 3 1 tT  ;■ 

f o : 1 1 > 

= X R£C  ( I»  f 'C  * L 

L) 

If  ( *M 

.03, : > ‘•oco 

= RtCC(I) 

00  iv 

)*l,  ^ryocs 

15 

o ( : , 

J > - CC  , J)  * 

^£0(1)  **  c ( 6 , J) 

■> 

c : : r : 

•c 

* » * 

I-  { ; 

wtch.nc.2  ) go 

TO  150 

: - ( 

29 

no  u 

1= If  VAX  OTA 

can  ( : ) 

= XCANt I,; C) 

!*■*<  vc 

.no,:  ) canid 

= CAM}  (I) 

co  12: 

J = ! • JT  YC£3 

123 

£\„  < : , 

J)  * £NL  ( I , J) 

C ( I • J ) • CANII) 

co  re 

150 

29 

00  4 3 

L = tf NOPSGN 

L L - L 

- 1 

122  = 

I 9 £ G I N ( L ) 

I “ r = 

ic?:C  (L) 

CO  j 3 

1*133, IdE 

C A ‘1  C I ) 

= XC  AM  ( I,  VC*  LL) 

Ir  ( VC 

. £G, 0 ) CAN ( £ ) 

= CANom 

00  3 L 

J * 1 * j types 

3C 

L'*«C  (I  » 

J)  = ENLC.J) 

♦ C(I»J)  • CAN (I) 

4C 

CONTINUE 

ISO 

RETURN 
Ef  10 

SYSTOLIC  R£F£R£NCC  MAP  (R=l) 


ENTRY  POINTS 
J EN4ST2 
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L 


This  subroutine  is  used  only  when  canvassers  are  in 
the  model.  For  specified  values  of  W and  a given 
constraint  index  MT  the  associated  canvasser  solution 


is  calculated.  The  sequence  of  W values  assigned 
by  the  Fibonacci  search  is  such  that  convergence  to 
the  requirement  specified  by  RSIDEC  (MT)  occurs. 

30  Error  message  if  fail  to  converge. 


SEVERITY 


I 


'/r«»  v “ ' / £ ( / « * > .0(65,3)  , RE C ( 65  ) , C AN  (o5),£N;.(35»d} 

c ::»icn/;  .a /w.\x  c:  *. , j it  p.e , < r :r.*s  ,ncpesn 
co-mo*i/cc:t/f  z>c;i  .reccst  ,cancst 

rsicc/  r::o:<5:) 

.✓*.£. :*\:  / : Jf :e , ::  :jnt,flagrm<5C> , fs:der(  su»  » rs:oec<53> 

Co  ■ :*j/? eto/i  ('5,5) 

: •**  to  solve  iterative.,*  fcr  recruiters  - given  lagrange 

C MULTIPLIER  .CLA.V3A  **• 

cPSI  = l.fc-3 
SUH  =0. 

XL  0 = I. 

co  l;  :=id,ie 

5 XL  = X L 3 

cc  z:  L=ifro 
su  ■:  = c.i 

SUM  2 z o.C 
SC  10  J=1,JTVP£S 

Tl-n  = £(7,  J)  *5(1  ,J)  *XL**  (£17,  j)  -l.» 

SJ*'i  = SU'-’I  ♦ TER  ’1 

1C  SL'^2  = CL ”2  ( £ ( / , J ) - i.)  * TERM  / XL 

I?  ( COM2  . £ G.  <js  ) / Ll  = E 3S I 

IF f SUM2 . EG . I . > Gw  TO  2 5 

XcL  = XL  - ( S'J VI  - 77  ) / C»JM2 

if! xll .lt. : .)  al : = xlo  / ? . 

IF(XLL.LT#u.)  GO  TC  5 
r:-L  = A«3<  XL l - XL) 

IF(  T LL.LT.EPSI  ) GO  TO  25 
23  Xw  = X-L 

W.R 1 T £ ( 6 , 52 ) I , X LL , XL • V 7 

5»2  FO-v  AT  W/i  :X,*ERRD>  MESSAGE  - FAILURE  TO  CONVERSE  IN  FVAL1  1=  *, 
1I3»*  VALUE  1 = *,£10.5,*  VALUE  2 = *,£10.5/2: X, *LAMCA  = *,£10,5) 
CALL  EXIT 

2 5 SU-  = SUM  ♦ XLL 
40  CONTINUE 

F V A L 1 = A33<SUM  - RSIJEC<MT>) 

K £ TU  K N 
ENO 


DETAILS  O I AG NOS I S OF  PR03LEM 


25  CO  32 
L 3 CO  32 
69  CO  32 


field 

field 

field 


WICTH  of  a 
WIDTH  CF  A 

wICTh  of  a 


CCNVERSIO N 
CONVERSION 
CONVERSION 


DESCRIPTOR 
CESCRI  FTOR 
CESCRI f TC  ^ 


SHOULD  CE 
SHOULD  BE 
SHOU-O  3E 


AS  LARGE  AS 
AS  L A ROE  AS 

AS  large  AS 


the  MINIMUM 
the  minimum 
Th£  MINIMUM 


s=ec:c:eo  fc.r  t-i: 
specified  for  That 

SPECIFIED  FCr.  That 


Y M 3 J l I C REFERENCE  map  (R=  1 ) 


ANNOTATION 


SUBROUTINE 

STATEMENT 

NUMBER 


30 


SUBROUTINE:  FVAL2 


This  subroutine  is  used  only  when  canvassers  are 
in  the  model.  For  every  value  of  WW  assigned  by 
the  Fibonacci  Search  the  corresponding  number  of 
recruiters  associated  with  constraint  (MT)  is 
evaluated.  Successive  values  WW  assigned  by  the 
search  will  be  such  that  the  return  gained  by  the 
recruiters  will  converge  to  the  adjusted  right  side 
RSIDER  (MT) 


Error  message  if  fail  to  converge 


i fval ? ?f/75  opt  -1 


FTN  4,5*R4C6  02/10/77  C9.45.A2  PAGE 


r .ctic*. 

u N A l 0 " / . ( /'  • 3 > , C ( <-  J v EC  ( 6 5 > . 0 A M G5 ) » i >i L i 05  • 3 ) 

C *•’ " V : / S ! 2 : / * - T A , J r y " - ; • < T ER  FS  • M OR£G  N 

z o • • ••  c ; / . ' w t / r : x c ; r , . : c is  r , c a -a:  st 
c :\/rTs:  )•/  rsiokso) 

OH'*:  :'/  I3,IS,:CCUNT,FLAGRM(5C)  ,RS:C£R<53)  t RSI  SEC  (5  C ) 

c '.'—  l-./ZZX  */£(S5*ai 

: •••  \-:»i „on  l t - ; : to  solve  iteratively  for  recruiters  - given  lagrange 

0 * • * f .•  w I Z 7 F L A M3  A *** 

2 ° S I = i.t-3 

:j“  = c. 

Xtfl  * 19, 

CO  40  1 = 19,  IS 
5 XL  = XL  3 
:c  20  l = i • o a 
Cv’U  = 0 • c 
So*'*  = G.) 

*0  12  J= 1 J TYPES 

TE'M  = £ (6,^)* 6(1,3)  *X  L**  (E  (6,J)-1.) 

S’jmi  = SLM  1 ♦ rE-** 

12  ZV’C  = C U * * ? ♦ (£C5,J>  - 1.)  • TERM  / XL 

/u.  = :3sr 

<sj  *2.=.:.c.)  go  re  25 

>LL  = > L - 1 Su-i  - ww  ) / S’JM 2 
I r ( < L W . L T . 2 . ) X u 3 - x L Q / 2. 

1 r < > . - . L T . ) GO  TO  5 

TO*.  = ATS<  <LL  - XL) 

IFJTOi  .L r»£*m  ) GO  TQ  25 
2:  >-L  = x LL 

w-ITE(6,S2>  :,y-L,XL,Fi.AMOA 

52  FIaMAT (//10V, 'ERROR  MESSAGE  - FAILURE  TO  CONVERGE  IN  FVAL2  X=  *, 

1 1 3 , * Vijji  l = * ,ElO  • 5*  ♦ VALUE  2 = * , E :0 . 5/22  X , *l  AMO  A - *,£13.5) 

C A L L EXIT 

25  :o  j;  j=i,jty?£S 

TERM  = 2(1, J)  * XLL  *»  E(6,J) 

z : sum  = sum  ♦ term 
4C  CONTINUE 

F V A L 2 - A6  S ( SUM  - RSIOER(MT)) 

KcTJRM 

cNO 


SEVERITY  DETAILS  DIAGNOSIS  OF  PR03LSM 

l 23  CO  32  FIELD  rfJOTF  OF  A CONVERSION  OESC=IPTOR  SHOULD  3E  AS  LARGE  AS  THE  MINIMUM  SPECIFIED  FOR  T^A T 
I 43  CO  32  FIELD  WI^TM  CF  A CONVERSION  CESCRIFTOR  SHOULD  3-  AS  LARGE  AS  Tit  MINIMUM  SPECIFIED  FCR  T-»AT 
I 69  CD  32  FIE LU  hIOTH  CF  A CONVERSION  DESCRIPTOR  SHOULD  9E  AS  LARGE  AS  Th£  .MINIMUM  SPECiFIc.0  FC.R  THAT 


SYMBOLIC  REFERENCE  MAP  ( R= 1 ) 


ANNOTATION 


SUBROUTINE 

STATEMENT 

NUMBER  SUBROUTINE: 


FVAL3 


This  subroutine  uses  the  Newton  Method  to 
solve  iteratively  for  recruiters. 


71/7  T Cr>T=l 


/ 1 0 / / 7 


TA  . JT  Yr> 


R<  5 J ) .RSIOEC  (50  ) 


YEN  LAGRAN 


sccsr  / flamo 


j-  :♦  j types 

£(<»J)  * H « I , J ) 'XL**  C (<i J) -it) 


1 : SUM  2 = SUM’  ♦ (Z(<.J>  - 1.)  * TERM  / XL 

I'(  SJ'«2,c3.0.)  Y L L = £ PS  I 
If  (Sc---2.ES.  :.)  CO  TO  2 5 
XL.  * <c  - ( SO* t - V AL  ) / SUM2 

:r(xLL.ur.:.i  >l:  = xlc  / 2. 

IC(XLL.LT.3.>  GO  TO  5 
T Cl  = A v S ( * l L - XL) 

If (T0L.LT.ZP3I  ) GO  TO  25 
22  XL  = x LL 

w:TE<t,52>  r,X  LL.Xt.  FLA  MO  a 

52  FO  VA  *(//::  X,*cftf.OR  MESSAGE  - FAILURE  TO  CONVERGE  IN  FVAL3  I 
1 1 3 * * VALUE  l = *#£10.5,*  VALUE  2 = * ,E  10.5/2: X , *L AMOA  * *,£l 
CALL  exit 

25  SO  3:  J=I.JTYPES 

TERM  = O(ItJ)  * XLL  •*  c ( K t J ) 

3C  SUN  = SUN  TERM 
40  CONTINUE 

XF<  riWTCH  • £Q • 1 ) GO  TO  50 
If (LSET.EO.l)  GO  TO  50 


GO  TO  4 

50  rVALJ  = AOS*  SUM  - KSIOcCMTJ) 
KETURN 
ENO 


ANNOTATION 


SUBROUTINE 

STATEMENT 

NUMBER  SUBROUTINE:  RE ADI N 


This  subroutine  reads  the  input  and  switch  values 
for  the  ORAM 


39,40  These  values  are  default.  Currently  set  for  5 different 
regions  NOREGN.  IREGSWjU  was  dropped  when  it  was  decided 
not  to  specialize  by  particular  cohort  type. 

41-48  Initial  environmental  factors  set  to  one  in  case  they  are 
raised  to  a power  but  not  read  in. 

49  READ  and  write  input  cards. 

51-53  Free  title  card.  This  title  will  appear  on  all  reports 
Use  it  Wisely. 

56  If  IRUN=1  full  optimization,  otherwise  evaluation  of  initial 

recruiter/canvasser  placement  only. 

59  Set  size  on  no.  of  DRCs,  no.  of  supply  limited  type,  no. 

different  terms  of  service. 
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MJ!  * 


f 





; 


ggj  ** 


SJ«*R0j7llJE  RiiOIK  72/71  OPT=t 


FTN  4.5*R436  0 2/  13 ///  G5.-,5.32  PAGE 


13 


23 


33 


35 


45 


53 


SL'iR  ) j T : :. r ua:;im 

C w,m:  i/c»-  - * :2  / I (13  3)  , N ,M,  H.)  1 , N'H 

C ' 'M  • /ST  A t V 0 <cv.  3)  , a <«,*>  , -O  , J d6  , i ) , A T ( o5 , 6 ) , 40  ( 65,  5) 

C2*:Hi  ./P>.:  ~n/a  TT  (65,  i)  ,5J  TM  (0i  , * , 3)  , KPCTCA<  65, 3, 2) 

c v.MC*:/?£f, rr‘j/  r-  gmc •; i , ? , i-nsc  51 

cj  *‘«:n/c:  :c fR/*i ♦ ic ,13 ,14 , i>  , iistcm,i  ;egsr 

;/:r^2-T/~CTv;it5^,-i,  Jl  ,PCT:K65,»*,2)  , NCAPC  T (65 

1 , J ) , C A PC  T ( 65 • i ) , SC  TS  T A ( 65 » 

c^'mc  ,/ciw sc*./;  (7,3) , : ( t; , j. ) , - :c  (o-; ) .can ( 65) , enl  ( 65,  a) 

CO  :*1  . : /SIES/^AxCT  *• , J TYPE$,<Tc.RKS,NOPSGN 
C ' ^-w../3UP3LY/  C[CN(fl  ) 

COM  1C  :/CCj  Y /Fixes  T,  RECCST.CANCST 
CC-:iC-,/-.TSiri/  RSIC£(50) 

CC*'’  ~ = f I£,ICCU\Y  ,Fw4G;K52>  ,RSI0SR<  50)  , P.siDcCC  50) 

C C •■*. « C ■ . / ACC  £ SS/  4 C Z 3 £ C ( o 5 , 4 , 1 5 ) , ACCC A N ( 6 5 » t , 1 5 ) 

C 2:'  * "/3-Cjl^/*“CC(6  5»15)  , X C •*•  N (66  * 1 5 ) «R£CQ(6?)  ,CAN0(65),RLAG*G(50) 

1 ,CL  A SR  :<  53  > 

co-*. ont/  fractt  ( 5 c i ,ca.nfra<5g> 

CCmmCN/J-\0  r :t/NM0?Cf  b6  ) 

C 0NM0.«/''lT/.\;hc-U  ( 3)  * ITITl£C13) 

CCOM  :/C'jT  I/IRUN 
integer  r.NAM'di) 

DIMENSION  valcm 
o; w£t:s :cm  cc;. ( 3> 

:c?:  nr2‘c(65> 

0:MioIC.N  p "i  GN  ( 1 2 ) ,FSE°V<5) 

real  ‘ Capct 

c : ■■  £ r :gm  signs c *) 

cat  a signs/:, . l.  #-i. • -i.t-i • 

CAT*.  I NAME /5>*;  NO  ,5hQvE  , j ->W  AGE  A , 5 H .JNj-MP , pH  A IT  ,5HACV  , 3 H p.  £ 0 R 
1 • SmCAN'/  »;'■£-  A 3 T , 5 n CONST  , 5h?CT  / 

o s ° . c : a l switches  — iswtch  anj  regs* 

C *•  IF  IS-'CH  l The  l mCSEL  IS  "OR  R£0  A.  j I T£R  S OfiLY 

C *•  IF  ISwTCH  NOT  EC.  1 Tm£N  l T£  R3  AN  3 CANVA;,SESS  A~£  VARIABLES  M* 

C***  IF  IP'fGSW  22.  1 Ti*H  REGIONAL  COh3TRAINTb  ARE  TOTALS  ONLY  • •• 

CM*  IF  IRSOTW  "0.  2 TntN  REGIONAL  CONSTRAINTS  ARE  3Y  COHORT  TYPE 
I***  IF  IR’SSW  £).  3 THEN  TH  £R£  ARE  NO  REGIONAL  CONSTRAINTS 
C***  IF  I-ESSW  NOT  EC,  l THEN  REGIONAL  CONSTRAINTS  ARE  3Y  COHORT  TYPES  »*• 
1REG5  w - 1 
NCRESN  = 5 

C *•*  IllirlALlCc  £NV  IRONHFNT  al  factors 
00  31  1=1,65 
00  31  J = 1 • 9 
Q(I*J)  = 1. 

U(I,J)  = l. 
w (I ,J)  = 1 , 

AT(t,J)  = 1. 

31  AO(ItJ)  = 1. 

WRITE (6,1^) 

14  FORMAT (1H1//19 <•-•) ,*IN?UT  CARCS  AS  RECEIVED* //> 

RE A 0 ( 5 ,28)  ITITLE 
23  FORMAT (3A13/5A10) 

W3ITE (6,29)  ITITLE 
29  FORMAT (9X, 3A1C/9X.5A 10 ) 

Z FORMA T ( 1 C 1 5 ) 

RE  AO  ( 5 » 2)  IRUN 
WRITE (6, 42)  IRUN 
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SUBROUTINE 

STATEMENT 

NUMBER 


64 


67 


71 

75 

79 

86 

89 

94,95 

96 

97 

98,99 


1 


ANNOTATION 


SUBROUTINE:  READIN 


ISWITCH=1  recruiters  only,  2 recruiters  and  canvassers, 

ISWPCT=1  proportions  by  type  set  by  historical  percentages 
otherwise  ISWPCT^l  will  accept  coefficients  individually 
by  type.  (The  latter  is  generally  not  the  case.) 

FUNEMP  is  a factor  (near  1 e.g.,  1.5)  which  modifies  the 
basic  unemployment  rates  at  which  the  x sectional  analyses 
was  made.  FACTOR  multiplies  the  coefficients  for  missing 
data — e.g.,  given  8 months  of  data  for  analysis  multiply 
by  1.5  to  get  yearly  total. 

SUPLY  is  desired  level  of  supply  limited  cohorts.  CANTOT 
is  total  canvasser  limit. 

RECCST,  CANCST  relative  recruiter  and  canvasser  unit  cost. 

Set  to  1.,  1.  for  best  convergence  where  1 may  actuahly 
represent  1 unit  of  $25,000. 

FSERV(k),  fraction  of  total  supply  with  form  of  service  type  k. 
Fraction  of  total  supply  in  combat  arms. 

Fraction  of  total  supply  in  region  k. 

Set  indexes  for  loops  later  on. 

Fixed  cost  is  set  to  zero  since  it  does  not  influence  the 
optimization. 

Initial  counter  (for  what  I am  not  sure). 

Initialize  counters  for  no.  of  DRCs  in  each  region. 
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A’;'’"  ,n'i  COPY 


OPT  a 1 


f t:j 


• //  /*Il  JN  a , I?) 

• « ; r a , jt  r«ts, < tk phs 
» ta . jr •» Pt*:,v te *«s 

..  a , ij,  3*,  ?7*no.  supply  limited  types  * , 

«.  wP  3 IF  F . jF  SERVICE  = .13) 

•" c . - 1 

is  r c *« • : c nfl c t 

, t : - * ,:2,3x,  rhic^ct  a ,i2> 

i * : :**,  jc*pcr 

t * . * * tU*P  ,f  1C  TO» 

. ^0*M*NT  FACTO9  R Su  . uASc  CASE  = ,F  3 . 5 , SX  , 2 6H-  A 

- - Ti  a #F5.3) 

. . < . '9v~~:t  r 

' , :&t.  * :r 

'5  ■ nn.  J P ° l T a ,F10.0  , 5x,  1 JhTOTAL  CANVASSERS  = ,F  10 

. j zj-li  . c a *i r : r 
*.c:sr  .cans:* 

• ,**r .C’JIT  COST  FACTOR  = ,F7.1  , 1C  X ,24hC AN VASSER  CCS 

.rM) 

* . , - . . . , :a\cst 

t?,-)  ir  JLV<  < \>  t<*i#<TERMS) 

r:  (3,^-i 

' ;•,  .-Fraction  Or  SUPPLY  LIMITED  WITH  TERM  OF  SERVICE  TYP 


a « 

» : ; : i 

, FCER  /(  <) 

*3 

5/ 

WRITE 

M-.X,:C,2X,F5,4) 
f 5 , - > FCO^TT 

( - ,49) FCO  MOT 

49 

” • “ 

T (//E  < ACTION  OF  SUPPLY 

<c,-)  (Fr£0.\  > > , < = 1 , NC.REG  M 

9 2 

**1  r: 

( , 5 2 ) 

52 

C-  i- 

’(//n,«.  “FRACTION  OF  SUPPLY 

^ - 1 | *.u'  » gN 

iu 

V?;:I£ 

(c,  5*)  < , F R £ GN ( < 1 
a MAXSTA  ♦ 1 

95 

>*7  7 CA 

a 2. • < T £ R MS 

FIXSCT  a 0. 

Ncounr  = o 
u3  9 T = 1 # f.OPEGM 

9 I-.Eo‘.(  I ) = c 

OC  50.'  L La  1*11 

5 , :cj)  k word ,12,  u,  <val  d ) , 1=1 , jtypes) 

IF(iOF (51 ) 1C-- ,15 
10  FORMAT <A5,I2,2x,Il,3Fd«C) 

12  FORMATC  9X,Ap,:2,2X,It,3(2X,Fll.5)> 

10  WRITE (&,12)  KWCPO ,12, II , (VAL (I ) , 1=1 , JTYPES) 

IF (KWORO.EC  .OhPCTVL)  GO  TO  403 
30  2C  1=1,11 

I VAC  a 3 

IF  C <WOP  O.EC.INAME  CD  ) IVAL  = I 
IF<kW0P3.£C.lNAM£UMG0  TO  25 

20  CONTINUE 

IF( IVAL.cQ. C)  wRITS(6,21) 
iF ( IV -L ,£C,Q)CALL  EXIT 

21  FORMAT  C//1CX,  1C  <*-*)  , ’UNRECOGNIZED  NAME  CARO  - EXIT*) 


ANNOTATION 


SUBROUTINE 

STATEMENT 

NUMBER  SUBROUTINE:  READIN 


100-175 


150 

153,154 


157 

158-160 

162-169 


Basic  input  loop.  Everything  entering  here  should  be 
preceded  by  a name  card  listed  in  DATA  INAME.  A card 
with  the  word  END  cols  1 to  3 should  end  each  data 
section  in  this  loop.  There  is  no  total  counter  check 
within  each  section  of  data  so  that  for  example  if  the 
user  states  he  has  63  DRCs  then  he  must  have  the  appro- 
priate information  for  all  DRCs. 

Formal:  10  is  the  primary  form  in  which  the  data  is  read. 

Currently  some  values  of  the  environmental  factors  are 
the  same  for  all  cohort  types.  This  may  change  later 
hence  right  side  VAL(l)  could  be  changed  to  VAL(J)  123, 
126,. ..,139. 

The  key  word  PCTVL  sets  up  the  branch  to  400  which  reads 
in  the  historical  profile  information. 

FORMAT  401  is  fundamental  to  the  historical  data.  The 
three  indexes  IR,  ICA,  and  ISER  identify  the  region, 
combat,  noncombat  arms  and  type  of  term  of  service, 

FVAL  is  the  total  supply  enlistments  by  DRC  (it  is  not  used 
in  the  calculations)  VAL  (J)'s  represent  the  discrete 
breakout  by  DRC,  by  term  of  service,  and  combat  and  non- 
combat arms. 

Exit  from  historical  data  input  section. 

First  card  for  each  DRC  establishes  name  and  tallies  on 
region  counter. 

Combat  arms  and  noncombat  arms  are  separated  within 
term  of  service. 
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pr.TT  • • »r>l  C PODV 
Dlji  «.  ..  -I  Luri 


*3/72  CPT  = l 


F T ■'!  4.;*R43C 


.2/1  ' /7 7 ■ . -.  j.  S2 


Ir(-.r.'-'0.t2.5HC!J[)  I GO  TO  53  0 

11  - :val  - i 

r-C0ui::  = mount  * t 

r • ‘Otv  , ID  < C,  12,  Ii  » (VAl  U>  f 1=  1.  JTYP£S) 

' i f i I 1 i ! i I Km. 

I F ( .< „ ; - ~ . ii ; . b n £ N 3 ) GO  TO  5 0 3 

o3  T'  ( 11C  *120  *133*140  , ly  0 * 16  J *17C,15C»19G)  » 1 1 
JO  lit  J-l,JTY->cS 
KLi.Ji  = VAl<1) 

GO  TO  27 

121  J=ltJTYP£S  . 
w ( : 2 , J ) = V-i  L ( l > 

GO  TO  27 

cc  iii  jsi.jrv^es 

U ( I 2 » J ) = VAL(l) 

GO  TO  2 7 

"■)  iwi  J = 1 , JTYP£S 
AT ( I 2 » J)  = VA L(l) 

S3  TO  27 

1 2 151  J=1 ♦ JTYP'S 
AOCIZ.J)  = VAl (II 

~ZZ(Z?)  - V AL  < 1 > 

GO  TO  27 
CA\fI2)  = VALC1) 

GO  TO  27 
iwASTLCIUES 

23  m J=1,JTY?£5 
E ( i 2 . J ) = VAL(J) 

30  TO  27 

CONSTANTS  FROM  X sectional  ANALYSIS  ••• 

CO  10!  J=1,JTYP£S 
CON ( J ) = VALCJ) 

GO  T r 27 

P-:  vGT  \ T AG£:>  TC  3£  C32£0 
CONTINUE 
20  *15  I = i , IPIG 
20  412  II=t,<TTCA 

R£A0(5.4CilKW0F0.ISTA,  IR , I C A , l SER  , (VAL  ( U > fJ-l  «<))  # F V A L 
FORMAT  U5.  I2.2X.3I1,  8F5.4,5X,F10.1> 

K.xli:  (c;,«,22>K«O^C,  ISTA  , IR,  ICA,  IS£S,  (V  AL(J)  ,J=1,3)  ,FVAL 
T ( i:xf  A 5,  2X  ,I2,2x  ,3  ( 1 1 ,2X  ) ,3  (2X,  F f.D  ,5X  , FlJ  .1) 

I F ( < * 0 *.  0*wC*5h£nJ  ) GO  TO  5 3 0 
:c(II.:C,ll  1R£G\ ( I < ) = IRLGNCIR)  ♦ l 
IF  ( H.SQ.l  ) TCT£,-iL<I)  = FACTOR  * F V AL 
IF ( I I . £0 . 1 ) NMCRC(I)  = <W0R0 
I F ( ICA  .£<3.  1)  GO  TO  438 

oo  *-;5  jn.a 

CMASUi  STATION  I N?£X  L ATT £R  »•* 

PCTNCA (I.J, IS£R>  = VAL(J) 

GO  TO  4 1C 

GO  4C<A  J = 1 fd 

PCTCA( I,J, iSERI  = VAL ( J) 

CONTI  \N/E 
CONTINUE 
Vi.RI  T£  ( 41  6) 

FORMAT (lOX.lCC*-*) ••  NO  END  CA F3  ENCOUNTERED  AFTER  PERCENT  COHORT 


SUBROUTINE 

STATEMENT 

NUMBER 


173 

177-186 

187-203 

204-212 

213-216 

217-221 

222-231 


I 


ANNOTATION 


SUBROUTINE:  READIN 


Error  exit  if  illegitimate  card  name  appears  on  list. 

Separate  fraction  by  term  of  service  for  each  category 
type  J and  each  DRC  I. 

Define  percents  of  cohort  type  J for  each  DRC  I. 

Normalize  fractions  for  cohort  type  J at  DRC  I so  that 
computerwise  they  sum  to  1. 

Write  out  total  fraction  (should  be  near  1)  before 
correction  as  well  as  corrected  fractions. 

Renormalize  total  fractions. 

Calculate  C(I,J)  coefficients  for  supply  functions  from 
environmental  factor  values  and  ISWPCT  switch.  Also 
calculate  enlistments  based  on  input  values  of  recruiters 
and  canvassers.  The  values  of  the  x sectional  parameters 
play  an  important  role  here. 
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OUT  /r-U?5  r rftr\w 

Dwi  >u  5;  LuPY 


»’3/7  3 OPT 3 i 


F T N 4 . 5 ♦ * 4 C 6 


i ?a  r a - exitm 
c:.l  £*:t 

GO  TO  vuC 

5cc  continue 
icj:  con t i-jl'-: 

00  131c.  I=l,MAxSTA 

oo  i: 1 6 j=  i , -j 

S*JM  = 3. 

co  ::::  < = i ,k  terms 

P0TT«C,J,O  = =CTr.CA  ( IfJtO  ♦ PCTCA(I,J,<) 

ir io  so*-  = yy*  * pct  t^mi,  j,  o 

:f  i ^y.zi.  :.o  > oo  to  ioio 
uo  «. :ir  ^-uk  terms 

i:i5  f c T t t : . j * < i = pcttm(i,j#o  / sum 
t:io  co:: r : .US 

JO  1C23  I - X tMAXSTA 

oc  ::c:  j=i,a 
cjm:  = i. 

jU-2  s 0 . 

00  1 J 1 J N=t,KT'£PMS 

5'JMi  = SJ*U  ♦ ° C T C A ( I , J , < ) 

1C  id  SUV  2 T SUV*  ♦ r*CTNC A ( I » J*  <) 

ca  = ct  ( : , j)  = suv i 
tc,j)  _ S’jv2 

T £ • * = SUM l ♦ SUV2 

: -it-:  w)i:i/,uiMn9 
i:i  7 p’ctca  c*j,  i)  - :.c 
PP3TCA  ( I * J » 2)  =0,3 
GO  to  :3  20 

i:  19  p p c r c *.  t : »j,i>  = suvi  / term 

P»C  TCA  ( ; ,J  , J)  = SUU2  / TERM 
13  23  P^'C.w)  = CAPCTCI.J)  NCAPCTUfJ) 

oo  i : j : i=i,maxsta 

S’JM  = o. 

CO  1025  J=  1 ,d 
1C  2 E SUU  = SUM  f PCT  < I * J > 

1 w 03  PCrSTKI)  = S'JM 

CO  11C2  1 = 1 t MA  XST  A 

ro  i:;i  j=i*a 

1331  pctci.jj  = pcrn.j)  / pctstacx) 

133?  CONTINUE 

00  lu4y  I r J.MAXSTA 

KRIT.*:ib»U41J  I .FCTET  A (I)  * (PCT  (IiJ),J-l,5) 

1041  FORMAT! 1CX. 12, 5X, 9<2XfF6,4> ) 

13C0  CONTINvE 

JO  I3a0  1=1 .maxsta 
fctsta  (II  = O.c 

30  1-40  J=1,JTYP£S 

13  49  PCTSTA(I)  = PCTSTACII  ♦ PCT(I,J) 

1350  CONTINUE 

CO  113-3  J=  1 » J T YP£S 

oo  i : : j i=i. maxsta 

c •••  special  input  changes  until  further  notice  ••• 

CC.J)  = ••£(2,  J)  MUC,J>  •F'JKcMP)  ••E(3,J)  1 

1AT(I,J)**£{4.J)*A0(I.JJ**E(5.J) 

IF ( I3WPCT.E0.1  ) C(I.J)  = C(ItJ)*PCT  (I,J)  / PCTSTAU) 

ccitJ)  = c<i,j)  • factor  • ExPtcoNun 


ANNOTATION 
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NUMBER  SUBROUTINE:  READIN 


232-235 

236-244 

245-262 

263-265 

266-270 

271,272 

273,274 

275 

276,277 

278-280 

281 

282 


Write  the  values  of  the  C (I,J)  coefficients 

Write  the  enlistments  total  and  by  cohort  type 
for  each  DRC. 

Set  indexes  for  constraints  and  indexes  for  begin- 
ning and  ending  DRC's  within  each  region. 

Set  RECO  (I)  to  the  initial  input  value.  Set  X (I). 

Set  initial  value  of  canvassers  CANO(s)  to  initial 
input  value.  Also  set  X (I)  for  canvassers. 

Set  signs  (the  instructions  may  currently  be  obsolete). 

Initialize  right  side  constraint  values 

Set  right  side  (1) — the  national  goal  to  supply. 

If  canvassers  are  in  the  model  right  side  (2) 
upper  limit  is  CANTOT. 

Set  right  side  constraint  values  for  term  of  service. 

Set  right  side  constraint  value  for  combat  arms. 

IREGSW  is  set  to  one  by  default.  At  one  time  it  was 
used  to  bypass  regional  constraints. 


283-285 

Set  right  side  constraint  values  for  regional  require' 
ments. 

286-292 

Write  the  values  of  the  right  side  constraints, 
canvasser  fractions. 

Set 

293 

Obsolete  possibly. 

294-296 

Set  canvasser  regional  fractions  proportional  to 
regional  goals. 

297-300 

Initialize  array  of  summary  solutions  as  used  by 
routine  SOLN(MT)  later. 

sub- 
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readim  is/?:  r.r* t = i fin  4.s**406  02/10/77 


:•:!■(!,  J>  = C(I.J)  • xE7(I)**E(u,J> 

If  ( liwTC-i.rjE.  1 ) EM.d.j)  = -t._  (I , J)  • C AMI)  **E(7,J) 

1 i ? a c"  ri> uf 

T£  (o.  1 jit  7) 

1547  fz-  f*i  r ( // 1 ox,  i c ,j  i ,3  nr  typei 
.3  i.'.'  I - i ,rs  >sn 

:2cc  w-.ite  ( . ,12 ig>  i.  (- a,j>  ,j=i,j types) 
mo  r r ( i o x , I 5 ,5 « i < , f t o . o > > 

<••=•:  n ( e , 1 2 n > 

1211  F 7-MT  (//  jx,20  I*-*)  ,*£NLIST:i£NTC*,  20(*-*)) 
m)  i‘2o  :>i,MAxsra 
= c. 

<77  1315  J = l ,jrro£s 
1215  S2“  = ;j-  * E.nl  C,  J) 

:rr ! 7, i2io)  ienlc.j),  j=i,jty=es> 

1220  C 2 * • 1 1 1 ( U £ 

••  : : .:r»;-,r  ingexeg 
11=1 

•••  CANVAST;'  CONSTRAINT  - IF  AMY 
12  - I 1 » l 

IF<  IswTGn.iC.  1 ) 12  = 11 

<rx^'-s 

:»  = n * i 

ii  S [1  * JTYFCS  • NOREGN 
I- CnEGSV..  EG.  l ) 15  = 14  ♦ NCSEGN 
lr(i-'.03-,EC.3)  15  - 14 

N - “AxSIA 

If  (I  rOf.Ne  .1 » M = 2 • M 

Ms  15 

PIOINIll  s 1 

I;  ‘.0(1)  s IP£  G '.'!!) 

oi  i£on  i=2.mc7£gn 

1:70(1)  s 0(1-1)  ♦ ISEGN(I) 

15  0 5 I 7 2 G I M I ) S I£  (.0(1-1)  * t 
00  1710  I=i,IA<STA 
r :C2  ( . ) = ICC  ( II 
1710  * ( I ) = 7.EC  ! I ) 

:f(io4t:-i.m-.2)  go  to  i7ig 
GO  1712  I=1,YAXSTA 
CAMO  (I)  s CAM  I) 

II  s -ZtSTi  ♦ I 
1712  x(II>  s CAM(I) 

171=  20  172C  J=l,JTYPES 
1720  StGN(J)  s SIGNS(J) 

00  1705  Isl(50 

toco  pcmu)  s o. 

iv3T2C(l)  = GUP  LY 
IMtSwTCH.Efi.ilOO  TO  1305 
R3I0E (2)  = CSNTOT 
13  05  LO  1U0  K=1,KT£RMS 
* s u , K 

1310  PSIO-Mi.)  = FG£FV(<)  • RSIOE(l) 

4GI  C£  ( 14 ) s FCC*GT  • o.StO£(l) 

IF  { IFEGSn.NE,  1 ) GO  TO  1830 

00  1320  X=  l.UCfCiG.N 

1 = r-  » K 

1320  AGIOS ( 1)  = FREON(K)  * RSIDE(l) 

1333  C-.TIM':: 

•.  I T £ 1 6 . I 7 3 9) 

1735  ‘ ; .‘A  T (5  ■ , • f IGmT  SICE  CONSTRAINT  VALUES*) 

03  174  1 1=  1 ,15 
JANF-.ACII  = 1.0 
1742  w = :TE(t,  1741)1,  P.SIOt(I) 

1741  f r ( i jx , 15, 5x,f  io . i ) 

I SC  T = 0 

1.1  = 14  ♦ 1 

" 19:2  1=1-1.15 

1333  CAMFwAC)  = RSIGE(I)  / RSIDE(t) 

05  144  2 I=1,*1AXSTA- 
GO  1940  .VT  = l,  15 
> .7  •■;.(:, -T)  = C .0 
19 4 2 X = £C(  I «“T)  = 0.0 

4 4 F ) -.MAT  ( 12,1X,A5,2X,5F10.5) 

.-.£  rup1; 

L'.J 
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ANNOTATION 


SUBROUTINE 

STATEMENT 

NUMBER SUBROUTINE:  READOU 


This  subroutine  is  the  Report  Generator 
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BEST  AVA1U3LE  CCPY 


SUBROUTINE  RZADO'J  7 3 / 7 3 CPT  = 1 


FTw  h.5*R4:6 


22/13/77  2i.-5.J2  =>ACc 


ir(<.£2.2»  W’lTil'i.HD 
ip<  ir  ;:,T  = i 

.3  I r ( < . Z fi  • 2 ) IPRINT  = 2 

95  3 fCrMA  T(  Or  *£*/I0E  - 7-"-  *E  YEARS*) 

*31  m six, •term  cf  s e.-.vice  - c"',.w  or  .ic.RE  years*) 

c r : -*.:l  speciality  constraints  ••* 

IF  ( :rw:r:  1 ) WRITE ( 6, 9:1)  (J,J=i,JT<=_5) 

o5  IF  C lu*  TCH.NI,  l ) w->aT£  (6*995)  <J  tJ-1,  JTYPES) 

DO  2 1 ^ ll= : .ncreon 
PSu/R  - C.O 
PS'J'u  = y.3  , 

PS'JMA  = L.O 

73  ILL  = LL 

IF (Ll.GT.il  LLl  s LI  M 
W ® I T E ( c » ■*  5 4 ) LLL 
= I OEG  * N (L  L ) 

IS  = IIVC(LL) 

73  CC  2: ? J=1 , JTYFES 

232  FCj^-3  (J)  - 0.3 

D*'  213  : = !•?,:£ 

S J ' A = i . ' 


C,,v3  = * rv  £ C ( I ) 

: S L v C - Sv'*  . ♦ CA  ,*:(!) 

= hSU"®  ♦ ®ZCCU 
= r --c  = ^ y.j^c  ♦ cam  i) 
c : 2 : - j=i,jty?es 

T £ T r P ( J ) = E:l(I.J)  * PCTT  M ( I * J v K) 
EjVTYP(J)  = SUM  T/P ( J ) ♦ T Z\ TYP  ( J ) 

S’J '!  A = SL.-.A  ♦ T £ “ T Y P ( J) 


P.S'VTPJJ)  = PSL.rTPCJ)  ♦ T£RTYP(J) 

WJ^A  = r.S:-PA  f T £ R T Y p ( J > 
ruv  2 sum  ••  t£tTyp(  J) 

}:  22  4 c:  4t:*;u? 

Ir  CSWTC-.EC.l)  WR I T £ ( 6 , 992 ) NKQ3C  ( I ) t I.  R£C(I  ) 
l(Jl  , J=  1.  J“YPCSI 

IF ( :_wtCH.Nc.1>  W.PxTE(6,592)  NHCRC ( l > * I » *£C <1 ) , CAN ( I ) 
1 ( J ) , J 2 1 , J T Y P£  S ) 

35  212  C 0*1  T I Ni'  Z 


. SUM  A • (TE.RTYP 
, SUMA *(T£RTYP 


133 


1C  5 


113 


1F(  IjWTCh.EQ,  1)  WRITE  (6*1075)  RSU  MR , PSL'f-A,  ( RSJ  M T?  ( J ) , J=  1 , J T Y *£  S ) 

IF ( Ir WTCH.NE.il  WRITE (6, 13  75) RSU MP  ,3SUMC,RSUM A, (RSUMTP(J)  , J=1 , JTYP 
1ESI 

P®S'JMR  ( LL  I •=  RSUMR 
°f.:j"A  (LD  = RSUVA 
P-\EUMC<LL>  = RSUMC 
30  21*.  J = 1 * JT  YP£S 

214  PSUM  T3  ( s.L  * J ) * RS'JMTP(J) 

215  CONTINUE 
GO  TO  1333 

330  L = v T - 13 
W R I T E ( 6 . 111 

IF<  IST^AT.EQ.C  3 WRITZ16,9S9)  ITITlE 
iF(IST-  iT.EO.t ) WRIT£(6*993I  ?SI CE CMC > , IT  ITLE 
IF(IST“AT.c3.2)  WRIT  • ( o,  E^C  ) .RSI  CC  ( "C  » ,1  7 ITlS 
IF(  1STPAT  . £0.3)  W.RI  T i ( o . 2 ? ? ) RSI  OE  ( VC)  * ITITLE 
IrfirTr.AT.CQ.4l  wp;TE<6.3->C)  R SIDE ( MC ) * I T I TL£ 

Ir ilSTPAT.fcQ.5l  WRITE (o,9;d)  RSIOt (1  I , I TITLE 
IF(L.ZC.l)  wrITE<6«902) 
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k 


1/4  M an? 


copy 


L 53 


133 


16: 


173 


73/75  CPTsl  FTN  4.5*R4Q6  02/1  Z/ 77  C 9 . ^ 5 . -S  2 


: " (l . : :.  o *'3.  r£(ov  d 

ie(l.  .i)  • 3 

5*:  fj  -;(  -c«  t:m:al  rjtCi:i:Tr  - ccmja;  arms*) 

Mo  r ' - * • r i i :>.* 'CC-Mat  : : ^=.c  ia.ity  - .;..-cum3A7  arms*) 

(:•>?  > -.v:rc<s>,M9i)  u,jsi,jty?es) 

I r.E.  1)  "RITE <6,5 95)  (J,J*1,  JTYPtS) 

J ~ 31  v L L = 1 *f.'  GN 

‘•yjf*®  - 

•:  :*j  M : = : . i 


3.2  J=t.JTY?ES 
3-.  RG'JHT*  I J)  s c.O 
L « L = lI 

I'tLL.GT.l)  LLw  s ll  ♦ 1 

»-ITC(c,H4)  ill 

m 

::  = :_‘0(ll) 
co  3i“  : = r-i , : t 
SJ'*A  = C.C 

■j.M  - S V" - »•  PiCCI) 

*2  - Su-c  ♦ r. .*.  *•  ( : > 

- « 36 “ ♦ -.£ C ( : ) 

kZJ  '3  = ?wJ  « CAM  I ) 

-0  3 * • J=l, JTYPZS 

t : • t r - c j ) s ; lc*j)  • p»crc.i < :» j,l> 
s iTYriJ)  T.RTYP(J) 

:S>vT-fj»  s fs->-tp<j)  t-rtyr<j) 

•IJ'A  = S M A * 7E?TY*<JI 
r'SJ^A  = RSJMA  ♦ r£;,TYP(J) 

:,1  = 3 .. M f T £ R T Y P < J > 

3C  <•  CMTIMjE 


Ir  { 

cZ  .1) 

WRITE ( 

6 , 

592) 

NMORC<I>  i 

,suma  .trs-UYP 

1 i J) 

. j = : , jrr:£s> 

: ( 

I £ <■.  T£ m • •:  £.  i ) 

WRITE ( 

0 , 

v 93) 

N MO  RC  ( I ) . 

i i.R£C<;).CAN(i).suMA.<rEsjYP 

1 1 j) 

* J- 1, J TYPES) 

3 1 C COM 

T l n l £ 

I r ( 

IGWTCK.EQ.  1 ) 

WRITE ( 

o » 

10  7 5) 

=,SL'UR,?S  JMA,  <*$J  MTP(  J)  , J=l,  JTY^iS) 

lF( 

ISMSH.NE.  l ) 

WRITE < 

6 » 

10/5) KSUMS,RSJ fC, rS JH  A « <.RSUMT?(  J>  , J=  1, JTYP 

1 E.S) 

P^UM3  (Ll>  = KE’JMR 
(L->  = RiUMC 
FK^UMA(lL)  = P.SUMA 
CO  314  J= 1 , JT  YPES 

314  P£U«TP(Lw,J)  = RSUMTP(J) 

315  CONTINUE 
GC  TO  1 G JO 

••*  rESIPSAL  QUALITY  CONSTRAINTS  ••• 


MTT  = 

wy  _ T {* 

IF  (MS 

• f,E.  y T T)  G 

O TO  23JC 

If  ( IF 

£GSw. £0.1) 

GO  TO  430 

mm  - 

(NTT  • 1 > 

/ JT Y P£S 

1 

M“M  = 

MM 

:c(  -m 

. GT.l)  MUM 

= MM  f 1 

A 3 S 

IDEGIN(MM) 

IE  = 

IEMQ(NN) 

J TV  P£ 

a MTT  - ( 

MH  - 1 » • 

JTYPES 

GO  TC  9CQ 


PAGE 
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BBT  AVAIL'D! E COPY 


fc/74  CPT  s 1 


F T N <*«  5*itC6 


C2/1C/  77  G9...5.32 


430  p:irr'!N(''ITI 

> « i s - r r 

: r < m r r . r . i > f'l  = h r t ♦ i 

g:  to 

• ••  LZC.ZZr  : ON  CONSTRAINT  •** 

90 * r. r.T :*:t  ? 

w : t c ( -3 , 1 1 > 


:c<  ::t-at.ec 

. 0 

) WRIT 

C(i,Ti  )l 

irrrLc 

: © c : : r **  *.  r . : i . 

t > 

vi  © . r i i 

9, 1 > 

3 

zizzr'2) 

, IT ITc£ 

:r« : ~r. it. :c. 

2 ) 

w s* : • : ( 

, C 93  l 

3 

. IT  I T L£ 

I- ( : ji , *.  -g. 

! 1 

* - :th 

, w ; 2 5 

K 

i: :e(mo) 

tITITLi 

: f ( I : r ••  : T . C 0 . 

- ) 

^:T£l 

? ,3  ?2  ) 

j 

JI7  2 <MC» 

» I T I Tc£ 

u:r-AT..o. 

5 > 

WRIT  £ ( 

6, 9i1> 

SICS  < 1 ) 

,iriTc: 

9®,  4 f ' ' M A T {/>*)>  «• 

A £ 

9 1 C : •. 

17) 

99?  =c  <-•*.  r 5 : * ,*4 

: ( 

• - • ) / 9 X 

, *:r-{A 

J i 

GY  - Ml  4 

coot  : 

G 1 J £ C T 

TO 

NAT  10NAL  TO 

:ta_  accession 

: : j : 3 s •» 

* *•  r • 

f : 

o/  i < » i 

JA 12/ IX 

* 133  <•- 

•)  ) 

99a  F0<vatcx  ,12 

: ( 

♦•‘i/jx 

, • 3 T © A 

r : 

jY  - v£ N 

LujT  j 

uGj£GT 

TO 

REGIONAL  Tj 

» r :r. 

£ -.  v r * - " * 

*;f  ♦ 

r t 

.,3/LM 

3A1C/K 

, 12:  (•- 

•)  ) 

39  r F 0 > s‘ A T ( *.  X , lw 

: < 

, *CT9A 

T£ 

GY  - v'\ 

cost  : 

G 3 J G T 

TO 

CCMOAT  ARMS 

i c::ji-£M-f4T 

• 

- 1 : . : / 1 

<*1311 

:/ 

:v , i::  <• 

-•) ) 

29:  f;;k:tc. a 

2 

*-*) 

i * - T t .*. 

T £ 

GY  - Vt\ 

GOST  3 

U3  J- CT 

TO 

3 YR.  TERM 

; 

•,F12 

, ' 

/ 1 < , 1 7 •.  1 

3 / 1 X , 14 

29  2 - "match  ,1a 

2 ( 

• - * l / 5 * 

, • S T H A 

7£ 

GY  - v£  { 

GOST  SUTJECT 

TO 

4 YP..*  TERM 

i :r  o - 

£ 0 

u: 

r •.Ft 

V • 

:/ :* , iii 

20 ( *-•)  ) 

: -*  ( : „ - r 

1 ) 

.>  >; : r £ ( 

a , ? 'j  i ) 

( J 

, j = i , j r y 

- £3 ) 

IF(  I ~ - . : . £ . 

1) 

« - : T r ( 

3,995) 

( J 

,J=1  ,JT  ^£S) 

991  c0..w:r{/*<,V 

• • 2 X , * N 

0.*  ,2X 

NO.  F£:9JIT£RS* 

, 3 X , * N 0 

. ACCESSIONS', 

i ox , * i£w£  - sc  c? 

t 

nr  : y 3£ 

• / 4 5 < * 

«,  ( *T  Y PE  *, 

11 ,3X)> 

DO  919  LL=l»f 

3 3 

£GN 

RHJMR  = C.C 
P.S'J.MC  - 0*0 
RSUVC  = C.C 
co  ?:r  j - i , jtypes 
902  \SJMTP(J)  ■=  O.C 
LLL  = LL 

IF (cc . Gf. 1)  L L L - wL  ♦ 1 
♦•HIT £ ( o , 944  > LLL 
13  = P'OINIU) 

I£  = x£.NO  (lL) 

00  911  I-ID  ,IE 
SUk*R  = SUMR  ♦ SEC  (I) 

su-:c  2 sums  ♦ can (i) 

RS’Ju-f  = ♦ R£C  C 1 1 

RSUk*C  = RS'JMC  ♦ CANd) 

SUVA  = 0,0 
00  934  J=i*JTYP£3 
T£*  T f 3 < J)  s Em(I.J) 

SJwTv«MJ>  - SUV  T Y 3 ( J ) ♦ TERTYP(J) 

©CUMT-MJ)  s PSUMT°(J>  ♦ TcHTYP(J) 

SUM  A - GUM  A ♦ TSRTYPCJ) 
kSUMA  = RSUHi  ♦ T£©TYP(J) 

SJM  = SUM  ♦ TERTYPCJ) 

9C4  CON T T NU£ 

IF(  ISWTCH.SO.I)  W.9IT£(6*992)  NMCKC  (II  , I • RSC  C ) * S'JM A . ( TCRT  f P 

1 < J) • Js t , JTY©£S> 

IF(  ISwTCH.Ni.  1 > WRITE  (0,992  > N VO.'; Cl  I)  , I,R£C(I>  , CAN  ( I > , SUM  A , (T£  RT  YP 
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7 2/7.'. 


cpt  = i 


r T'J  4.  v*HC6 


C £/  1C/ 77 


1 4 ;>  , ; • i , jt  i p.ts  ) 

>»:;  rr.vji 

:f  ( ; ...  'Cn.  e s. : ) w-i  tz  (6, : «7  5>*  cjm,  rc  j.^a,  iksjmtm  j)  • j=i,  jt  y^es) 

:c  ( tc*i.  **.:Tit6,;o75)  icump ,->si*c ,,^u,  usu.m tp  <jj  , j=i,  jtyp 

: £ ; i 

F-  ;u  MLl  I = ?sumr 
'**  » J • - ( LL ) = RSU :A 
r . -u  ill)  = ^SLMC 

: •'  1 **  J-\  * J”  Y ^F. > 

9:-  f ...t“(Ll, J)  = tp ( j) 

919  CvM  . -ui 

go  r?  tea) 

995  F£-.-  • r t /5x  , ?'<,  ,2X,  *f.O.  c.£C3UZri«S»  ,3X,  *kO.  CANVASSERS* , 

uj>.*  vo.  access: : *Z9*  i *, ‘accessions  .9*  r*?c*/o.;x.*(*TYr»E  *,n,3xn 


952  L'  * * :?  X , A E 

959  F : ' r ( 14,1 

to  . ;\plt  • 
i:o : cc-:t:-.je 

: - 1 r'* i tch. £: 
: - : ; - r . •:  L 


1 2 , * : : t . * i l . i > ) 

: (* -M /5X  ,»:j3  OPT  If*  1’AT  ION  Sd-UCSY  - cyALUA  T 20 .4  3ASc 
L’JlS  :F  REJR'JlTiNo  F0i-.3i*/ixfij;  13/ IXfi  3-3  (♦-•)> 


:r  l I'-xTCM.E:.  1)  /MITE  (6,  135:>  SCV.R.SUM,  (SJ'-TfM  J)  , J=l  ,JTY=>ES> 

1 :-  t ;>.*-.:)  - :te ( bf  1;; :»  , 3wm: ,s-jm,  ( su.mtypi j)  , j=i > jtypesj 

i.52  F . 'M  * l?Xf£- /-iX,  • TOTA.S*  * MFi:.  1,  2X)I 
i:75  F*  v * (£Y,  ;"(•-*> /EX,*  SUBTOTAL  S*,  MFl3.t,2X)) 

c ’ ;c  l - 1 , • 

p-  . < l ) - • p3.S'jm(l)  / s *jy* 

P . _ r:  ( . ) = 1 : : . • / £ (..*•* 

: r ( : *c • »c  • 3 • : ) 9=:u*c(l>  = ica.  * prslmccu  / sjmc 
sj  > 1-15  j- : ,j r r»:s 

( j)  = 


1,95  PJ.- 
150  0 C~.J 


p ( t * j ) = t c : • : * 3$jmp(l.j>  / sum 

.ut 


v.  - . r : < 5 , 1 9 : ; > 

no:  F„>  M /..?  <,  10  t •-*)  , *p.>.Cl*:T  j'JI MARY  TAbLw*,  l3(*-*>  ) 

ire  isr-AT.  e:.  : > wRiTEd,  939)  itttlz 
(rrTvAT.cO.i)  w>;ri(cv93  j)  rsicf  cmo.itike 

IF  (ISTRAT. EC.2)  rM.I  T£  (6,29) ) <ST.0E  (MO  , IT  ITwi 

if< istrAT, £ j)  Rsn£<MCi  ,ititl£ 

IF(ISYQAT.iS,U»  WP.ITK5, 393)  *S  10  £ ( M C)  , IT  IT  L£ 

IF(  ISTRAT.  £C.v)  WRITE  (5,-993)  3SICEU)  dTITLE 
IF (*■". £:.Il ) wPI T£<  6, 1675) 

IF < If  dNT.2  ••  1 ) ^ - .Tin, <950) 

IFf IP?INT.cC.2)  WN;T£Cto,93i) 

IF  C'tlfiT  • c2«.3)  n-;IT£(6,9«)2) 

1675  FC?**.M(  50  X,*  TOTAL  ACCESSIONS*) 

S4  - 3. 

SC  = c. 

= 3.0 

Ir  ( r.^TCH.tM)  *9 1 T£  (6,1650  I ( J,J*  1 , J TYPES) 

IF  ( I3.4TCh.Ni.  1)  T£  ( o,  15?:)  ( J,  J~1  ,JTYP£S) 

1660  FC3r. *.T  (/17Y  ,*PtC>\U:T£9i*,  IX,  *ACCcSS  I 0NS*,UX,7  (•TYPE*  ,t2,6X)) 

1 67C  FORMAT  ( / 17  X *•  RfcC^vJ  i Tc  »<S*  , 1 X , • C ANVA  SStRS* , 1 X,  * ACC£SSI  ONS  • , 4X,  6C  *TYP 
l£  • , t ? , 6X  ) ) 

OO  1 7 J C L=i #NJ9£GN 
LL  = 1 

iF(L.Cr.l)  LL=L  *■  1 
W=*:t£(6,93-)LL 

:p<  n^rcH.cio.  i )wRir£<o  ,to75)PRSuM^a> , prsu ma(  l ) • ( psumtpu.  , j)  , j*i, 

1 JTYP£S) 


I 
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ANNOTATION 


SUBROUTINE 

STATEMENT 

NUMBER SUBROUTINE:  REQMTS 

This  subroutine  performs  two  functions.  If  KT=*1  the 
right  side  values  for  the  recruiters  and  canvassers 
is  modified  when  ISWITCH=1. 

If  KT=2  then  the  coefficients  C(I,J)'s  and  fractions 
for  combat  arms,  terms  of  service  are  appropriately 
modified.  Some  right  side  modification  is  also  made. 

20,21  Set  beginning  and  ending  indexes. 

22  Set  sum  total  to  zero  initially. 

23-27  Branch  to  appropriate  class  of  constraints  for  given  index 
MT. 

28-30  Error  exit  for  index  out  of  range. 

32-42  Evaluate  total  cost  for  recruiters  and  canvassers. 

44-46  Readjust  right  side  values  for  canvasser  total  and 

recruiters. 

48-59  Modify  right  side  or  B(I,J)'s  for  term  of  service  depending 
on  whether  K=1  or  2. 
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pacsDim  a-**-"0*  rtuaB.  _ 


S' 

i r.:  ? 

:oh  r 

5 7 5 ✓ / J OPT=i 

FIN  54^436 

1 

su-'irt  njTir.-:  '>iCMr  :r-*r,<n 
Cr  ni^  ::  / * ( : Is  ,fJPl,M*n 

j ir°.  , , r w »*c , 

"0*1 

C ' • .../►v  - O'  ■:/'  «:  ( ; > , J ) ,‘--T  T*(  c5  ,3  , / 

> , :>'>CTCA( 

6=.J,2) 

5 

c *>  : ( •-  > . : ^o i *. t •*>  , : f* 

0 ( 5) 

c;  -"5  ./-  ■ "!  -v u ;n :iv.  t.>. 

n £5SW 

(c-r , *• ) ,<•:“<  c5> 

cc,,,,c../2 ; siomh) 

co  imw*:/  :r-t  /f  i « r , * ; :C'-r,  c a host 

> 0 AN ( o5 ) * 

c.N.  lop,  3) 

ta 

: . :v/r.7-;  ./  «cr;f  <30 

;.i  /•.  1 -j  . ; £ , ic  clm  ,r  l n.-tus: 

> ,5'i3_*<  5 J ) ,.<SI0£0(5C) 
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61-72 

74 

75 

76-88 

86-88 

89,90 

91-102 

104-114 


SUBROUTINE:  REQMTS 


Modify  right  side  B(I,J)'s  for  combat  arms  depending 
on  whether  KT=1  or  2 

Set  basic  value  of  index. 

IREGSW=1  by  default 

This  branch  is  not  intended  for  official  use  but 
is  left  in  in  case  of  future  specialization  to 
specific  cohort  type  optimization. 

B(I,J)  transformation  in  proportion  for  regional 
constraints. 

Set  beginning  and  ending  indexes  for  regional 
constraints. 

Modify  canvasser  and  recruiter  constraint  values  or 
B(I,J)'s  depending  on  KT=1  or  2. 

Modify  total  accessions  remaining  requirements  or 
set  B(I,J)’s  if  KT=2. 
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SUBROUTINE 

STATEMENT 

NUMBER SUBROUTINE:  SEARCH 

This  subroutine  is  mathematically  the  most  efficient 
search  minimization  known.  The  search  is  made  over 
interval  TL,TU„  Three  different  functions  are  called 
depending  on  the  value  of  IFUN.  The  coefficients 
B(I,J)  used  in  the  functions  FVAL1,  FVAL2,  FVAL3 
are  set  by  REQMTS  (MT,KT) 

8,9  Set  the  constants  for  the  "golden  section." 

10-46  Switching  of  end  points  and  function  values  over  a contract- 
ing interval  TL,TU  which  decreases  by  a factor  .618... per 
iteration. 

48,49  Error  message  if  convergence  criteria  are  not  met  within 
60  iterations. 
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SUBROUTINE:  SOLN 


This  subroutine  records  all  solutions  and  writes 
out  solutions  as  they  occur  for  each  constraint 
or  requirement  MT 
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This  subroutine  is  called  only  when  canvassers  are 
in  the  model.  The  W value  calculated  is  a guess 
at  an  initial  value.  It  is  used  to  set  end  points 
for  the  Fibonacci  search. 
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This  subroutine  is  used  to  calculate  an  approxi- 
mation WW  to  the  solution  to  be  found  by  the 
Fibonacci  search.  The  values  IL  and  TU  for 
the  end  points  of  the  search  interval  are  set 
using  WW. 
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8.0  AMRRO  SYSTEMS  IMPLEMENTATION 


GENERAL 

Under  the  AMRRO  Study,  GRC  has  conducted  research  and  development 
activity  toward  the  establishment  of  an  Army  recruiting  Market  Informa- 
tion and  Analysis  System  (MIAS).  The  system  development  portion  of  the 
Study  has  been  conducted  under  Task  V of  the  contract.  The  MIAS  is 
comprised  of  a number  of  subelements  which  are  integrated  into  an  overall 
concept  with  the  principal  focus  oriented  toward  providing  market 
information  and  conducting  recruiting  analyses  at  the  Army  District 
Recruiting  Command  (DRC)  level. 

MIAS  SYSTEMS  DESCRIPTION 

General 

The  research  efforts  associated  with  the  development  of  the  MIAS 
elements  were  conducted  under  Tasks  I -IV  of  the  AMRRO  study  and  are 
described  in  Sections  4,  5,  6 and  7 of  this  report.  The  flow  diagram 
shown  in  Fig.  3-1  depicts  the  basic  elements  that  comprise  MIAS.  MIAS 
can  be  divided  into  three  separate  modules: 

• A Market  Information  Preparation  Module 

• A Market  Information  Data  Base  Module 

• A Market  Analysis  Module 

Each  element  of  the  system  is  dependent  upon  the  other  elements.  The 
market  analysis  activity  cannot  be  conducted  without  established  market 
information  data;  in  like  fashion,  market  information  data  serve  no  purpose 
unless  utilized  in  a meaningful  manner.  The  interrelationship  of  each 
element  in  the  system  with  regard  to  which  element  feeds  another  element 
can  be  seen  in  Fig.  8-1. 

Market  Information  Preparation  Module 

The  primary  function  of  the  Market  Information  Preparation  Module 
is  to  process  raw  market  information  data  into  a form  suitable  for 
analytical  usage.  Since  the  scope  of  market  information  related  to 
Army  non-prior  service  recruiting  is  so  large,  by  necessity  a large 
number  of  data  sources  must  be  accessed  and  a considerable  amount  of 


data  processing  must  be  performed  to  obtain  usable  information.  In 
particular,  there  are  three  primary  processing  flows  which  comprise  the 
Market  Information  Preparation  Module.  They  are  as  follows: 

• Processing  DOD  and  Army  accession  data 

• Development  of  the  Qualified  Military  Available  (QMA)* 
inventory 

• Creation  of  the  Basic  Recruiting  Data  file  containing 
information  on  the  utilization  of  recruiting  resources 
and  socio-economic  conditions  which  affect  recruiting 
success 

Each  of  these  flows  are  discussed  in  the  following  subsections. 

Processing  DOD  and  Army  Accession  Data 

USAREC  with  the  use  of  its  Mechanized  Reporting  System  (MRS)  has 
the  responsibility  of  maintaining  DOD-wide  data  files  on  accession 
information  for  every  non-prior  service  enlistee  into  each  of  the 
Services  (i.e.,  Army,  Navy,  Air  Force  and  Marine  Corps).  These 
data  files  on  accessions  contain  a considerable  amount  o-f  information 
regarding  the  numbers,  locations  and  types  of  individuals  entering 
each  of  the  Services.  In  order  to  utilize  this  information  as  a part 
of  MIAS,  the  basic  USAREC  files  are  first  edited  and  compressed  into 
a smaller  file  called  the  DOD  MVA  (Modern  Volunteer  Army)  File  where 
the  superfluous  or  unneeded  data  elements  are  screened  out  of  the 
individual  accession  records.  The  next  step  is  the  creation  of  a 
"contract”  file  from  the  consolidated  DOD  MVA  (accession)  File. 

A "contract"  is  defined  here  as  either  a straight  accession  where 
the  enlistee  does  not  enter  the  Delayed  Entry  Program  (DEP)  or  a 
DEP  accession  where  the  enlistee  enters  the  DEP  prior  to  actually 
accessing  into  one  of  the  Services.  For  the  DEP  enlistee,  two  data 
records  are  contained  on  the  USAREC  files:  a DEP  record  and  an  accession 

record.  In  order  to  avoid  double  counting,  these  two  records  are  combined 
into  a single  record.  The  result  of  this  processing  step  is  to  provide 
a DOD  MVA  Contract  File  which  contains  a record  on  every  enlistee  and 
the  actual  date  when  the  enlistee  signed  the  "contract"  to  enlist  in 
one  of  the  Services.  The  actual  "contract"  date  is  important  because 

* 

Subsequent  to  the  conduct  and  preparation  of  this  report,  the  QMA 
processing  responsibility  was  transferred  to  MARDAC  for  annual  operation. 
Therefore,  this  element  of  the  MIAS  will  eventually  be  deleted  from 
USAREC 's  responsibility. 
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that  is  when  the  individual  makes  his  decision  to  join  a Service  rather 
than  when  he  actually  accessed  into  that  Service. 

The  next  step  in  the  processing  sequence  is  to  create  a DOD  Contract 
Summary  File  by  Army  DRC  boundaries.  This  process  involves  the  consolida- 
tion of  the  individual  contract  records  into  a summary  file  containing 
contracts  by  Service,  age,  sex,  race,  education  level,  mental  category 
and  DRC  location.  A county  Fips  code  to  DRC  cross-reference  file  is 

: k 

required  to  identify  the  DRC  location.  The  final  step  in  the  accession 
processing  flow  is  to  create  a Historical  Accession  File  by  DRC  to  include 
accession  breakouts  by  cohort  categories,  occupational  specialities  (i.e., 
combat  arms  options/non-combat  arms) , and  term  of  service  enlistment 
options.  This  is  accomplished  for  any  prescribed  period  of  time  from  the 
DOD  MVA  File  by  utilizing  only  the  accession  records.  It  should  be  noted 
that  the  ''contract"  record  cannot  be  utilized  because  it  does  not  contain 
term  of  service  data. 

Development  of  the  Qualified  Military  Available  (QMA)  Inventory 

The  QMA  inventory  represents  the  estimated  number  of  people  in  the 
country  within  the  17-21  year  age  group  that  is  both  available  and 
qualified  to  enlist  in  the  military  as  of  1 July  of  a given  year.  For 
a detailed  discussion  of  the  QMA  processing  cycle  refer  to  Section  4.0 
of  this  report. 

Preparation  of  The  Basic  Recruiting  Data  File  (BRDF) 

The  third  component  of  the  Market  Information  Preparation  Module 
involves  the  creating  and  updating  of  the  Basic  Recruiting  Data  File  (BRDF). 
This  file  contains  time-series  data  describing  allocation  of  recruiting 
resources  (e.g.  recruiters  on  station,  canvassers  on  station,  etc.)  and 
socio-economic  conditions  affecting  recruiting  success  (e.g.  unemployment 
rates  and  attitudes,  toward  the  military).  The  BRDF  includes  both  monthly 
and  yearly  data  elements.  It  is  easily  expanded  to  accommodate  additional 
data  items.  The  remainder  of  this  sub-section  will  describe  the  design, 
contents  and  maintenance  of  the  BRDF. 


The  process  of  allocating  other  service  contracts  to  Army  DRC 
boundaries  has  not  been  performed  because  the  inclusion  of  the  Fips  code 
is  a recent  change  in  the  content  of  the  accession  records. 
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BRDF  Design 

The  Basic  Recruiting  Data  File  shown  in  Figure  8-2  is  a hierarchial 


file  containing  four  levels  of  data: 

• Level  1 — DRC  identifier 

• Level  2 — monthly  data 

• Level  3 — yearly  data  (non-cohort  oriented) 

• Level  4 — yearly  data  (cohort  oriented) 

For  ease  of  processing,  the  length  of  each  data  record  is  taken  to  be 
80  characters. 

Further  data  may  be  accommodated  in  the  file  by  defining  additional 
fields  in  the  current  structure  or  by  defining  additional  levels  and/or 
cards  per  level. 

BRDF  Contents 

The  data  fields  currently  defined  are  described  in  Tables  8-la 
through  8- If.  The  precise  definition  for  most  fields  is  given  elsewhere 
in  this  report  and  merely  referenced  in  the  tables. 

A group  of  data  elements  not  discussed  elsewhere  involve  other 
service  recruiters  on  station  in  an  Army  DRC.  The  methods  used  to 
estimate  these  elements  are  discussed  below. 

Air  Force:  Data  available  from  the  Air  Force  showed  the  number  of 

recruiters  assigned  to  each  of  seven  regions  by  month  for  fiscal  years 

1972-5.  The  problem  was  to  determine  the  percentage  of  the  Air  Force 

regional  staff  operating  in  each  Army  DRC,  since  every  service  has  a 

* 

different  geographic  organization.  An  Air  Training  Command  manual 

gave  the  location  of  each  Air  Force  recruiting  office  and  the  number  of 

recruiters  assigned  to  it.  Thus, it  was  possible  to  calculate  by  hand 

the  number  of  Air  Force  recruiters  whose  offices  were  located  within  the 

boundaries  of  each  Army  DRC  and  to  determine  each  DRC's  "share"  of  the 

Air  Force  regional  recruiting  staff.  Eight  Army  DRCS  straddled  Air 

Force  regional  boundaries;  therefore,  the  number  of  Air  Force  recruiters  for 

*ATC  Manual  23-3 
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Table  3-lb 


SASIC  RECRUITING  DATA  FILE  LAYOUT  — LEVEL  2 


Columns 

Data  Element 

Desc rip cion /Definition 

1-? 

DRC  Code 

3-4 

Fiscal  Year 

5-6 

Fiscal  Month 

7 

Level 

2 

8 

Card 

1 

9-13 

Cohort  Codes 

Zero-filled 

14-17 

Recruiters  on  Station 

As  described 

f 

18-21 

Canvassers  on  Station 

As  described 

22-25 

Unemployment  Rate 

As  described 

26-31 

Recruiting  Objective 

32-36 

Total  other  Service  Recruiters 

[4 

37-40 

Navy  Recruiters 

1 

41-44 

Air  Force  Recruiters 

45-48 

Marine  Recruiters 

49-50 

P0B1 

As  described  in  Table  8-lf 

51-52 

P0B2 

" 

53-54 

P0B3 

It 

55-56 

P0B4 

II 

57-58 

P0B5 

If 

59-60 

P0B6 

II 

I 

61-62 

P0B7 

undefined 

i 

63-64 

P0B8 

undefined 

65-66 

P0B9 

undefined 

67-68 

POBIO 

As  described  in  Table  8-lf 

69-74 

Advertising  D ata 

undefined 

75-80 

Advertising  Data 

undefined 
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Table  8-lc 


Columns 


14-17 


18-21 


22-25 


26-31 


32-36 


37-40 


41-44 


45-48 


49-50 


51-52 


53-54 


55-56 


57-58 


59-60 


61-62 


63-64 


65-66 


67-68 


69-74 


75-80 


BASIC  RECRUITING  DATA  FILE  LAYOUT  — LEVEL  3,  CARD  1 


s 


Data  Element 


DRC  Code 


Fiscal  Year 


Fiscal  Month 


Level 


Cohort  Codes 


Recruiters  on  Station 


Canvassers  on  Station 


Unemployment  Rate 
Recruiting  Objective 


Total  Other  Services 
Recruiters 


Navy  Recruiters 


Air  Force  Recruiters 


Marine  Recruiters 


POBIO 


Advertising  Data 
Advertising  Data 


Descrip  t ion/ Definition 


Zero-filled 


Sum  of  monthly  data 
Sum  of  monthly  data 
Average  of  monthly  data 
Sum  of  monthly  data 


Average  of  monthly  data 


undefined 


undefined 


undefined 


Average  of  monthly  data 


undefined 


undefined 


Table  8-ld 


BASIC  RECRUITING  DATA  FILE  LAYOUT  — LEVEL  3,  CARD  2 


Columns 


1-2 

3-4 

5-6 

7 

8 

9-13 

14-19 

20-25 

26-31 

32-38 

39-80 


Data  Element 


Description /Definition 


DRC  Code 
Fiscal  Year 
Fiscal  Month 
Level 
Card 

Cohort  Codes 

Total  Other  Service 
Quota 

Navy  Quota 
Air  Force  Quota 
Marine  Quota 
Undefined 


99 

3 

2 

Zero-filled 

National,  yearly  figures 

II 

II 


II 


Table  8-le 


Columns 


14-17 

18-21 

22-25 

26-30 

31-35 

36-39 

40-80 


BASIC  RECRUITING  DATA  FILE  LAYOUT  — LEVEL  4 


Dace  Element 


DRC  Code 
Fiscal  Year 
Fiscal  Month 
Level 


Cohort  Codes 
Cohort /Population 
Unfavorable  Attitude 
Unemployment  Rate 
Military  Pay  (MP) 
Civilian  Pay  (CP) 

M P/CP 
Undefined 


Description/ Definition 


As  described  in  Table  8- 
As  percent 

As  percent  of  cohort 
As  percent 


Average  for  cohort 
As  fraction 


Table  3-lf 


ARMY  POLICIES,  OPTIONS  AND  BONUSES 
Identifier  Description 


P0B1 


POB2 


P0B3 


POB4 


POB5 


POB6 


POBIO 


This  variable  is  positive  for  months 
when  recruiter  credit  was  not  given  for 
category  III,  non-high  school  graduate 
enlistees  and  0 for  all  other  months.  It 
was  initialized  at  1.0  and  scaled  down  to 
0.3  to  represent  the  gradual  withdrawal 
of  the  policy. 

This  variable  is  positive  when  a maximum 
high  school  graduate  recruiting  policy 
is  in  effect  and  0 otherwise.  The  with- 
drawal of  the  policy  is  modeled  by  a ramp. 
This  variable  is  set  to  1.0  for  July  and 
August  of  1970  during  which  time  the  Army 
restricted  two-year  enlistments.  The 
variable  is  set  to  zero  otherwise. 

This  variable  is  set  to  1.0  during  the 
time  when  the  2-year  travel  and  training 
option  was  offered  to  potential  enlistees. 

The  variable  is  set  to  zero  otherwise. 

This  variable  is  set  to  1.0  when  the  $1500 
bonus  was  in  effect  for  combat  arms  enlistees 
and  to  zero  for  all  other  periods. 

This  variable  is  set  to  1.0  when  a $2500 
bonus  was  being  paid  combat  arms  enlistees 
who  were  high  school  graduates.  The 
variable  is  set  to  zero  otherwise. 

This  variable  is  set  to  1.0  for  the  first 
quarter  of  FY76.  During  this  time,  no 
17-year  old  enlistees  were  accepted.  The 
variable  is  set  to  zero  otherwise. 
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them  was  computed  by  adding  fractions  of  two  different  recruiting 
groups.  Beckley  DRC,  for  example,  had  6.4%  of  USAF  Group  3502  recruiters 
and  1.6%  of  USAF  Group  3503  recruiters.  It  was  assumed  that  these 
percentages,  based  on  a November  1973  publication,  were  constant  over  the 
last  few  years. 

Examp le : 43  recruiters  listed  from  USAF  Recruiting  Group  3503  were 

assigned  to  the  19  offices  of  Detachment  303,  with  headquarters  in  Coral 
Gables,  Fla.  Of  these,  13  offices  with  19  recruiters  were  inside  the 
boundaries  of  Army's  Miami  DRC  and  6 offices  with  14  recruiters  were 
located  inside  the  Army's  Jacksonville  DRC.  Also  operating  inside  the 
boundaries  of  the  Jacksonville  DRC  were  26  Air  Force  recruiters  from 
detachment  302.  So  Jacksonville  DRC  had  a total  of  40  recruiters  of  the 

268  listed  for  the  region,  or  14.9%.  This  share  of  the  total  strength  of 

USAF  Group  3503  gave  Jacksonville  55  Air  Force  recruiters  in  September 
1971,  69  in  April  1975,  etc. 

Navy  and  Marine  Corps:  Mere  assumptions  about  the  geographical 

distribution  of  Navy  and  Marine  recruiters  had  to  be  made  than  were 
necessary  for  the  Air  Force  recruiters,  since  the  detailed  cross-  sectional 
data  on  the  office  address  of  each  recruiter  were  only  available  from  the 
Air  Force.  The  figures  from  the  Navy  and  Marine  Corps  were  more  detailed 
than  those  from  the  Air  Force  in  one  sense  (the  Air  Force  gave  the  staffing 

levels  for  each  of  seven  regions) , while  the  Marine  Corps  gave  levels  for 

each  of  48  recruiting  stations  and  the  Navy  for  each  of  44  recruiting 
districts.  The  problem  and  the  solutions  for  handling  the  Navy  and  Marine 

Corps  figures  were  the  same:  how  to  "map"  48  or  44  stations  into  64  Army 

* 

DRCs.  An  additional  problem  was  posed  by  the  fact  that  the  Navy  recruit- 
ing district  boundaries  changed  every  year,  there  being  37  districts  in 
FY  1971  and  44  in  FY  1975. 

In  rare  cases,  the  boundaries  of  different  service  recruiting 
districts  coincided:  USMC  RS  Los  Angeles,  for  example,  covered  the  sane 

area  as  USAREC  DRC  Los  Angeles.  In  most  cases,  though,  districts  overlapped. 

If  recruiters  were  identified  on  a county  basis  the  problem  would  be 
eased  considerably. 
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The  amount  of  overlap  was  estimated  by  dividing  the  number  of  counties 
in  the  USMC  district  which  were  inside  the  boundaries  of  a particular 
USAREC  DRC  by  the  total  number  of  counties  in  the  USMC  district.  The 
resulting  percentage  was  taken  to  be  that  DRC's  share  of  the  USMC  district's 
recruiting  staff. 

Examp le : USMC's  Denver  district  covers  107  counties  in  4 states. 

Of  these,  11  were  in  the  Army's  Sioux  Falls  DRC  and  96  (approximately  90%) 
in  the  Army's  Denver  DRC.  The  Denver  DRC  also  covers  territory  assigned 
by  the  Marine  Corps  to  their  Omaha  and  Kansas  City  districts:  the  Denver 
DRC  share  was  6%  and  10%,  respectively.  Therefore,  the  number  of  USMC 
recruiters  operating  in  the  Army's  Denver  DRC  during  FY  1975  was  obtained 
by  adding  90%  of  the  USMC  Denver  staff,  6%  of  USMC  Omaha  staff,  and  10% 
of  USMC  Kansas  City  staff  (33  + 2 + 5 * 40). 

The  important  assumption  which  had  to  be  made  was  that  Marine  and 
Navy  recruiters  were  distributed  fairly  uniformly  over  the  counties  they 
served.  In  cases  like  the  Denver  example,  it  was  felt  that  the  large 
number  of  counties  and  the  addition  of  figures  from  several  districts 
would  tend  to  even  out  the  error  introduced  by  this  necessary  assumption. 

Some  of  the  Navy  districts,  though,  had  a relatively  small  number  of 
counties,  with  suburban  and  rural  counties  located  in  a different  Army 
DRC  from  the  urban  counties.  In  these  cases,  an  assumption  that  the 
rural  or  suburban  county  in  one  Army  DRC  had  as  many  Navy  recruiters  as 
the  urban  county  in  another  DRC  would  be  a significant  source  of  error. 

To  avoid  underestimating  the  number  of  Navy  recruiters  in  certain 
metropolitan  areas  (including  Boston,  New  York,  Philadelphia,  Washington, 
Detroit,  Chicago,  and  Los  Angeles),  urban  counties  or  boroughs  were 
weighted  double. 

The  figures  which  resulted  for  the  numbers  of  Navy  and  Marine  Corps 
recruiters  in  the  Army  districts  were  felt  to  be  accurate  enough  to 
reflect  changing  force  levels  and  geographic  distribution  of  recruiters. 

It  was  not  felt  that  the  precision  implied  by  monthly  figures  was  warranted; 
therefore,  annual  average  numbers  of  recruiters  were  used  in  these  calculations. 
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BRDF  Maintenance 


The  Basic  Recruiting  Data  File  Maintenance  System  is  depicted 
schematically  in  Figure  8-3.  It  consists  of  a primary  processor 
(BRDF-FMP)  and  a supporting  processor  to  aid  in  data  preparation 
(BRDF-MUG) . The  functions  of  these  programs,  their  input  require- 
ments and  their  outputs  are  discussed  below. 

BRDF-FMP:  File  Maintenance  Program.  This  program  performs  the 

following  maintenance  functions: 

« Edit  input  data  (in  excess  of  100  edit  functions  are  provided 

0 Add  edited  data  to  file 

0 Create  summary  data  elements  (totals  and  quotients)  and 
summary  data  records  (Level  3,  card  1 records)  as  available 
data  permits 

0 Insure  consistency  of  data  by  revising  summary  data  elements 
and  summary  data  records  to  reflect  a change  to  any  component 
thereof  and  by  prohibiting  inconsistent  specifications  of  a 
summary  item 

0 Delete  individual  records  (and  automatically  any  summary  record 
constructed  therefrom) 

0 Delete  all  records  for  a given  year 

0 Provide  edit  reports  describing  any  errors  detected  in  input 
data  and  detailing  all  changes  made  to  BRDF 

Yearly  data  are  prepared  manually  in  the  appropriate  format.  Monthly  data 

are  prepared  separately  for  each  monthly  data  element  and  submitted  to 

BRDF-MUG  for  editing  and  formatting. 

BRDF-MUG:  Monthly  Update  Generator.  This  program  accepts  up  to 

12  monthly  values  for  a data  element  in  the  format  of  Figure  8-4,  edits 
the  data  and  outputs  the  data  in  the  level-2  record  format  required  by 
BRDF-FMP.  When  data  for  more  than  one  element  are  provided,  the  data  for 
all  elements  appear  on  the  same  output  record.  For  example,  suppose  that 
data  are  individually  prepared  as  per  Figure  8-4  for  recruiters  on  station, 
unemployment  and  recruiting  objectives.  When  these  three  data  sets  are 
submitted  to  BRDF-MUG,  a single  level-2  record  containing  the  appropriate 
values  for  recruiters,  unemployment  and  recruiting  objective  is  produced 
for  each  DRC  and  each  month  in  the  period  under  consideration.  BRDF-MUG 
produces  an  edit  report  which  describes  any  errors  in  the  input  data  and 
lists  the  output  records  so  as  to  facilitate  verification  of  input  to 
BRDF-FMP. 
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USAREC  Basic  Recruiting  File  Maintenance 


Market  Information  Data  Base  Module 


As  shown  in  Fig.  8-1,  there  are  four  data  files  which  comprise  the 
Market  Information  Data  Base  Module: 

• The  DOD  Contract  Summary  File  by  DRC 

• The  Historical  Accession  File  by  DRC 

• The  QMA  Summary  File  for  DRC 

• The  Basic  Recruiting  Data  File  by  DRC 

Each  of  these  files  must  be  maintained,  updated,  and  accessed  on  an  as- 
required  basis.  Each  file  with  the  exception  of  the  Historical  Accession 
File  will  be,  for  the  present,  maintained  on  magnetic  tapes.  The  Historical 
Accession  File  will  be  maintained  in  punched  card  form.  Subsequent  usage 
could  be  made  of  a generalized  data  base  package  which  would  place  these 
files  on  a computer  disk  pack. 

The  first  three  files  will,  most  likely,  be  updated  only  on  an  annual 
basis.  The  Basic  Recruiting  Data  File  will  most  likely  be  updated  on  a 
quarterly  basis.  The  process  of  accessing  these  files  for  analytic  pur- 
poses will  be  on  an  as-required  basis. 

Figures  8-5  through  8-7  present  the  format  layouts  of  the  first  three 
base  files,  respectively.  Table  8-2a  describes  the  cohort  groupings  dis- 
played in  Figs.  8-5  and  8-7.  Table  8-2b  describes  the  cohort  groupings 
displayed  in  Fig.  8-6.  The  Basic  Recruiting  Data  File  layout  was  pre- 
viously shown  in  Fig.  8-4. 

Market  Analysis  Module 

The  Market  Analyses  Module  of  MIAS  is  comprised  of  three  computer 
models: 

• A Cross-Sectional  Regression  Model 

• A Time-Series  Regression  Model 

• An  Optimal  Recruiter  Allocation  Model  (ORAM) 

Each  of  these  models  draws  upon  a combination  of  the  Market  Information 
Data  Base  files  designated  in  the  previous  section  for  their  basic  data 
inputs. 
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Table  8-2a 

COHORT  CODING  KEY  FOR  BASIC  RECRUITING  DATA  FILE, 
QMA  SUMMARY  FILE 


Code 


Group/Class 


1 

2 

1 

9 


1 

2 

3 

1 

2 

3 

4 

5 

1 

2 

3 

4 

5 

6 
7 


Sex  Group 

Male 

Female 

Age  Group 
17-21  years 
22-24  years 

Racial  Group 
White 
Black 
Other 

Educational  Categories 
High  school  graduates  - in  school 
High  school  graduates  - not  in  school 
Non-high  school  graduates  - in  school 
Non-high  school  graduates  - not  in  school 
All  including  unfit 

Mental  Categories 

I 

II 

♦ iIIA 
IIIB 
IVA 
IVB 

All  including  V 


i 
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Table  8-2b 


COHORT  CODING  KEY  FOR  CONTRACT  SUMMARY  FILE 


Code 


1 

1 

2 

1 

2 

1 

2 

1 

2 

3 


Group/Class 

Sex  Group 
Male 

Age  Group 
17-21  years 
All  others 
Racial  Group 
All  except  black 
Black 

Educational  Categories 
High  school  graduate 
Non-high  school  graduates 
Mental  Categories 
I - IIIA 
III  B 

All  others 
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The  first  two  models  are  utilized  to  conduct  analyses  to  determine 
which  factors  (i.e.,  recruiting  resources  and  socio-economic)  affect 
recruiting  success  and  the  respective  influences  of  these  factors  in 
achieving  Army  accessions.  The  cross-sectional  model  is  utilized  to 
examine  variation  in  factors  among  Army  DRCs  for  a specified  time 
period,  and  the  time  series  model  is  utilized  to  examine  variations 
in  factors  affecting  recruiting  over  time. 

The  third  model  is  utilized  to  determine  the  optimal  placement  of 
recruiters  and  canvassers  among  DRCs  based  upon  projected  socio- 
economic conditions  and  specified  accession  requirements  and,  subse- 
quently, to  determine  DRC  objective  allocations  based  upon  the  optimal 
placement  of  recruiters  and  canvassers. 

Cross-Sectional  Regression  Model 

The  cross-sectional  regression  model  is  a mathematical  regression 
technique  which  is  utilized  to  examine  the  variations  in  Army  accessions 
versus  those  factors  which  are  believed  to  influence  Army  accessions 
across  geographical  areas  of  the  country  for  a given  time  period.  For 
the  purposes  of  Army  market  analyses,  the  geographical  boundaries  have 
been  specified  to  be  the  established  Army  DRC  areas.  Typical  recruiting 
factors  to  be  analyzed  will  include: 

• Recruiters  on  production 

• Attitudes  toward  the  military /Army 

• Size  of  the  QMA  inventory 

• Unemployment  rate 

• Advertising 

• Others 

Each  of  these  factors  and  others  will  be  analyzed  to  determine  which 
factors  explain  the  variations  in  accessions  among  DRCs.  The  supply 


equation  which  best  explains  the  geographical  accession  patterns  (adhers  to 
valid  statistical  regression  requirements  and  is  historically  consistent 
with  empirical  recruiting  results)  is  selected  as  the  cross-sectional 
supply  model  to  determine  enlistment  propensities  by  DRC.  Additionally, 
the  elasticities  of  the  independent  variables  are  utilized  as  a fundamental 
input  to  the  recruiter  allocation  model. 

The  data  input  requirements  for  the  cross-sectional  analysis  are 
derived  from: 

• the  DOD  Contract  Summary  File 

• the  QMA  Summary  File 

• the  Basic  Recruiting  Data  File 

The  DOD  Contract  Summary  File  provides  historical  data  on  Army  accessions 
by  DRC  which  is  the  dependent  variable  in  establishing  the  cross-sectional 
model.  The  QMA  Summary  File  and  the  Basic  Recruiting  Data  File  provide  the 
data  for  the  independent  variables.  No  other  external  data  sources  are 
required.  If  a new  independent  variable  is  to  be  analyzed,  the  historical 
data  stream  must  be  added  to  the  Basic  Recruiting  Data  File. 

The  cross-sectional  regression  model  utilized  the  generalized  time- 
series  processor  (TSP)  computer  software  package  which  is  designed  speci- 
fically for  regression  modeling. 

Time-Series  Regression  Model 

The  time-series  regression  model  is  much  the  same  as  the  cross-sectional 
model  with  the  exception  that  variation  in  accessions  vs  recruiting  factors 
are  analyzed  over  a series  of  time  periods  rather  than  for  a given  time 
period.  The  objective  of  this  analysis  is  to  identify  the  influence  of 
recruiting  factors  over  time.  One  of  the  advantages  of  this  type  of  analysis 
is  that  shorter  time  periods  can  be  examined  (i.e.,  monthly,  quarterly  or 
semi-annual  time  periods). 

The  data  sources  for  the  time-series  model  are  the  same  as  those  for 
the  cross-sectional  model.  In  addition,  the  computer  software  package 
utilized  in  the  time-series  modeling  is  the  UCLA-developed  BIO-MED  package. 
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It  is  anticipated  that  the  time  series  model  will  be  utilized  for 
monthly  forecasting  and  that  actual  accession  data  will  be  matched 
against  forecasted  accessions  as  an  operational  capability  of  the  time 
series  model. 

Optimal  Recruiter  Placement  and  Objectives  Allocation  Model 

The  final,  and  perhaps  most  important,  computer  model  of  MIAS  is 
the  optimal  recruiter  allocation  model  (ORAM).  The  optimization  model 
consists  of  a series  of  nonlinear  equations  which  must  be  solved  to 
achieve  an  optimum  solution.  The  primary  equation  relies  on  a form  of 
the  enlistment  supply  equation  developed  from  the  cross-sectional 
regression  analysis.  The  supply  equations  are  the  objective  function 
of  the  model  with  the  projected  USAREC  recruiting  budget  being  minimized 
by  varying  recruiter  and  canvasser  levels  at  each  DRC,  subject  to 
overall  national  accession  requirements.  The  solution  results  in  the 
optimal  placement  of  recruiters  and  canvassers  among  the  DRCs  that 
yields  the  largest  number  of  accessions  at  the  least  cost.  However,  since 
other  factors  such  as  term  of  service,  occupational  specialty,  and 
cohort  specific  requirements  are  essential  considerations  in  establishing 
recruiting  objectives,  the  objective  function  described  above  is 
constrained  by  additional  equations  to  ensure  that  the  solution 
satisfies  these  additional  goals. 

DEVELOPMENT  FEASIBILITY  AND  INITIAL  USAREC  ADP  IMPACT 
General 

The  research  effort  conducted  under  the  AMRRO  Study  has  centered 
around  the  development  of  each  of  the  elements  of  MIAS.  During  the 
course  of  the  study, preliminary  methodology  reports  were  written  which 
described  the  technical  approaches  to  be  followed.  Subsequent  activity 
was  directed  toward  the  actual  data  collection,  computer  programming 
and  testing  of  each  of  the  MIAS  elements. 

MIAS  Development  Feasibility 

The  development  of  MIAS  has  shown  that  each  element  of  MIAS  is 
technically  sound  and  achievable  within  practical  limits  of  development. 
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However,  data  collecting,  data  processing,  and  computer  modeling  processes 
associated  with  recruiting  activity  can  undergo  continuing  refinement  and 
improvement  over  time;  as  new  information  is  uncovered,  new  opportunities 
arise  for  refinement  in  the  systems  design.  Since  resources  and  time  are 
not  unlimited,  it  has  been  necessary  to  place  limitations  on  the  scope  of 
GRC  activity  and  to  work  with  controlling  assumptions  to  proceed  with  the 
development  of  the  basic  MIAS  elements.  In  this  regard,  the  AMRRO  study 
team  has  placed  primary  emphasis  on  achieving  a fundamentally  sound  overall 
MIAS  capability  which  will  be  of  immediate  use  to  USAREC  management  and 
which  is  additionally  suitable  for  refinement  and  expansion  over  the  long 
term. 

At  the  initiation  of  the  AMRRO  study,  the  feasibility  of  conducting 
market  analyses  in  the  form  of  regression  and  optimal  allocation  modeling 
at  the  DRC  level  was  unknown.  GRC' s work  indicates  that  working  at  the  DRC 
level  is  both  feasible  and  practical  for  the  cross-sectional  analysis  model 
and  ORAM.  In  the  case  of  the  time-series  analysis,  given  the  availability 
of  accession  data,  regional  analysis  appears  to  be  the  lowest  practical 
level  of  analysis  at  this  time.  The  usefulness  of  conducting  analyses  at 
the  DRC  and  regional  levels  is  readily  apparent,  because  it  allows  USAREC 
management  to  examine  recruiting  activity  at  a finer  level  of  detail  than 
was  available  previously. 

In  summary,  the  feasibility  of  developing  each  of  the  MIAS  elements 
has  been  successfully  established.  It  is  further  recognized  that  the  system 
under  operational  development  will  undergo  refinement,  improvement,  and 
probably  expansion  in  succeeding  years;  however,  the  fundamental  design 
of  MIAS  will  provide  a technically  sound  capability  for  immediate  USAREC 
recruiting  management  analysis. 

Anticipated  Impact  on  USAREC  ADP  Operations 

Some  elements  of  MIAS  are  very  straightforward  and  will  have  little 
impact  on  USAREC  ADP  operations;  other  elements  are  complex  and  will  require 
implementation  planning  and  training  to  achieve  an  effective  USAREC  opera- 
tional capability. 
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Impact  of  the  Market  Information  Preparation  Module 

As  described  earlier  there  are  three  basic  data  flows  associated 
with  this  MIAS  module.  The  first  data  flow,  the  processing  associated 
with  the  creation  of  the  DOD  Contract  Summary  File  and  the  Historical 
Accession  File,  is  relatively  straightforward  and  is  comprised  of  four 
to  five  computer  programs.  The  data  processing  associated  with  this 
data  flow  involves  processing  large  tape  files  with  a substantial 
amount  of  input/output  processing;  because  of  this,  the  computer 
programs  have  been  written  in  COBOL.  It  is  anticipated  that  this 
processing  cycle  will  be  performed  quarterly  and  will  require  approxi- 
mately one  hour  of  CPU  time  and  eight  hours  of  wall  clock  processing 
time  per  cycle.  Currently,  this  flow  is  being  run  on  the 
CDC  6400  computer  at  GRC.  Only  minor  problems  are  anticipated  in 
converting  this  processing  cycle  to  a UNIVAC  1108  computer;  however, 
substantial  documentation  will  be  required  to  permit  USAREC  personnel 
to  independently  operate  this  processing  cycle. 

The  second  data  flow  is  concerned  with  the  development  of  the  QMA 
inventory  report  output  and  summary  file*  As  shown  in  Fig.  4-1,  this 
processing  flow  is  complex  and  lengthy.  In  fact.  Fig.  4-1  is  a simpli- 
fication of  the  actual  processing  requirements.  In  addition,  there 
are  three  data  sources  associated  with  this  flow  which  should  be  up- 
dated periodically: 

• The  ten-year  Census  data,  1/100,  5 & 15% 

County  Group  and  5%  State  Public  Use  Samples 
(PUS) 

• The  yearly  P-25/26  series  Census  migration  reports 

• The  yearly  "Cause  of  Death"  summary  tapes  for 
mortality  data 

Each  of  these  data  sources  is  not  complex  in  themselves;  but  a thorough 
understanding  of  the  entire  processing  cycle  will  be  required  to  keep 
the  QMA  processing  cycle  updated  and  operational.  It  is  envisioned  that 
this  processing  cycle  will  be  run  annually. 

*Note  that  the  responsibility  for  this  processing  flow  has  been 
transferred  to  MARDAC. 
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The  computer  programs  associated  with  this  cycle  are  written  in 
FORTRAN  IV  and  are  operational  on  an  IBM  360/370  computer.  However, 
some  of  the  subroutines  within  the  computer  programs  are  written  in 
machine-oriented  language  because  of  the  complexity  of  the  processing 
routines.  For  example,  in  order  to  avoid  excessive  round-off  errors 
due  to  the  IBM  word-size  characteristics,  double  word  boundaries  had 
to  be  used  in  many  cases.  In  addition,  when  processing  the  multiple 
reel  Census  tape  inputs,  a special  machine-oriented  subroutine  was 
utilized  to  reduce  excessive  processing  times. 

The  QMA  processing  cycle  is  comprised  of  approximately  thirty-six 
computer  programs  which  will  utilize  approximately  eight  hours  of  CPU 
time  and  approximately  twenty  hours  of  wall  clock  processing  time.  It 
is  believed  that  substantial  UNIVAC  1108  conversion,  personnel  orientation 
and  detailed  documentation  will  be  required  for  USAREC  representatives  to 
operate  this  data  flow  on  an  independent  basis. 

The  final  data  flow  of  the  Market  Information  Preparation  Module 
involves  the  creation  and  updating  of  the  Basic  Recruiting  Data  File. 

This  process  is  relatively  straightforward  and  is  not  extensive.  There 
are  two  computer  programs,  written  in  COBOL,  associated  with  this  flow. 

They  are  currently  operational  on  GRC's  CDC  computer,  and  conversion  to 
the  UNIVAC  1108  should  be  relatively  easy.  However,  documentation  is 
required  to  achieve  an  independent  USAREC  operational  capability. 

It  is  anticipated  that  this  processing  flow  will  be  run  on  a quarterly 
basis  with  a small  additional  amount  of  processing  required  at  the  end  of 
the  fiscal  year.  The  yearly  CPU  time  is  anticipated  to  be  approximately 
30  minutes  and  wall  clock  processing  time  about  2 hours.  The  data  inputs 
to  this  processing  flow  would  be  prepared  manually  on  punched  cards  by 
Market  Studies  and  Analysis  Office  representatives. 

Impact  of  the  Market  Information  Data  Base  Module 

This  module  consists  of  four  data  files,  three  of  which  are  maintained 
on  magnetic  tapes  and  one  of  which  is  maintained  in  punched  card  form.  The 
ADP  impact  of  this  module  will  be  essentially  nonexistent.  It  is  anticipated 
that  up  to  10  years  of  historical  data  will  be  maintained  in  each  file;  a 
given  year's  data  can  either  be  contained  on  a single  tape,  or  separate  tapes 
can  be  created  yearly. 
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If  a decision  is  made  to  implement  a generalized  data  base  software 
package  which  utilizes  computer  disk  files,  there  would  be  a more 
significant  ADP  impact;  however,  at  this  time  such  a decision  is  not 
planned. 

Impact  of  the  Market  Analysis  Module 

The  Market  Analysis  Module,  consisting  of  two  regression  models 
and  a nonlinear  programming  model,  will  have  a larger  functional  impact 
on  the  Market  Studies  and  Analysis  Office  operations  than  on  USAREC's 
ADP  operations.  The  two  regression  models  utilize  the  generalized  TSP 
and  BIO-MED  software  packages,  and  ADP  usage  will  depend  upon  the 
frequency  with  which  these  models  are  run.  It  is  believed  that  the 
TSP  and  BIO-MED  packages  are  compatible  with  the  UNIVAC  1108  computer. 

Data  inputs  to  TSP  and  BIO-MED  can  be  either  in  punched  card  or 
magnetic  tape  form.  A series  of  regression  runs  can  be  made  on  a 
single  processing  pass,  and  computer  processing  times  are  minimal.  The 
total  estimated  CPU  time  over  a given  year  is  about  three  hours  CPU  time 
and  about  five  hours  wall  clock  processing  time. 

The  major  impact  of  the  regression  modeling  would  be  experienced 
by  the  functional,  operations  research  representatives  who  would  be 
required  to  understand  regression  and  statistics  theory  and  to  prepare 
the  appropriate  data  inputs  for  the  software  packages.  These  data 
inputs  are  accessed  from  the  data  files  maintained  in  the  Market 
Information  Data  Base  Module. 

The  ORAM  is  more  complex  than  the  regression  models.  This  model 
is  currently  operational  on  GRC's  CDC  computer,  and  has  been  specially 
programmed  for  USAREC's  application.  The  computer  language  utilized  is 
FORTRAN  IV;  conversion  requirements  to  make  this  model  operational  on 
the  UMTVAC  1108  are  unknown  at  this  time.  However,  it  is  believed  that 
conversion  problems  would  only  be  minimal. 

Substantial  documentation  will  be  required  for  USAREC  representatives 
to  modify  and/or  operate  this  NLP  model.  All  data  inputs  are  currently 
in  punched  card  form  and  computer  processing  times  are  not  extensive. 


E-31 


Ic  is  estimated  that  on  a yearly  basis  two  hours  of  CPU  time  and  four 
hours  of  wall  clock  processing  time  will  be  sufficient  to  operate  the 
model  on  an  as-required  basis. 

As  in  the  case  of  the  regression  models,  the  ORAM  will  require 
that  operational  research  representatives  on  the  functional  side  have 
a thorough  understanding  of  the  model  in  order  to  prepare  the  proper 
data  inputs  and  to  interpret  output  results. 

IMPLEMENTATION  SUMMARY 

The  Market  Information  and  Analysis  System  (MIAS)  developed  by 
GRC  is  judged  both  feasible  and  practical.  It  is  acknowledged  that 
refinements,  improvements  and/or  expansion  in  selected  areas  could  be 
expected  over  time.  However,  the  fundamental  MIAS  design  will  prove 
to  be  a very  useful  and  viable  system  to  recruiting  management. 

Impacts  on  USAREC  operations  will  be  two-fold:  ADP  impact  and 

functional  operations  research  impact.  The  ADP  impact  is  summarized 
on  Table  8-3;  it  is  estimated  that  approximately  17.5  CPU  hours  and  approxi- 
mately 63  wall  clock  processing  hours  yearly  would  be  required  to  operate 
MIAS.  Conversion  problems  in  utilizing  the  UNIVAC  1108  computer  require- 
ments will  be  substantial  for  the  DOD  Contract  File  processing  flow,  the 
QMA  processing  flow,  and  the  NLP  optimization  model. 

The  functional  operations  research  impact  will  result  from  two  areas 
of  MIAS:  basic  data  collection,  and  the  operation  and  understanding  of  the 

computer  models. 

In  summary,  Implementation  of  MIAS  will  require  an  understanding  of 
unique  elements  of  the  system  by  both  the  ADP  and  functional  organizations 
at  USAREC.  The  ADP  impact  of  MIAS  on  USAREC  operations  should  be  relatively 
small  in  comparison  with  other  information  system  requirements;  however,  it 

♦ 

is  believed  that  initial  implementation  will  require  conversion  activity, 
substantial  documentation,  and  user  orientation.  In  order  to  achieve  an 
annual  operational  capability  it  is  recommended  that  a small,  identified 
group  of  both  ADP  and  functional  representatives  be  established  to  maintain 
and  operate  the  system  on  a continuing  basis. 
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PRELIMINARY  IMPLEMENTATION  PLAN 

Since  all  of  the  MIAS  has  essentially  been  developed,  each  individual 
portion  or  all  portions  could  be  implemented  either  concurrently  or  se- 
quentially. However,  to  permit  early  utilization  of  the  Market  Analysis 
Module,  it  is  recommended  that  the  existing  Market  Information  Data  Base 
Module  be  implemented  first;  the  Market  Analysis  Module,  second;  and  the 
Market  Information  Preparation  Module,  last.  The  sequential  schedule  will 
permit  early  utilization  of  the  analytical  models  with  existing  data  files 
at  minimal  impact  on  other  USAREC  ADP  operations.  Such  an  approach  should 
be  advantageous  to  the  ADP  operation's  representatives,  because  it  is  anti- 
cipated that  they  will  be  faced  with  a difficult  schedule  in  converting 
existing  operational  activities  to  the  newly  installed  UNIVAC  1108  system. 

Given  the  implementation  of  the  existing  data  base  and  analytical 

models,  the  three  basic  data  streams  associated  with  the  Market  Information 

Module  can  be  implemented  concurrently.  In  this  regard,  responsibility  for 

the  QMA  data  processing  sequence  has  been  transferred  to  MARDAC  (OSD)  and 

implementation  of  this  sequence  will  be  independent  of  the  other  elements 

of  MIAS.  To  facilitate  implementation  and  operations,  it  is  recommended 
t 

that  USAREC  identify  specific  individuals  to  have  cognizant  responsibility 
for  each  element  of  the  MIAS  operational  and  maintenance  activity. 

Figure  8-8  presents  a suggested  implementation  schedule  for  planning  purposes. 
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Fig.  8-8  — Preliminary  Milestone  Chart  for  MIAS  Implementation 
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1 

□ 

658 

White  Sands  Missile  Ronge 

1 

f—1 

FI 

US  ARMY  FORCES  COMMAND 

Hq,  F0RSC0M 

\ 

c 

F2 

First  US  Army 

1 

LJ 

F3 

Fifth  US  Army 

1 

□ 

F 4 

Si'Xth  US  Army 

1 

— 

F5 

Eighth  US  Army 

1 

Required  IAW  AR  5-7. 
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Address 

code 

Agency 

Number  of 
copies 

□ 

G1 

UNIFIED  COMMANDS 

Commander  in  Chief,  Alaska  (CINCAL) 

i 

□ 

G2 

Commander  in  Chief,  Pacific  (CINCPAC) 

i 

□ 

G3 

Commander  in  Chief,  Europe  (CINCEUCOM) 

i 

□ 

G4 

Commander  in  Chief,  US  Readiness  Command 

i 

□ 

GS 

Commander  in  Chief,  Atlantic  (CINCLANT) 

i 

□ 

G6 

Commander  in  Chief,  Continental  Air  Defense  Command  (CINCONAD) 

i 

□ 

G7 

Commander  in  Chief,  US  Southern  Command  (USCINCSO) 

i 

□ 

HI 

US  ARMY  COMMANDS/AGENCIES  (CONUS) 
US  Army  Air  Defense  Command 

i 

□ 

H2 

US  Army  Communications  Command 

i 

□ 

H3 

US  Army  Computer  Systems  Command 

i 

□ 

H4 

US  Army  Recruiting  Command 

i 

□ 

H5 

US  Army  Safeguard  Command 

i 

□ 

H6 

US  Army  Health  Services  Command 

i 

□ 

H7 

US  Army  Intelligence  Command 

i 

n 

H8 

US  Army  Combat  Developments  Experimentation  Command 

i 

□ 

HI5 

US  Army  Security  Agency 

i 

□ 

H16 

US  Army  Logistics  Evaluation  Agency 

i 

□ 

HI  7 

US  Army  Concepts  Analysis  Agency 

i 

□ 

H18 

US  Army  Management  System  Support  Agency 

i 

□ 

H19 

US  Army  Operational  Test  and  Evaluation  Agency 

i 

□ 

H20 

US  Army  Missile  Intelligence  Agency 

i 

□ 

H21 

US  Army  Nuclear  Agency 

i 

□ 

H22 

US  Army  C*E  Engineering  and  Installation  Agency 

i 

□ 

H23 

US  Army  Communication  Systems  Agency 

i 

□ 

J1 

MISCELLANEOUS  ARMY  (CONUS) 

US  Army  Intelligence  Threat  Analysis  Detachment 

i 

□ 

J2 

Engineer  Studies  Group 

i 

□ 

J3 

US  Army  Research  Institute 

i 

□ 

J4 

US  Army  Military  Personnel  Center 

i 

□ 

K1 

US  AIR  FORCE 

Hq,  US  Air  Force  (SAMI) 

□ 

K2 

Hq,  US  Air  Force  Systems  Command 

i 

□ 

K3 

Assistant  Chief  of  Staff  for  Intelligence 

i 

□ 

K4 

Foreign  Technology  Division  (AFSC) 

i 

□ 

K5 

Air  University  Library 

i 

□ 

LI 

US  NAVY 

Chief  of  Naval  Operations,  OP-96 

i 

□ 

L2 

Chief  of  Naval  Operations,  OP03EG-CNO 

i 

c 

L3 

Naval  Intelligence  Command 

t 

□ 

L4 

Naval  Intelligence  Support  Center 

i 

□ 

L5 

Naval  War  College 

i 

□ 

L6 

Naval  Postgraduate  School 

i 

□ 

Ml 

US  MARINE  CORPS 

Hq,  US  Marine  Corps,  G-2 

i 

□ 

M2 

Hq,  US  Marine  Corps,  G-3 

i 

□ 

M3 

Marine  Corps  Development  and  Education  Command 

i 

□ 

N1 

US  GOVERNMENT  AGENCIES 

Director,  Central  Intelligency  Agency 

i 

□ 

N2 

Central  Intelligence  Agency  (Publications) 

i 

□ 

N3 

Central  Intelligence  Agency,  Director  for  S&T 

i 

w 

1 

r 

■ 
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X 

1 

Address 

Number  of 

I 

code 

Agency 

copies 

1 

□ 

N5 

Director,  Office  of  Science  & Technology,  Executive  Office  Bldg 

i 

□ 

N6 

Department  of  State,  Foreign  Affairs  Research  & Doc  Center 

i 

DEFENSE  CONTRACTORS 

a 

01 

Center  for  Naval  Analyses 

i 

a 

02 

Human  Resources  Research  Organi  zation,  Inc. 

i 

a 

03 

institute  for  Defense  Analyses 

i 

a 

04 

RAND  Corporation 

i 

□ 

05 

American  Institute  for  Research 

i 

□ 

06 

Stanford  Research  Institute 

i 

^3 


Headquarters,  Department  of  the  Army 

Army  Research  Institute  1 

US  Navy 

Naval  Personnel  Research  and  Development  Center  (Govindan)  1 


